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We review the so called A-scaled Nambu-Jona-Lasinio (NJL) model both in the vacuum 
and at finite temperature and density in the case of 3 flavors (it, d, s). Starting from QCD 
and integrating out the gluons, we show how to reintroduce some of their effects in order 
to take into account the trace (or scale) anomaly of QCD. The axial Ua{1) anomaly is 
also present in the model in order to discuss the rj' particle. In the present work we treat 
the pseudoscalar and scalar sectors, pointing out the defects related to the omission of 
the vector particles. We use our model to study the restoration of chiral symmetry as a 
function of temperature, the thermodynamical functions and the equation of state, and 
we survey the mixing between the scalar isoscalar particles and the glueball. 



Contents 



> 
o 

o 

o 

00 



X 



Introduction 



A Symmetries 



B From QCD to GCM 



3 

3 
6 



C Fierz identities, Feynman-like gauge, covariant gauge 



D 



Fierz identities in color space) 9 

Fierz identities in flavor spa ce] 10 



Fierz identities in Dirac space] 11 

13 

Th e 3-flavor scaled NJL models and anomalies 14 

BosonizationI 15 

A and B scaled models 18 



E Parameters 24 



Gap equations (or Schwinger-Dyson equations^ 18 

Equations for meson masses (or Bethe-Salpeter equations)] 22 



Drawbacks of the models 27 



II The A-scaled NJL model at finite temperature and density 



A Thermodynamic! 



1 Pressure, energy density, entropy density, bag constant 



2 Bag constant 



3 High temperature zero density limit (T > 



29 

29 
29 
32 
33 



Low temperature zero density limit 34 



B Results: condensates and thermodynamical functions at flnite temperature and density 34 



I 



2 Critical surface T, 



Condensatetj 34 

39 



Thermodynamics 45 



C Anomalous decay ttq 77 



1 Deflnition of the model above the threshold tt 



Comparison with lattice QCD] 50 

52 

53 



2 Disintegration ttq 77 



qq 



A Hybrids: from unphysical pure nn, ss and glueball fields to physical particle; 



B 



56 

60 

61 

Mixing angles and coupling constants as an eigenvector problem 63 



III Mixing betvyeen scalar isoscalar mesons and the glueball] 



C Two-photon decay of scalar isoscalars 65 



IV Conclusions 



68 



APPilNDIXES 



69 



A Convention^ 



B 



E 



Color quark current-current interaction 



Propagator matrices 



D 


Thermodynamics 








1 


Pressure . . . 










a High temperature zero density expansion 




b Low temperature zero density expansion 






c Finite density, zero temperature! .... 






2 


Energy density! 






a Higli temperature zero density expansior 







b Low temperature zero density expansion 






[; Finite density, zero temperature! .... 






3 


Entropy density 






a Higli temperature zero density expansion 




b Low temperature zero density expansion 




z Finite density, zero temperature 



69 
69 
72 
74 



Masses and on-shell coupling constants for mixed particles 

1 Generalities 



79 
80 

80 

First order Taylor expansion, as a function of q , around the respective zeroes of the 
eigenvalues 81 



Bibliography 



84 



Pictures 4-7, 9-17, 20-23 are reprinted from 



• Nuclear Physics A582 , M. Jaminon, B. Van den Bossche, " SU{3) scaled eifective Lagrangians for 
a hot and strange system", 517-567, Copyright 1995 

• Nuclear Physics A592, D. Blaschke, M. Jaminon, Yu. L. Kalinovsky, P. Petrow, S. Schmidt, B. Van 
den Bossche, "Anomalous pion decay in effective QCD at finite temperature", 561-580, Copyright 
1995 

• Nuclear Physics A598, J. Cugnon, M. Jaminon, B. Van den Bossche, "Phase transition and ther- 
modynamics of a hot and dense system in a scaled NJL model", 515-538, Copyright 1996 



with permission from Elsevier Science. 



2 



I. INTRODUCTION 



Quantum chromodynamics (QCD), the theory of strong interaction, is weU suited to perform calculations 
at high energy-momentum. The non-linearity in the gluon sector (three- and four-gluon vertices) together 
with the high value of the coupling constant make it difficult, if not impossible, to apply standard 
calculational techniques for the description of the hadronic sector. One has then to choose between two 
strategies: either to put QCD on a lattice or to search for simplified models believed to mimic QCD in a 
given range of energy-momentum. Needless to say, due to the structure of QCD (high number of degrees 
of freedom (d.o.f.)) and problems related to the treatment of fermions on the lattice (fermion doubling, 
fermion determinant), lattice gauge theory is restricted to a very small number of space and time slices. 
Although improved discrete actions can lead to the hope that these small lattices already give stable 
outputs, it can be advantageous to turn to the second strategy. It is not, of course, solving QCD but it 
can shed some light on the way the true theory behaves. 

It is not our purpose here to describe the different models which exist. We shall mainly focus on an 
extension of the Nambu Jona-Lasinio (NJL) model, while keeping in mind that it is a simplification 
of the Global Color Model (GCM). The NJL model was invented by analogy with superconductivity. 
The physical phenomenon underlying the BCS theory is the creation of an energy gap, resulting from 
interactions of electrons with phonons, near the Fermi surface in the spectrum of one-fermion excitations. 
However, the gap is not sufficient to realize superconductivity whose essence is related to the appearance 
of the gap through the dynamical breaking of U{1) gauge symmetry. 

In the absence of a current quark mass term m, QCD possesses an additional symmetry named the chiral 
symmetry which is related to the fact that the two inequivalent representations of the Lorentz group 
for spinors ((1/2,0) and (0,1/2)) transform independently (if m — 0). As soon as this symmetry is 
dynamically broken, one might expect that the fermions will acquire (dynamical) masses. For a review 
of the dynamical breaking of chiral symmetry in quantum field theory, see |jl| . 

A. Symmetries 

In this section, we discuss the symmetries connected to the QCD Lagrangian which are of concern for 
the construction of effective theories modeling it. 

The basic theorem for the examination of symmetries is Noether theorem. It states that for any invariance 
of the action under a continuous transformation of the fields, there exists a classical charge Q which is 
time-independent and is associated with a conserved current, d^J'^ = 0. 

A deviation from this theorem implies that either the symmetry is only an approximate one (for example 
the chiral symmetry for light quark masses) or that quantum corrections have broken it: we enter the 
world of anomalies. Moreover, there still exists the possibility that the symmetry is spontaneously broken 
(this is the case of the chiral symmetry). 

Although gauge invariance is of prime importance in QCD, because it is the basis to justify the QCD 
Lagrangian, it is of little interest in this work since we shall work within the frame of the NJL model 
which is not locally gauge invariant. (It is however globally color invariant.) The relevant symmetries in 
the present context are shown below, where we indicate the corresponding transformations of the fields, 
and the status of the conservation^. We restrict ourselves to three fiavors of quarks u, d, s and we take 
the convention to denote by m the current quark mass matrix diag(TO„, m^, m^), and by q the vector 
representing quarks in the flavor space q =diag((7„, qd,qs)- 

• Global gauge symmetry: as stated above, NJL is not locally gauge invariant but only globally 
invariant. The color enters only through the number Nc of each quark flavor. In the GCM model, 
the important features of local gauge invariance of the color group SU(3)c such as confinement and 
asymptotic freedom can be parametrized in the form of an effective two-point gluon propagator 
(see below). However, this is not the case for the NJL model where an unphysical meson to quark- 
antiquark pair threshold is a direct consequence of the lack of confinement. 



In addition to these symmetries, QCD is also invariant under Lorentz transformations, P, C, T and G— parity, 
which are also symmetries of the NJL action. 
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Scale symmetry: 

g - x'/\ 

is exact in the limit of vanishing current quark mass m = 0. However, this symmetry is broken by 
quantum effects (see below). 

• Vector U{l)v symmetry or Beiryonic number conservation: q exp{ia)q. 

• Quark number conservation: Qi — > exp(ia,)gi {i = u,d,s): each flavor has its own conserved 
number. 

• Axial U{1)a symmetry: q exp{ij5a)q is exact in the limit of vanishing current quark mass 
TO = 0. This symmetry is broken by quantum effects, which explains why it is not seen in the 
spectrum of physical states. 



Isospin symmetry: 

^ exp^iT"^ aa)q , a = 1, 3. 



u 
d 

Contrary to axial symmetry, isospin symmetry is exact in the limit where light quark masses are 
equal to„ = nid {Ta are the Pauli matrices). 

• Vector SU{^)v symmetry: 

exp(zY"y)9 ' " = l>-,8. 

The isospin symmetry can be generalized to the three-flavor case and is exact if to„ = to^ = (Ao 
arc the S'C/(3) flavor GcU-Mann matrices). This approximate symmetry explains why the hadrons 

are ordered into multiplcts. 

• Axial SU(2)a and S'[/(3)a symmetries: 

exp(i75YaA)9 , a=l,...,8 (5'{7(2)^ : A„ ^ t„). 

The quark flavors can also be mixed in the axial case. This symmetry is exact as long as to„ = 
nid = rris = 0. This approximate symmetry is not seen in the spectrum. Axial transformations 
alter the parity that is associated with a state. A manifestation of SU{2,3)a in nature would 
require that each isospin (or SU{3)v) multiplet be accompanied by a mirror multiplet of opposite 
parity. In the same way, since we do not observe opposite parity partners to all hadrons, the U{1)a 
symmetry cannot be realized directly by QCD. While the axial SU{3)a symmetry is realized in the 
Goldstone mode through the dynamical breaking of chiral symmetry, U{1)a is never realized, being 
completely broken by quantum effects (the anomaly) . 

Chiral symmetry: vector and axial SU{3) symmetries can be combined to realize transformations 

1 ~F Ts 

on the left and right parts of the quarks (^r = — - — q): 

L Z 

5[/(3)y SU{^2>)a ^ SU{3)l SU{3)r, 

where 

SU{3)l =^qL^ exp(iyOgL, SU{3)r gjj ^ exp(iya?j)to. 

Under chiral symmetry, left-handed and right-handed quarks transform independently. This sym- 
metry is broken by the quark mass matrix m. Apart from this explicit breaking, chiral symmetry 
is also spontaneously broken down to SU{3)v- 



4 



Spontaneous symmetry breaking and anomalies 



As mentioned above, a symmetry can be manifested in several ways. 



It may remain exact. This is the case of the electromagnetic gauge symmetry (except for supercon- 
ductivity) and the color gauge symmetry of QCD^. 

It may be explicitly broken. This is the case of the isospin symmetry if rriu ^ md or the axial SU{2) 
or SU{3) symmetry if m 7^ 0. 

It may be hidden. It is an invariance of the action but not of the ground state: the symmetry is 
not seen in the spectrum of physical states. Two types of mechanisms are possible. 

1. The symmetry can be spontaneously broken. It is the case of the SU{2)l symmetry in elec- 
troweak interactions. The breaking requires the presence of scalar fields (Higgs fields in elec- 
troweak interactions) which give rise to vacuum expectation values. The order parameter 
associated with the breaking is based on an elementary operator. In the case of the Higgs 
mechanism, the order parameter is the vacuum average of the scalar field. 

2. The symmetry can be dynamically broken: it does not require any scalar field. This is the 
case of the chiral SU (2, 3)l (g) SU{2, 3)^ symmetry in QCD, or the U{1) electromagnetic gauge 
symmetry in superconductivity. The order parameter is given by the vacuum average of a 
composite (and not elementary) operator. In the case of dynamical chiral symmetry breaking, 
the order parameter is the quark condensate^ < OlgqjO >. 

In the following we shall use the term "spontaneous symmetry breaking" to describe both cases 
of hidden symmetry, making the distinction when appropriate. To the notion of hidden symmetry 
can be attached the Goldstone theorem: if a theory has a continuous symmetry of the Lagrangian 
which is not a symmetry of the vacuurr^, there must exist one or more massless bosons (Goldstone 
bosons) . This means that spontaneous or dynamical breaking of a continuous symmetry will entail 
massless particles in the spectrum. The Goldstone theorem, associated to the dynamical breaking 
of chiral symmetry, explains the small pseudoscalar nonet mass (ttq, tt±, Kq, Kq, rj, 77'), except for the 
rj' . In fact, in the limit of vanishing current quark mass, both pions, kaons and rj have a vanishing 
mass in accordance with the Goldstone theorem. Their nonzero value is just the reflection of the 
explicit breaking of chiral symmetry by current quark mass. The fact that the strange quark mass 
is heavier than the up and down quark mass implies a greater corresponding expli cit br eaking and 



thus a heavier mass for the kaons (and rj' without the axial anomaly, see section [ D 4 ) compared 
to the pions. The breaking of chiral symmetry can be summarized by 

SU{3)l(E> SU{3)r ^ SU{3)v 



• The symmetry may have an anomaly. An anomaly is a symmetry of the action (classical level) 
which can be broken by quantum effects: the Noether theorem is no longer valid. In QCD, this 
is for example the case of the axial U{1)a and scale symmetries: even if the quark masses are 
identically equal to zero, the corresponding current does not vanish. For the axial U{1) symmetry, 
we get {F, F are the SU{3) field strength tensor and its dual, respectively, /3qcd is the Callan- 
Symanzik /3-function of QGD and 7^ is the mass anomalous dimension) 

d'imY'l^ - 2zg-75m°|^g + ^^F^^Pr, (1) 



^See however the works [^-^ and where the authors, taking care of both the quark and diquark d.o.f., show 
that the global color symmetry might be dynamically broken by the formation of a < > condensate of quark 
Cooper pairs. 

^ Other order parameters are possible. The quark condensate is the simplest one. However, it should be stressed 
that although < 0|gg|0 >^ implies dynamical chiral symmetry breaking, the reverse is not true. 

''This is the Goldstone mode of the symmetry, by opposition to the Wigner mode. Goldstone mode can only 
occur in the case of an infinite number of d.o.f. 
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while the scale anomaly reflects itself in the form 

d,j^ = e^^ = {i + 7™) E ^^^^-^^ + ^K^Pr- (2) 

,:=l -9 

An other example of anomalyn concerns the decay of neutral pions into two photons ttq — > 77. 
The isovector axial current = "07^75?/ — d'y^j^d is not conserved in the quantum world: its 

divergence leads to a term FF, where F is now the electromagnetic field strength tensor. The 
strong axial anomaly (|l|) is believed to give its high mass to the ry' particle compared to the other 
members of the pseudoscalar nonet. 

With this small introduction on symmetries, dynamical breaking and anomalies, we are now in a position 
to understand the basis of the NJL model which is described below as a limiting case of the global color 
model. 



B. From QCD to GCM 

Our purpose is to discuss the NJL model as an approximation of QCD and its extensions. Our starting 
point is the global color model (GCM) of which NJL is a particular case. 

QCD is the theory of quarks and gluons. However, gluons are often neglected in favor of an effective 
interaction between quarks. This is the case for non-relativistic models such as the quark model as well as 
for relativistic models such as the NJL model. For models incorporating the dynamical breaking of chiral 
symmetry, one argues that the interaction between quarks and antiquarks, coming from complicated 
processes of gluon exchange, is attractive and leads to an effective interaction between quarks. This 
interaction is responsible for quark-antiquark condensation in the vacuum, when the interaction exceeds 
a critical strength. Although still a symmetry of the Lagrangian (if m = 0), the condensation is responsible 
for the fact that the vacuum is no longer symmetric under SU (3)^ transformations. This leads to massless 
Goldstone bosons as well as a dressing of quarks by qq pairs. 

The shape of the effective interaction is important to model QCD and its choice is a competition between 
mathematical simplicity (still retaining some important aspects of QCD such as dynamical symmetry 
breaking and, more generally, the relevant global symmetries) and conceptually more appealing ap- 
proaches, at the price of loosing simplicity. While NJL belongs to the former, the latter corresponds to 
the more general GCM-type model. 

In this section, we summarize the main steps which enable going from QCD to GCM. 
The QCD 0| Lagrangian is given by: 

Cqcd = 9(*7"9m - - + gq^A^q, (3) 

where q is the quark field in flavor and color spaces (in the fundamental representation) and g is the 
coupling constant. The quantity is a shortened notation for the eight gluon fields = yl°A"/2 in the 
adjoint representation, where A" (a = 1...8) acts in the color space. The second term is the Yang-Mills 
Lagrangian, constructed from the gluon field strength: 

F^. - d^^A: - d^Al + gr'^A^Al. (4) 

The last term of eq. ^ gives an interaction term between gluons (three- and four-gluon couplings) and 
is due to the non-Abelian structure of the theory (the totally antisymmetrical coefficients f"^'"^ are the 
structure constants of SU{3)). The last term of the Lagrangian (||) gives the coupling between the gluon 
fields A1^ and the color quark currents 

Ja^ = g Y7^g- (5) 



^The reference Q is a good introduction on anomalies in the standard model. Its chapter VII-3 deals with purely 
hadronic processes through the study of the Wess-Zumino-Witten anomaly action. 
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In quantum theory, a given Lagrangian has to be suppHed with a quantification procedure which is 
chosen here to be the Feynman prescription: solving the quantum theory is equivalent to search for 
the vacuum- vacuum transition amplitude defined through the following functional integral in Minkowski 
space: 



Zqcd = j VqVq J PA^lexp l^i J cI^xCqcd j ■ (6) 

The integration over gauge fields has to be carried out only over gauge inequivalent orbits, which is 
performed by fixing the gauge and using the Faddeev-Popov trick ||^. In the following, this is supposed 
to be included in the gluon measure VAI^. Several procedures have been investigated to manage the 
functional integral (H). For example, Alkofer and Reinhardt [p|-p^ use the field strength method to 
obtain an effective action for an auxiliary field, describing the non-perturbative vacuum already at tree 
level. Making a quadratic expansion around this auxiliary field, these authors are able to obtain a color 
current-current interaction, which is local in the limit of small momenta, while being identical to one-gluon 
exchange at high momenta. 

Another, more usual, approach is possible and described for example in ||ll|,|l^ . We shall however follow 
the notation of The integration over the gauge fields in (||) is not Gaussian. It cannot be evaluated 
exactly^!]. Equation (^) can be rewritten in the form 

Zqcd = J VqVqexp (^t J d^xqii^^d^ - m'')q + tW[j]^ , (7) 

where 

tW[j] =lnj VAl exp (^-'- j d'x{Fp^ + ig j d'xAl^^^ . (8) 
The idea is to expand (|^) in powers of quark currents 

iW[j] = iW[Q]+g j W^^\xi)lj^{xi)d^xi 

+ W(^\xuX2)l\%\j^l{x,)j^l{x^)d^x^d^x^ + ... 
+ 9l f win) ^ ^ ^ a;„);:i::-: j,^^^ {x^)... jH: . . . 



+ ... (9) 

The coefficients . ■ ■ Xn)l\::!tl are the connected n-point functions of the gluons without quark- 

loop contributions. The second term is 

^-'W^'Hx^,X2r,\% = (Al^li^MZi^^)) - (^mU^i)> (Kli^^)) , (10) 

which readily shows that W''^\xi, X2) is the gluon correlation function^. The GCM model consists in 
keeping only W^'^\ modeling the effect of the suppression of the higher n-point Green functions {n > 3) 
through the use of a particular shape for W^*-^-* [^15-1^. This shape can be chosen to reproduce 



important features of QCD such as asymptotic freedom and confinement |18|-p0|. However it can be 
shown that the inclusion of non-canonical quark-gauge boson vertices is still necessary. 



^In the field strength formalism, Alkofer and Reinhardt Q| get a Gaussian form for the auxiliary fields they 
introduce. This can be integrated exactly, leading to an effective functional for the auxiliary fields 
^An exact calculation can also be performed for the approximations QCD2 and Abelian QCD4, jisljl 



*We have defined 



_ / PA°Oexp(-»/4 / d^xjF^,,)^) 
J VAl exp (-1/4 / d^xiF^^y) ■ 
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Each function W*-""-* is separately invariant under Lorentz and global gauge color transformations. How- 
ever, only the total sum (y) is invariant under local gauge transformations: its truncation to a given order 
is then gauge dependent. 

The model is finally described by the Lagrangian 

Hgcm = q{il^d^ - m)q + ^jii{x)W^^\x,y);tf,{y), (11) 

where W'-^^x, y^Jl has the interpretation of an effective gluon propagator. 

This model is very appealing in the sense that, while mathematically tractable, it 

• gives the correct high momentum limit (through the asymptotic freedom form of the propagator) ; 

• leads to quark confinement (although giving no insight on the mechanism of color confinement, 
choosing a form of the gluon propagator so that the quark propagator has no pole on the real axis 
is a sufficient ~ albeit not necessary - condition to ensure the absence of propagation of colored free 
quarks fl^); 

• is renormalizable (the non-locality of the interaction makes a natural cut-off on high momenta 
implied in loops and makes the theory renormalizable. We shall however not discuss the renormal- 
ization procedure further on in this paper). 

However, the model has the drawbacks associated to its qualities, apart from the lack of local color 
invariance: 

• the first limitation is of mathematical nature: although tractable, the model is quite heavy to 
handle; 

• a second limitation is of conceptual nature. To ensure the convergence of calculations - see the 
i factor in the exponential of (||) - the model has to be rewritten in Euclidean space|^ which has 
none of the problems due to the indefinite norm inherent to Minkowski space. However, since the 
transcription of non-pcrturbative equations from one space to the other is not equivalent (except for 
simple models) to the analytic continuation of the solutions, the field theory has to be considered 
as defined in Euclidean space. This means that a solution of eq. ( pT| ) written in Minkowski space 
is not a solution of eq. (|lj) written in Euclidean space. Indeed, the essential singularity at infinity 
which is encountered when the quark propagator is an entire function implies that the Wick rotation 
cannot be used to justify a change of space (Euclidean — Minkowski) by just transcripting the form 
of the equations pH ; 

• one should not forget that there is no reason to expect that the effective interaction should have 
the form of a gluon exchange interaction. The "derivation" we have given is at best an attractive 
intellectual game (which, however, gives a beautiful frame to the model). 

C. Fierz identities, Feynman-like gauge, covariant gauge 

As we saw in the previous section, removing the full set of gluon Green functions except for the 2-point 
one implies that i) we loose the local gauge invariance of the theory and ii) we have to model the 2-point 
Green function. The modeling depends not only on space-time but also on internal indices. Without loss 
of generality, the gluon propagator can be written]^ 

Df,{x, y) ^ g^W'^^\x, y)f, = SabD^.{x, y). (12) 

In this section, the summation over repeated indices is implicit. We work in Minkowski space and the 
notational conventions are from Itzykson and Zuber ||2^ , see appendix ^ In particular, we have the 



Note also that results from lattice QCD and from studies through Schwinger-Dyson equations are obtained in 
Euclidean space. 

'^"Xranslational invariance implies D{x,y) = D{x — y). To derive general relations, we shall however keep the 
first form up to the end of the derivation of the Fierz identities. 
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metric 5'''' = diag(l, —1, —1, —1), 7^ = 75 = i7°7^7^7'^, and the totally antisymmetric Levi-Civita tensor 
is given by (with e^^p„ = -e^"'"") 

+1 if {/i, /5, cr} is an even permutation of {0, 1, 2, 3}, 

— 1 if it is an odd permutation, (13) 
otherwise. 



Fierz transformations are useful identities allowing to rewrite the color-color interaction (11) in terms of 
physically observable quantum numbers. This is a standard procedure when dealing with mesons and 
for a 7'' (g) 7^ structure such as in the Feynman gauge. This procedure has also been recently applied 
to the case of diquarks, taking them as building blocks for baryons. However, there has been little 
effort to derive these identities in other gauges. This is not important in the Schwinger-Dyson (SD) 
and Bethe-Salpeter approach, where Fierz identities are not needed. (For a review of the SD approach 
and its application to hadronic physics, see [^3| .) However, the choice of the gauge is important in the 
GCM model, where one introduces the physical fields in terms of the corresponding quantum numbers 
as early as in the GCM action coming in the path integral formalism. In fact, this choice leads the 
appearance and the coupling between these physical d.o.f. For example, there is no tensor particle in a 
Feynman-like gauge. (We denote by Feynman-like gauge |^l[| a gauge whose Dirac structure is the same 
as the Feynman gauge for the perturbative gluon propagator: D^^(x,y) = g^yD{x,y). It is clear that, 
in non-perturbative studies, the Feynman-like gauge is not a conceptually good choice, since it does not 
yield the Slavnov- Taylor identity ki^^k^D^'^ — kf^ki/D'^'^ (-Dq" being the free gluon propagator) which is 
valid in arbitrary covariant gauges.) However this "could be" gauge is often used in GCM-type studies 
because it simplifies considerably the use of Fierz identities. Although the NJL model is also based on 
the Feynman-like gauge, we treat below the case of an arbitrary covariant gauge. 

In the following, we shall first focus on the color, then the flavor Fierz transformations (both for mesons 
and diquarks) before turning to the spin structure. We shall then summarize our results. 



1. Fierz identities in color space 
qq sector 

As long as one is not interested in the diquark sector, one can use the standard identity 

whose derivation starts from the obvious decomposition of a Hermitian matrix A = AqTL + AcX'^. Taking 
the traces tiA and tiX'^'A allows to obtain Aq and A^. 

A=i-tr(A)l + itr(AA^)A^ (15) 

We then have 

Aaf} = —AssSalS + -{A^sXs.y)K)3: (16) 
Aa0 = AjS^-jaSsp — A^si — SapSjS + '^Xg^X'^p): (17) 

i.e. 

1 c c 1 

l^^SjK/S = <^7"'^'5/3 ~ -^^apS-fS. (18) 

Permuting the [3 and 5 indices and injecting the factor 5af}5~fS obtained in this way in (|l8|) leads trivially 
to eq. (|l|). 

We note that the sign of the color octet term in the r.h.s. of ( p^ ) is the opposite of the sign of the 
color singlet. If channels are attractive for the singlet, they will be repulsive for the octet. This is what 
happens within the NJL model (for the GCM model, see below) where it is usual to neglect this octet: 
at the tree level it cannot contribute to give a colorless object while, anyway, the colored object would 
be unbound. 



or 
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qq and qq sectors 

It has been shown by CahiU and coUaborators ||l2|,^,^ that an alternative bosonization to meson modes 
qq is possible. It amounts to rewrite the interaction to build qq and qq diquark modeg|^ (which, when 
combined with a third quark, can give rise to baryons and anti-baryons). Although not conceptually 
necessary in the NJL model, the situation is different within the GCM model: since it is based on a model 
gluon propagator which makes the quark interaction non-local, bosonization implies the introduction of 
bilocal fields. With the Fierz identities of the previous section, this would amount to the presence 
of bilocal fields describing singlet and octet states. Physical mesonic states are obtained through the 
expansion of these bilocal fields into local ones. However, since the gluon exchange is repulsive between 
color octet states, the corresponding local expansion cannot be performed. We are then left to another, 
alternative, bosonization: now, the bilocal fields describe qq singlet Ic color states, as in the previous 
section, and qq antitriplets 3c color states, i.e. diquark states. This result is important because the gluon 
exchange between quarks in the 3c states is attractive. An expansion of bilocal diquark fields into local 
ones is then possible. These qq 3c states play an important role in the baryonic structure: a baryon is a 
color singlet formed of three quarks. Two of the quarks have then to be in a 3c state. This comes from 
the fact that, since the quarks are in the fundamental representation 3c of the color group 5'C/(3), two 
quarks can only form antitriplet 3c or sextet 6c states^ (3 (8) 3 = 3 ® 6). The third quark can only be 
added to the antitriplet and not to the sextet (3 ® 3 = 1 ® 8 as for the mesons, while[]3® 6 = 10 ©8). 
An easy way to make diquarks appearing in the formalism is to add and subtract 5a-y5sj3 in the r.h.s. 
of jl^), leading to 

(^)., (^).. ^ " k^'"'''' ^ ^WW, (19) 

where we used (t'')a5{^'')ii3 — Sa^ySsg — S^pS^s- {^pafs is a possible representation of 3.) 
The singlet color part 6ajSsi3 of (|l^) will be associated to flavor and Dirac Fierz identities connected 
to qq modes while the color octet part tpaSf-p-yp will be associated to flavor and Dirac Fierz identities 
connected to qq and qq modes (see below). 



2. Fierz identities in flavor space 



qq sector 



We have to go from singlet flavor d.o.f. to any flavor d.o.f. We just rewrite eq. (fi^) in flavor space, with 

1/3 '^57- 



the last term on the r.h.s. written as (1/2)A°^AJ5 



5^J5M=Y.^G%l{G^)u,, (20) 

e 

where (with F meaning flavor space) 

We have ordered in the way singlet + octet. 



These Fierz transformations leading to diquarks have also been used in the NJL model. See for example works 
by Reinhardt Hi and Weise [|||. 
^^For a recent review of group theory with applications to nuclear physics, see |28| . 

Moreover, CahiU et al. have shown that in GCM type of models - current-current interaction - Fierz trans- 
formations do not yield color sextet diquark states. 
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qq and qq sectors 

We reorder ( po| , ^ ) to correspond to diquarks in the (antisymmetric) and 6f (symmetric) representa- 
tions: 

5.Ai=Y.^G%^{G%,, (22) 

e 

where 

G =(G^,Gs)-(^,^,^,^l_^,^,^,^,^,^). (23) 



3. Fierz identities in Dirac space 

Since we do not yet specify the chosen gauge, we work exphcitly with D^^{x^y) instead of the more 
famihar choice corresponding to the Feynman-hke gauge g^,^D(x,y). 

Any 4x4 matrix can be expanded (we recaU that we use the notation of |^^) on the basis F" with 



scalar vector 


tensor 


axial pseudoscalar 


pS pV 


F^ 


pA pP 









With these matrices, the forms qT^q are hermitian. If we define F^ = (F") ^, we then have 

Tr(F"F/3) = 4(5;^, 1 < a,/3 < 16. (24) 
The matrices F^ are obtained through their definition: 



scalar vector tensor 


axial 


pseudoscalar 


Fs Fy Ft 


Fa 


Fp 


1 7p f^fii^ 


-7Sp 


-17^ 



Since any 4x4 matrix X has the expansion 

X = x„F" = iF°tr(XF„) = iF„tr(XF"), (25) 

one can then show that 

SabScd = ^F^,F„,<j6. (26) 



Four-quark interaction 

From &(4, 2; 3, 1) = g(4)F"g(2)g(3)FQ,/(7(l), we deduce (taking properly into account the Fermi field anti- 
commutation) 

6(4, 2; 3,1) = -qa,ma,m^na,a,iTa')a,a^qa,{2)qa^{l) (27) 
= -qaA'^)qaA^")a,aA'^a')a,a,qa',mqai{l)Sa',aJa',a,^ (28) 

With Sa'^a2 = ^Q2Qi and eq. we have 

6(4,2;3,1) = -ig„,(4)9-„3(3)(F"),,,,(F'5)„,,,(F„0a3ai(F/3),,,-<z,^(2)g,, (1) (29) 
= -ig-(4)F"F^(l)g(3)F„,F^g(2). (30) 
This is the basic formula in the mesonic sector. 
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Color quark current- current interaction 



The complete demonstration of Fierz identities in Dirac space is given in appendix ^ (where we explain 
the notation and physical contents related to 5). Here, we give only the main result, that we mix 
with sections I C 1 and I C 2 to obtain the complete Fierz-transformed interaction part of the Lagrangian: 



2 



2 



(l{x)—Yq{x)D f,^{x,y)q{y)—q{y) = --D^^(a;,y)- 



+ 



- ^ (q{x)G'^K'y{y)) [q{y)G^Kq{: 

a,e 

(q{x)GHrq{y)) (qiyjCHj-^qix)^ - 2(q{x)GH^^rqiy)) [q{y)GH^^rq{x) 
(^q{x)G^-^q{y)) [^q{y)G^ — q{x)) 
g^'^Y. {q{x)G\/-iK'^Cf{y)) (q^ {y)CG\/-^Kq{^) 



V2 

2(q{x)G%ir'-^Cf{y)) (q^{y)CG%'-^i^'^q{x)) 



(q{x)G\'-^^^^rCf{y)) (q^ {y)CG\'-^^^''^^ q{x)) 



{q{x)G% '-^f^Cf {y)) (q^ {y)CG% ^ 



^/2 V2 



V2 

p 

V2 



(31) 



with K"" = (l,^7^^7''',^7^7''',(T^'''7^/2). 

We can introduce meson (diquark) fields in the usual way. Because of the part of the propagator propor- 
tional to fc^fci/ (see below), we have tensor mesons as well as non-diagonal couplings between vector mesons 
(diquarks), between axial mesons (diquarks), and between tensor mesons (diquarks). For example, we 
get, for the meson part, 



1 



1 



-Z?^,(x,2/)-|g'"'r(a;,y)0<,(y,x) - 2r{x,vW{y.x) - 2^^(x, 2/)</)^(y, x) 



(32) 



or, separating each contribution, 

^ii^^,(x,y)i|g^''0(a;,y)0(y,a;) -|-g'"'05(x,y)05(y,a:) -I- (s^''^ (x, y)0„(2/, a;) - 20''(x, y)0''(y, x)) 
+ (g''>^ (x, y)^5,a (2/, x) - 2<^^(x, y)<l>l{y, x)) 
+ (5^"^C^(x, y)^.,o.p{y. x) - {x, yW^y {y, x)) | , 

where we have defined (f>''{x,y) = (0(x, y), ^5(2;, y), ^''(a;, y), 0^(a;, y), 0^''(x, y)) with 

(l){x,y) = qix)G''q{y), 
05 (a;, y) = qix)G''ij^q{y), 
r{x,y)=q{x)GHrq{y), 
4>^{x,y)=q{x)GH^^-l^^q(y), 



(33) 



(34) 
(35) 
(36) 
(37) 



r„''{x,y)=q{x)G^^q{y). 



(38) 
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4- Conclusions 



• From the previous section, it is clear that Fierz identities are a very useful tool to rewrite the 
interaction in terms of d.o.f. having the good low energy mesonic quantum numbers. It is also clear 
that, as soon as the most general form of the gluon effective propagator is retained, all quantum 
numbers are allowed (scalar, pseudoscalar, vector, axial- vector, tensoip^ . The above derivation 
has also shown an extra interaction, compared to the Feynman-like gauge, among each family 
(vector-vector, ...). Indeed, in the Feynman-like gauge, the gluon propagator takes the simple form 
D^i,{x,y) ~ g^yD{x,y) and, leaving aside diquark contributions, ( p2|) is replaced by 

- -^D{x, y)i (4(02 + 02) ^ 2(0^0^ + 0^05,^)) ■ (39) 

One effectively sees that there is no longer tensor term and that the interaction inside each species 
is simpler. In a covariant gauge where the gluon propagator writes (because of the translational 
invariance, the propagator is only a function of [x — y), i.e. of one variable k in Fourier space) 



D,Ak)-[9,.~^)^+^^, (40) 

it is clear that, together with the g^jy term g^^d{k'^)/k'^ which gives rise to the usual Feynman- 
like gauge result, there is an additional term depending both on the gauge parameter ^ and the 
particular shape of d{k'^), leading to interaction terms of the form k'^(j)"'(j)a, {k.cf))^ , ... 

Note also from eq. ( |39| ) that we recover the usual factor 1/2 when comparing the strength vector- 
pseudovector versus scalar-pseudoscalar. 

The NJL model can be seen as a limiting case of the GCM class where the interaction is i) written in 
the Feynman-like gauge and ii) is a contact interaction: 13^^ {x^y) cx g^i,S^^^ (x — y). The interacting 
part of the Lagrangian is then 

A„* cx i[(g-A„g)(9-A"g) + (qT^Cq^^q^CT'^q)] = L^t + C^t. (41) 

where we have defined (with K"^ — (1, 175, Z7^/V2, «7^75 V2)) 

A" = Ic^G'^^if", 

F" = ^e" (g) S «) X^ (42) 
v2 

It is useful to note that, in order to reproduce the pseudoscalar and vector spectra, different coupling 
constants have to be chosen as free parameters (i.e. the factor 1/2, founcf^ when comparing 
pseudoscalar and vector strength, is not suitable). 



When A'c 7^ 3, the factor 1/3 in the interacting Lagrangian (41) is replaced by [N^ — l)/(27Vc) 
for and XjN^ for li\nf When Nc ^ 00, Cl^^ is reduced by a factor l/Nc as compared to the 
mesonic part C^^^. We are then left with a pure mesonic interaction, in agreement with model 
independent large Nc analyses, so that the baryons have to be exclusively described as solitonic 
solutions of the Lagrangian (see for example the works by Reinhardt and collaborators (l^j3^ , and 
the works by the Bochum group |33-35j). In the last reference, the interested reader will find results 



concerning mesons, baryons, and links with chiral perturbation theory and the Skyrme model. 



^^This has also been seen in the works by Shrauner |2^] and Roberts and Cahill |30| 

^^This factor is not constrained by the symmetries of the model, so that choosing different coupling constants 
still respect them. 



13 



• Apart from solitonic solutions, baryons can also be obtained from (pl|), treating them as a quark- 
diquark system. When a third quark is added to the diquark, one has to consider states formed by 
the product (3 (g) 3)c «) (3 «) 3)f, which leads to (1^ » If) ® (Ic » 8f) © (8c «) If) © (8c «) 8f). 
Calculations performed by Cahill have shown that for the {TLc^^f) and (8c(8i8f) states, quark 
rearrangement is repulsive, while it is attractive for the remaining terms. The state (8c (E" If) has 
then to be suppressed by hand. This defect, present both in NJL and GCM studies, can be traced 
back to the lack of local color gauge invariance. 



D. The 3-flavor scaled NJL models and anomalies 



We showed in the previous section how to rewrite the color current-current interaction as a function of 
bilinear forms having the physical quantum numbers seen in the spectrum. This procedure amounted to 
use Fierz identities in color, flavor and Dirac (spin) spaces. We also mentioned that the NJL model can 
be viewed as a limiting case of the GCM model, where the interaction is local (4-quark point interaction) 
and a Feynman-like gauge is used: eq. ( pi] ) together with Fierz identities limited to mesonic scalar and 
pseudoscalar modes leads to 

Cnjl = m - m)g + Gs Y}^q^^qf + (9*75^g)']. (43) 

The NJL model has already been extensively studied by several groups of which we can only quote a few. 
Our purpose being here to describe its scaled version, the reader is referred to the following literature]^ 
(and references therein) to more general studies about the model. The most recent work has been done 
by Ripka in a book which contains an in-depth analysis of regularization procedures and symmetry 
conserving approximations; Klevansky |3^ and Hatsuda and Kunihiro ||3^ give a general introduction 
to NJL, both in the vacuum and at finite temperature and density; Alkofer, Reinhardt and Weigel |^ 
mainly applied the model to discuss baryons as chiral solitons, as also done by Goeke and collaborators 
[Q who discuss, using heat kernel and gradient expansion techniques, the link between NJL and fully 
bosonized approaches of the Skyrme type on the one hand and the chiral sigma model on the other hand; 
Bijnens discusses chiral perturbation theory within NJL; Alkofer, Ebert, Reinhardt and Volkov | |lO| , [T^ 
hadronize the NJL model (both mesons and baryons p^ th rough the use of the alternative Fierz identities 
leading to diquarks; Vogl and Weise and Weise [^Tf review several of the above mentioned topics: 
bosonization and hadronization, finite density and temperature effects. 

It is clear that the Lagrangian (^3|) describes the interaction between flavor carrying colorless objects. 
In the chiral limit where the current quark mass matrix m is identically zero, it is invariant under 



SU{2,)v ® SU{i)A ® U{l)v ® U{\)a (see notation in section I A). However, the U{1)a invariance is 
broken by quantum effects (strong axial anomaly) which cannot be taken into account since the true 
QCD Lagrangian has been replaced by an effective form. The Lagrangian ( p^ ) is used at the tree level. 
If we want this order to implement the full quantum aspects of QCD, it is necessary to supplement it 
with a term which, while still invariant under the true symmetries of QCD, has to break the axial U{\)a 
invariance. This anomaly can be related to the formation of instantons p8|-p0t and yields anomalous 
contributions to the rj and r]' masses. The U{1)a symmetry is explicitly broken by instanton induced 
interactions which have the form of a 't Hooft determinant. The phenomenological term which has to 
be added to (E^) can thus take the form of a 't Hooft determinant. Other forms can be chosen (see 



[ }48| , [f9| , pl| -|54|]). For example, working on the bosonized (see section ID 2) version of ([43|), one can show 
that the Schechter iQ propositioiiPf 

Cu(i)^ (X (trF ln(o-<,A'^ + iTTaX") - trF ln(craA'^ - iTTaA"^))^ , a 0, 8, (44) 



^^The number of published papers in this field of research is enormous and still growing. We apologize if 
important contributions are omitted here. Since our main purpose is to make an introduction to the scaled NJL 
model, we mainly mention review papers and school proceedings. 

^^In this way, baryons are no more pure solitonic solutions but are also described by a three valence quark 
Faddeev equation. See also p^ , ^o|j4l| for this approach and p2[-|46[ for numerical results. 

^^TTa and da will be introduced in section IDS. They are related to pseudoscalar and scalar d.o.f., respectively. 
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is also compatible with the symmetries of QCD, while breaking the U{1)a one. This Lagrangian is 
interesting, being simplerpi than the 't Hooft determinant, while containing it as a limiting case. In 
particular, ( |4^ ) does not modify the gap equations. Up to second order expansion into the mesonic 
fields, it yields extra mass terms for the 77 and 77' particles. We shall use here an even simpler approach 
which consists in just adding a mass term for the 770 particle (i.e. the ttq in the above notation). In 
the chiral limit, all the pseudoscalar particles would then be massless (in agreement with the Goldstone 
theorem applied t o the dynamical chiral symmetry breaking mechanism) except for the rj' (770 in this 



limit), see section ID 4. This procedure is related (although not equivalent) to the approach of Frere 
and collaborators [ 35|,^6[ : in these works, there is a coupling between the 770 and r^g through the quark 
mass difference and the anomaly whose net effect is to modify the mixing angle from its ideal value. 
It can be shown that, in the mean field approximation, our approach only explicitly breaks (apart for 
quark masses) the U{1)a symmetry, then being a possible parameterization of the axial anomaly. The 
expression for the divergence of the axial current is given in (^ where the last term is the anomaly. It 
has the quantum number J^'" = 0^^ of pseudoscalar isoscalars. It is then tempting to identify it as an 
interpolating field for the pseudoscalar glueball and to incorporate it into the formalism so that the 77 and 
?/ have to be searched from a three state (77, rj' , glueball) diagonalization procedure^. We shall however 
restrict ourselves to the r^g-mass procedure. 



So far we only discussed the axial anomaly. As mentioned in section I A, another symmetry at the 
classical level of the QCD action is broken by quantum effects: scale symmetry. In fact, this symmetry is 
broken in three ways: explicitly, through the quark masses which introduce a dimensional parameter in 
the theory (see the first term on the r.h.s. of eq. (|^)); dynamically as for chiral symmetry down to vector 
symmetry; because of the anomaly, see the last term on the r.h.s. of eq. (|^). As for the axial anomaly, 
the effect of the trace anomaly has to be added by hand. In that way, QCD quantum effects are present 
to the lowest order of the effective theory we are interested in. 

Several steps are necessary in order to construct the modified NJL model in the perspective of the 
symmetries and anomalies as described above. In view of the way of treating the axial anomaly (mass 
term for the 770), it is better to work on the bosonized version of the model. The main points to bosonize 
the theory are sketched below. 



1. Bosonization 

For a review of hadronization techniques, the reader is referred to ]60|j6l[ |. 
We start from 

Zqcd~ JvqVqcxp^i j d^xC^jL^ = Znjl, (45) 

with Cnjl given in (^3|). Unity is introduced in (^), through auxiliary fields and (p°' 

1 ^ J V(f>''S{(f>'' - qV^q) ^ J P^'^P^" exp (^i J d'^ xip" {(j)'' - qV'q)^ , (46) 

where cjf — ((T",7r'^) and = (A'^/2, i75A''/2), a = 0, ...,8. (The first component is relative to scalars, 
the second component to pseudoscalars; to simplify the notation, we have removed the fiavor index F .) 
Eq. (^) is ~ apart form unimportant constant factors - the functional equivalent of a possible definition 
of the Dirac distribution. From f{x)S{x — xq) — f{xo)5{x — xq)^ we get 



Znjl = j VqVqV^'D^" exp (^i J d^x[q{i^ - m")g + Gs(</)")' + - f"q^"q] 
With the identity 



(47) 



Gsirf + = ^5(0"^ + ^^")^ - ^(^")' = Gs{rf - (48) 



^^This form does not produce any fiavor mixing, while the 't Hooft determinant implies maximum flavor mixing. 
^"This is similar to the approach of Shore and Veneziano Ist). See also the works of Nekrasov pqJSQl 
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the (/)'° field integration is Gaussian and can be performed exactly. (The Jacobian of the transformation 
0° — > 0'" is unity.) Since its value contributes only a constant term it can be omitted. The effective 
functional takes the form 

Znjl = j VqVqD^'^ exp J d\[q{i^ - m° - v^T'^)? - ^(v'")']) ■ (49) 

Then, the quark Gaussian integration has to be performed (for Grassman variables q and q), so that the 
model describes only mesonic d.o.f. 

We already mentioned the fact that the GCM model (and studies performed with SD equations) is best 
expressed in the Euclidean space. This is also the case here, although not for the same reason: since we 
are interested in finite temperature effects, we found that the Matsubara (or imaginary time) formalism 
is well suited for our purpose so that the Euclidean space will be used from now on. The functionals (^) 
and ( ^ ) become partition functions describing the thermodynamics of the systems. Formally, we go first 
to Euclidean metric, then making the substitution]^ t = —it with real t. We have 



2 2 I ^2 

X = T + X , 

7^ - (70 =«/3,7), (50) 

where (3 is the Dirac matrix^ diag(l, 1,-1,-1) ~ see appendix ^ for our conventions -. 
We define ||,|| 

TrO = tr J d'^x{x\0\x) , J d'^x = J^\t J d^x, (51) 

where tr is the trace w.r.t. internal d.o.f. (Dirac, color, flavor), f3 = 1/T is the inverse temperature and 
57 is the volume of the system. We get 

J V^^exp (^-[-Trln(-^^ + m" + (^"r'') + J d^^^i^'f]^ , (52) 



-■NJL 



where we have defined the NJL effective action in the Euclidean space 

1 



/e// = -Trln(-i^ + m" + ^'^r'^) + / d^x-^i^"")''. (53) 

4:Gs 



Following the standard procedure |62|-|64[| , density effects are considered through the introduction of 
a chemical potential for each flavor (Lagrange multiplier for the quark number conservation). In our 
formalism, the equivalent relation to H — fiN, where H is the Hamiltonian of the system and /i the 
Lagrange multiplier associated to a conserved quantity N, consists in subtracting the quantity^ q/J-Pq 
from the Lagrangian. In Euclidean metric, this is equivalent to the introduction of a vector Wi, = 
{—ifj,, 0,0,0), where /i is a notation for diag(/iu, /x^i, ^<,)- Each fii (i — u,d,s) is a chemical potential 
associated to the corresponding flavor. When all chemical potentials are equalFj, or in the case of a single 
flavor, the usual result nqPq is recovered. 

With these prescriptions, the effective action describing a system of quarks at finite temperature and 
density is 



-Trln(-z^ + m" + v^'F'^ - W) + f d^^^iv'f, (54) 

J 4Gs 



which still has to be modified in order to implement the QCD axial and trace anomalies. 



^^In the NJL model, this substitution does not cause any of the problems mentioned at the end of section 
^^There should be no confusion with the inverse temperature. 
^^/3 is the Dirac matrix and /i a diagonal matrix in flavor space. 

^^Note that different chemical potentials break the invariance under axial transformations. 



[B. 
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We already mentioned in the beginning of this section that a mass term for the rjo (or ttq in our notation) 
is well suited and particularly easy to implement. We have thus 

leff = -TrH-t^ + m^ + ^ara~W)+ J d^^Y^^l + J d^x^^Trl (55) 

where, to make the notation similar to the original papers ]6^ , |6^ , |65t , we have introduced the coupling 
constant = 2/(Gs) and have redefined T"- = (A", i75A°), a = 0, 8. 

The scale anomaly requires several steps to be included in our formalism. First of all, we have to make 
the NJL action scale invariant. This comes from the fact that without anomaly, the QCD action is scale 
invariant except for the quark mass term. The mass term in NJL has the same scale dimension as in QCD 
so that this term has not to be modified. However, the coupling constant in (|5^) has a dimension and 
has then to be modified. The second step is connected to the anomaly itself. It is generally modelized 
by a scalar dilaton field. This has been used to modify different models of QCD at low energy. For 
example the Syracuse group [ p6| , |67[ studies in this way scalar mesons and chiral solitons. More recently, 
this dilaton field has been used to render non-linear a models scale invariant for examining the effect of 
density on meson masses |68 and the deconfining phase transition |^,^. The authors of |71-74| used a 



modified linear sigma model to look at the mesons and nucleons at finite temperature and density. Cohen 
et al. studied the behavior of quark and gluon condensates w.r.t. density in a model independent 
way (at low density) and within NJL. Kusaka and Weise Jt^ study some hadronic properties, mainly 
as a function of density, and indicate conditions which have to be satisfied in order to get a scaling of 
these properties. They show in that work that the universal scaling as proposed in p8[ | can only occur 
for small glueball mass (compared to the mass of the chiral scalar mesons). They show, and it has been 
confirmed in ||7^ , that the low-density behavior of hadronic properties is insensitive to the specific form of 
the interaction. We shall comment further on this point later. The condition that the glueball potential 



must satisfy to reproduce the anomaly is described in 66|. In the works p9| , [70[ , the dilaton field is 
an order parameter associated with gluon confinement. Wc shall adopt this point of view, mainly to 
facilitate the language, while keeping in mind that it is questionable [ p^ . 

At the NJL model level, the dilaton field is put by hand in the interaction term in order to restore the 
scale invariance. Let x be the dilaton field, of scale dimension 1. The interaction term 

Jd'x^x'vl + jd'x^^x'-l (56) 

where the coupling constant has been replaced by a^x^, of scale dimension 2, is then scale invariant^. 
With this procedure, the NJL model is scale invariant^ except for a mass term of scale dimension 3, as 
in QCD. 

One drawback of the NJL model is its lack of renormalizability, which leads to a subtlety in the scaled 
model: the regularization parameter A, chosen is this work as a 4-dimensional cut-off, cannot be elim- 
inated from the theory. In this way, it is a fundamental parameter of the model. This cut-off A enters 
the game when one wishes to evaluate the fermionic determinant ( p5|) and breaks its scale invariance. 
The work of Ripka and Jaminon fz^ is of prime importance in the sense that the scale invariance of the 
fermionic determinant is kept by introducing the dilaton field in the regularization factor: A — s- A^. The 
effective action is then 

hfS = -TrA^ H-i^ + m° + ^aTa - W) + J d^x'^x^l + / d^x'^ix^i^l (57) 

where we have explicitly indicated by TrA^,- the fact that the regularization parameter depends on the 
dilaton. 

Up to this point, all we have done is to make NJL scale invariant. We still have to consider the anomaly. 
It can be introduced by adding the Lagrangian 



^^The new coupling constant has now no dimension. 

^^In the field strength method described in [ p"o|Jl4| , the first term of eq. ( |5^ ) can be derived in the NJL limit, 
i.e. the scale dimension is correct and in relation to the gluon condensate. It is then legitimate to argue that it 
is the elimination of the gluons from the theory which gives rise to the efi'ective interaction (pd) . 
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with 



yx-^{x'H^)~ix'-xh)]- (59) 



The trace of the energy-impulsion tensor associated to the Lagrangian (El) is 



2 



2g ^ t^v^ a dx 4: 

where the first equality comes from (|^). (In the last equality, we omitted the subtraction of (&^/4)xg; 
coming from the last term of ( |59| ) , included there only to make the potential equaling zero at its minimum 
X = XG-) 

2. A and B scaled models 

Finally, the bosonized, scale invariant (except for the quark mass term) NJL action is 

r ^2 

hff = -Ttax ln(-«^ + m" + ipaTa -W)+ d^x—x\l 
+ I d'-^^^x'nl + '-{d.xf + I {x'H^} - ix' - Xh)) }, 



(61) 



taking into account axial and scale anomalies. This actionp^ called A-scaling model, is the simplest one 
verifying the above mentioned conditions of symmetries and anomalies. In the works ||6^ , |63| , |65| , [77|j79[| 
several models have been studied, differing each other by the way the quadratic part ip^ is made scale 
invariant. For example, the B-scaled model respects all the symmetries of QCD, and in this sense is as 
valid as the A-scaling model (pi]). The 3-flavor extension of the B-scaling model writes 



2 



Je// = -TrAx ln(-*^ + m" + ^.F, ~ W) + J d^xi^-i^i + ^n'^ - ^x'Y 

+ \{d,xf + ^ (x'ln(^) - (X' - XS)) }• (62) 

The justification of this model was to try to reproduce the scaling behavior of the hadronic properties as 
proposed in [|8| . However as was noted in (in the NJL model) and in ^ (model independent study) , 
universal scaling is not expected to occur. This can be traced back to the high value [Mcl ^ 1500 MeV) 
of the glueball mass. In the following we shall restrict ourselves to the A-scaling models. For results and 
comments upon the B-scaling version, the reader is referred to [ |63[ . 

3. Gap equations (or Schwinger- Dyson equations) 

From now on we shall use the model (|6l|). We work in the spirit of the effective model: the tree or 
classical approximation is used throughout the study, which corresponds to the lowest order in a 1/Nc 
expansion. (For a study at next order in Nc, see [|80|-p6t.) We work in the saddle-point approximation 
which consists in evaluating the effective action at the point in (Ta,7ra space where its first derivative 
vanishes identically. More generally, the saddle-point approximation consists in Taylor expanding the 
effective action, putting the first derivative equal to zero and evaluating the higher order terms in the 
standard way. For example, the second order is Gaussian and can be evaluated exactly. Note however 
that, to get the effective potential, one should reactualize the expanding point at each order |^^. We 
shall work with the lowest order. This means that we expand (^ij) around its stationary points 



Without the ^-term, the model described by this action has also been studied by Kusaka, Volkov and Weise 
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51 

5<Ta 
SI 

SI 

Sx 



= 0, 



0, 



(63) 
(64) 
(65) 



These equations lead to a parity invariant vacuum tt* = 0. The first one is called the gap equation (by 
analogy with superconductivity) or the Schwinger-Dyson equation for the quark self-energy. The last 
equation, which is not present in usual NJL models, is a consequence of the introduction of the dilaton 
to modelize the trace anomaly. It can also be shown that, out of the nine equations summarized by (|63|), 
only a = 0, 3, 8 leads ^ 0. Moreover, if the isospin limit to„ — nid is considered, then we also have 
0*3 = 0. In the following, we shall work in this isospin limit so that we have to manage three coupled 
equations for dg , cTg , Xs ■ 

Working at tree level means that the field stationary values will be used in order to 



determine the thermodynamical properties of the system (see section [I A ) 



determine some static meson properties (see section ID4) 



study some mesonic processes such as the neutral pion decay into two photons (see secti on [II C ) or 
scalar isoscalar meson decay into two photons, including the scalar glueball (see section [II C| ) . 



The last two mentioned points require introducing the electromagnetic field in the formalism and expand- 
ing the action up to third order or, equivalently, working with triangle Feynman diagrams (see sections 
il^ and pT^ ). 

The two (x9) gap equations (63 6^) lead to (the calculations are well documented in | |63| , |88t ) 



with 



9Mi,l3,tj.i 



0,a = 0,...,8, 

0, a = 1, 7 (ct| = in the isospin limit), 

8NcMsgM,,i3,tj.,, 



-y- 



Mi 



u,s, 



(66) 
(67) 

(68) 
(69) 

(70) 



and k* = (fcp — 'ifJ-i,k),kQ = (2n -|- l)7r/3, the odd Matsubara frequencies being related to the anti- 
commutating nature of the fermions. The constituent quark masses M„ and Mg are related to cro''''8 

by 



1 



Mu = -=(cr^ -Hma) 



"^1 + \/ a'^o^ 



2 2 
Ms = -■^('^s + "^3) + mi + y -erg. 



(71) 
(72) 



As already mentioned, the NJL model is not renormalizable, so that the regularization parameter cannot 
be removed. Several procedures are available p6t . Here we shall focus on a 4-dimensional regularization 
scheme. One of its advantages is that it does not break the Lorentz invarianceQ of the theory at zero 



^^Regularizations with a 3-momentum cut-off A break it. A meson, described as a qq excitation, cannot liave 
a center of mass momentu m | fc| exceeding 2A. Tliis is unpliysical and leads to distortions of vr-mesons thermal 
distributions. See however |B9| for a way out of this difBculty. 
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temperature and density. Another one is that it gives the correct (since the temperature and density 
dependent part is not regularized) high temperature Umit^ corresponding to a plasma of free massive 
quarks. The regularization procedure chosen here is the following: 

Let A{f3,fi) be a temperature T — 1/(3 and a chemical potential /i dependent quantity. We can always 
write 

A{(3, fi) = A(/3, fi) ~ A{oo, 0) + A{oo, 0). (73) 

All the terms are diverging quantities. However, the subtraction between the first two gives a finite result 
which does not have to be regularized^] (in fact, it cannot be regularized with a 4-dimensional cut-off 
because of the heat bath as a preferred referential). The last term is a zero temperature and density 
quantity which is Lorentz invariant and can be regularized with a 4-dimensional cut-off (in Euclidean 
space). We write 

A^i(3,fi)=A{P,fi)-A{^,0) (74) 
as the Fermi part of the quantity A, while 

A^(oo,0) = yl(oo,0) (75) 

is the Dirac part. 

This regularization procedure will be applied to every diverging quantity we shall encounter in the follow- 
ing. Note however that the pressure (and hence the action) has to be specified with respect to a reference 
point. We choose it to be the non-perturbative vacuum. 
With these prescriptions, ( [70| ) can be evaluated and gives 

gMi,l3,tj.i = gMi,f3,fii — 5Mi, 00,0 + 5Mi,oo,0 (76) 

with 

d^/fc 1 



2 



(27r)4p + M, 



167r2 fc2 + M,^ 167r2 



(Ax.f-M.^ln(l + ^ 



2 



(77) 



and 



1 ^ r°° 



{2Try [{2n+l)TT/(3^itMY +E? 
dko f°° d^k 1 



(2^) Jo {2nf fco' + Ef 
where we have defined the single-particle energies E., = [k^ + Mi^f/^ and 



(78) 



'"'^^ l + exp[/3(i?,±M0]' ^^^^ 
The gap equations are intimately related to the quark condensates since they are given by 

I dl^ff 1 
{uu) = j^-Q^ {Mu. Ms,Xs) = - (Mu - m„) = -AN,AUgM^,p,^^, (80) 

1 Orrr 1 

- (M«, M,,xs) = -Tra'x' {Ms - ms) = -m,MsgM^,p,^^. (81) 



pn dills V--"'— '^^^ 2 



^^For a 3-momentum A cut-off, the temperature part decreases as \/T as T is increased, see, e.g., | pl| . 
^°The subtraction in ( [73|) also leads to diverging results when higher order loops are taken into account. 
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Note that these condensates are not independent, contrarily to the pure NJL model: they are related 
explicitly (see Xs in ( ^oUsi] )) and implicitly by thegluon condensate (implicitly because the value of M„ 
and Ms is extracted from the gap equations (^,|9|) which are coupled to the equation for the x field - 
see below). 

These quark condensates show that a quark condensate can be considered as a good candidate for the order 
parameter of chiral symmetry breaking. When this symmetry is explicitly broken, the quarks acquire 
constituent quark masses Alt > rui, solutions of the gap equations. Because of the explicit breaking of 
the chiral symmetry by the current quark mass mu,ms, it is better to set them to zero to study its 
dynamical breaking. We then see that the up quark gap equation (|6^) (and similarly for the strange 
quark) admits two solutions: one with M„ = ruu = 0, the other one with M„ ^ 0. When the solution 
My. ^ is energetically favored, the chiral symmetry is dynamically broken and the quark condensate 
does not vanish. The solution Mu = can be reached at high temperature or densities: the quark 
condensate vanishes and the chiral symmetry is restored. If the chiral transition is not discontinuous 
(what we shall call second order in the following), the critical exponents are well known since we work at 
the mean field level. It can be shown that there is a correspondence between the quark condensate and 
the magnetization, and between the current quark mass and the magnetic field. This means that when 
rriu = 0, the shape of the quark condensate near is like (Tc — T)^/^ (where is the critical temperature 
at which the quark condensate vanishes) while there is a long tail extending from Tc if m„ ^ 0. 
However, the phase structure of our extended NJL model is richer than in the usual NJL model. Tech- 
nically this is due to the fact that the gap equations are coupled to each other and with the equation for 
the X field (pa), which reads: 



Mr 1 



(Axo)VA/,f 



62/2 
Xo 



exp ■ 



2x. 



MO 



2JVc A"^ 
8x2 



2 I M? 



(Axo)' + Mf 



"T-ti\^ M'i f ^ TO, ^ ^ 

M^ 1 -] H ^ 1 : 

M,J 2 I M, 



(Axs) +M„^ 



X. 



62/2 



{KXsf + Ms 



(82) 



where Xs has been elimina ted in favor of y p, t he vacuum expectation valuePjof 
We show in section [I B 1 (see also section |^) that the coupling between (]68|,p9|J 
first order transitions, depending upon the free parameters of the model. 
As a conclusion to this section, it is useful to note that 



at T = ^ = 0. 
\} enables also to get 



• in field theory, gap equations are called Schwinger-Dyson equations for the self-energy. Defining 

I] = m°-Ai/3 + (^^ra, (83) 
one can effectively show that (E^) is equivalent to 



1 



{aXs 



Tr 



Fa 



(84) 



gap equations, written in the form (68) and ( |69[ ) or (p4D, show clearly that the coupling constant 
can be redefined by absorbing the color factor Nc- This suggests that, in the limit Nc — > oo, the 
effective coupling constant is 



This is related to the above mentioned fact that the pressure is only defined w.r.t. 
point is chosen to be the vacuum at zero temperature and chemical potential. 



a reference point. This 
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(85) 



In term of the strong coupling constant g, we find 



lim g Nc = est. = g 



(86) 



which is in agreement with the high Nc analysis pO| ; 

• high Nc arguments can be used to justify the saddle-point approximation used to derive the gap 
equations. If we admit (p5|), we have, leaving aside the dilaton contribution, 



leff =-Nc^^TT'^H-i^+m° + ipara-W)- J d^^^^^l^^ 



(87) 



where Tr' indicates that the color trace has to be omitted, allowing a factor Nc to factorize. Besides 
justifying the saddle-point approximation, this indicates that Nc is a relevant factor to organize 
the higher order correction series [ p0| , pl| , ^3| -|8^ . However, to recover the low energy theorems, 
a particular treatment has to be applied to these higher orders |8^ , related to the symmetry 
conserving approximation program l3q,pq] . 



4- Equations for meson masses (or Bethe-Salpeter equations) 

As far as we are interested in meson masses, wave function normalization and coupling constants, we 
have to expand the action up to second order in the fields. Expanding up to second order means that 
the operator sandwiched between meson fields can be interpreted as a meson propagator. This can be 
cast in the form of a Bethe-Salpeter equation which, in the NJL model case, is tractable because of the 
inherent simplicity of the model: both the kernel and the vertices are known, so that we only have to 
search for the eigenvalues which are the masses entering the meson propagators. (For a discussion of the 
Bethe-Salpeter formalism in the GCM and NJL models, see pl|.) 
Expanding the action (|6^) up to second order gives^ 



with 



1 

2ph 



t2 _ t2 
^eff - 



ps 



[^0q,^8q,Xq,aQq,KQ^^ 



I pOO p08 pOx 

p80 p88 p8x 

pXO px8 pXX 



\ 







paoao 










p/f-x; 



O'S-q 
X-q 
So-q 
V K*Oa-q 



(89) 



and 



Ips - X! ^Sg: T^qiKaq) 




(90) 



which shows that, in the pseudoscalar sector, ttq and TTg, corresponding to ?7o a-nd "q^, are coupled in 
agreement with the fact they have the same quantum numbers (pseudoscalar isoscalar J 



PC 



= 



while, owing to the glueball, the scalar sector has three coupled particles^ (scalar isoscalar J^^ = 0^^). 



^^The quantities P and F are given in appendix ^. 
Three coupled particles can also be described in the pseudoscalar channel by allowing the presence of the 
pseudoscalar glueball J^^ — 0""*". This is for example the approach of the composite operator formalism of 
Nekrasov |S8|j59| . See also the works of Schechter |5^ , p^ . 
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The meson masses can be obtained from (B9 



Only pseudoscalar pions and kaons, and scalar K* 
and ao can be directly extracted. Masses for 77,77' and the three scalar isoscalars have to be obtained 
from a suitable diagonalization. In the isospin limit and for symmetric matter, these masses are given by 

(A^ = pLu- Us) 



ml{q) = 



9 9 



2 



2 2 



2ZK±{f3,fi) \Mu 



Ms 



(91) 
(92) 



for the pion-kaon sector [Kq has the same mass - except for isospin breaking corrections that we do not 
consider here - as while Kq has the same mass as K^), while the 77 — 7/ sector is given by 



■2. I \ 2 2 

772^,,^ (g) = a X, 



X 



1 



3 cos 29 q 
3 cos 29 a 



-1 



3e 



(93) 



where 9q is the solution of 



tg261q 



2P' 



08 



poo 
2^2. 



3C(aXs 



g2 (Z« - Z^) + {arf (f^ - ^ 



In these relations, the functions Z are the field normalizations. Defining 

^^^^""-^'^'^ - h ? {kf + M^) [(^; - qy + M/] 



(94) 



(95) 



these normalizations are 



Z^±(/3,/i) = 47Vc 



{/3,fi) = m,F^p{M^,M^), 
F^p{Mu,Ms) 

Zl EE Z,(/3,m) = 4A^,i;^^^(M„,M„), 
= 4A^,i;^^^(M„M,). 



(96) 

(97) 

(98) 
(99) 



For the scalar sector, the masses are (there is the same correspondence between the charged and the 



neutral mass as for the kaon) 



2 2 



(77i^..± (g) ± A/i)^ = (Af„ + M,)" + 



Z,(/3,//)M„ 

2 2 

2 , a Xs 



2^k±(/3,m) 



777^ 

Ms 



(100) 

(101) 



The solution for the isoscalar sector is not illuminating (the masses are solutions of a cubic e quatio n) so 
that we do not write it here. A numerical evaluation of these masses will be given in section III A| . 
These formulae can be used to study the mass variation as a function of temperature and density (chemical 
potential) and are the starting point for studying meson propagation in hot and dense matter. 
Formulae (pl|-p3|) are interesting because they readily show the Goldstone nature of the pseudoscalar 
mesons, and the role played by ^ as a modelization of the strong axial 1}(X)a anomaly. Indeed, working 
in the limit fJ^u ^ fJ'd = fJ's, we note that when ^ 0, we get 



777,^ 



2 _ Q^Xs rus 
" Ms- 



(102) 
(103) 
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Moreover, in the vector symmetry limit rriu = = m^, there is a complete mass degeneracy 



(104) 



In the chiral limit all these mesons have a vanishing mass, as it should (Goldstone theorem). 
If the strong axial C/(1)a anomaly is taken into account, we have (in the chiral limit) 



9 2 

a Xs 



(105) 
(106) 



The particle r] is still a Goldstone boson, but the t] has lost this property (in a world with three colors Nc 
= 3). Large Nc studies [^,^ indicate that the axial symmetry is restored when Nc oo. This implies 
that the ^ parameter behaves like l/Nc- This can can be seen by writing 



(107) 



where we have used eqs. (|85|) and (|9^). It is clear that a vanishing of to,,/ at large Nc requires ^ = ^/Nc- 
In figure |l|, we show the mass variation of 77 and 77' as a function of the 1^ parameter (for Af ° = 400 MeV, 
Xo = 350 MeV). It is clear that, without the axial anomaly, the 77 is degenerate with the pion while the 
77 mass is identical to that of tt with the replacements ruu ttls , Mu Ms . 




FIG. 1. 77,77' mass variation as a function of the ^ axial anomaly parameter. 



E. Parameters 



The main aim of the model is, starting from zero temperature T and density p, to study the evolution 
of the system with respect to these thermodynamical external parameters. The model has intrinsic 
parameters that we can fix at T = p = so that the evolution is a prediction. In the isospin limit, 
the 3-flavor version of the model has seven parameters: the current quark masses mm iris, the coupling 
constants a^, 6^, the vacuum gluon condensate yq, the cut-off A and the 770 mass parameter ^ modelizing 
the strong axial anomaly. The gap equation ( |68| ) in the vacuum can be used to eliminate the coupling 
constant in term of the constituent up quark mass M°. We choose and xo as free parameters; the 
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five remaining parameters are adjusted to reproduce, in the vacuum, the pion mass^ rnvr = 139 MeV, the 
kaon mass Mk = 495 MeV (same remark as in footnote 34), the ry' mass m^' = 958 MeV, the glueball 
mass]^ moB = 1300 MeV and the on-shell pion decay constant. These parameters are determined in the 
following way: once M° and Xo are given, the cut-off is obtained from the approximate^ equation for 
the pion decay constant 

U^Mlzl'^{^,Q)^2=_mi. (108) 

The pion mass (^ fixes the product a^m„ and a? is obtained from the gap equation ( |68| ) in the vacuum. 
Since we are not interested in the scalar sector, the parameter is then obtained from the glueball mass 



in an equivalent sigma model: TOq^ ~ Xo(^^ ~ SNcA^ /2tt''') (see section [11 for its exact determination 
from the mass value of one of the /o hybrids). Mg and are obtained from equations ( |69| , p^ and, 
finally, ^ is fitted to reproduce the r]' mass. 

The values of the parameters, together with the quark condensates, are given in table |, adapted from 
|63|, for four sets of input (M°,xo)- The first two sets have been studied in [|62|j7^. The choice xo = 350 



MeV is estimated from QCD sum rules applied to charmonium data [ |95| , |96| . The coupling of the quark 
and gluon condensates is then weak and the dilaton field is frozen to its vacuum value^j (we can say that. 



with this value of the gluon condensate, QCD does not scale, see |75 ). Lower values such as xo = 125 
MeV (which gives a bag constant B'^^^ « 200 MeV), have to be used to enhance the coupling. Low values 
of the gluon condensate are required to obtain a chiral phase transition at relatively low temperature 
(T « 140 MeV, as suggested by lattice calculations JM). This justifies the use of the last two sets of 
parameters of table |[ which were first studied in | ]6^ , |65| . 





TABLE L Parameters 


of the A-scaling model (moL 


= 1.3 GeV; 




xo (MeV) 


350 


125 


80 


90 


M° (MeV) 


450 


600 


300 


800 


niGL (GeV) 


1.3 


1.3 


1.3 


1.3 


nin (MeV) 


9.1 


8.7 


7.7 


7.7 


TUs (MeV) 


205 


196 


186 


177 


M° (MeV) 


645 


761 


539 


929 


A 


2.06 


5.85 


10.13 


8.59 


~ < UU >Q 


207.9 


211.4 


219.8 


220.2 


— < SS >q' 


207.8 


208.2 


234.1 


216.4 


2 

a 


0.333 


2.044 


11.359 


3.327 




22.05 


642.6 


5062 


2691 




1.09 


0.639 


1.3 


0.385 



Note that a special treatment has to be reserved to ^ for the two sets of parameters Mu ~ 450 MeV and 
M„ = 300 MeV because of the threshold problem due to the lack of confinement of the NJL model (the 
meson can unphysically decay into two free quarks). We explain this problem, and the way to still extract 
results, in section II C. In figure we show the variation of the critical^ temperature as a function of the 



It is well known (see for example ^]) that the mass difference between neutral and charged pions is of 
electromagnetic nature. Isospin effects account at most for a five percent difference. This justifies neglecting 
these corrections. In order to improve our calculations, it is necessary to look at these electromagnetic corrections 
rather than at the isospin breaking ones. Note also that we use, as always in the literature, the value m,r = 139 
MeV for the pion mass in the vacuum, even if we are interested only in neutra l pio ns. This will have to be 
corrected for effects which are strongly dependent on m,r . See for example section II C . 

There is some debate in the literature concerning the glueball mass. It is now believed that the scalar glueball 
mass IS 1500 MeV 11 or 1710 || (see also Q). As we are not interested in the scalar sector, the exact value is 
not important for the analysis. To make contact with our previous works we shall keep the value 1300 MeV until 
section 

36 



[II 



We have neglected a term q^dZdq^ evaluated at 
With such a value c 
By critical temperat 
of the chiral symmetry. 



^'^With such a value of the gluon condensate, the model is almost a pure NJL. 

^^By critical temperature and/or density we mean the temperature and/or density associated with the restoration 
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two free parameters in the two degenerate flavor version of the model. This picture is instructive since it 
shows that 
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FIG. 2. Critical temperature as a function of the parameters M„, xo- 

the scaled model is able to reproduce a low critical temperat ure 140 MeV) contrarily to a pure 
NJL model (T « 190 MeV) - see figures0| and | of section |lIB ll ~; 



• the model allows for first order transitions (this is not possible with the pure NJL model which 
only gives rise to second order transitions^ when /i„ = /is = 0). Physically this comes from the 
fact that quark and gluon condensates are not independent. Mathematically, this comes from the 
fact that we have to deal with coupled equations (68 69| , |8^ ). An example of first order transition 
is given in figure ^, section [I B 1 , It is clear on this picture that quark and gluon condensates are 
strongly coupled: the chiral phase transition is abrupt and felt by the gluon condensate; 

• in scaled models, a transition temperature as low as T « 140 MeV can only be attained if the 
vacuum gluon condensate is small. Then, there is the problem that the chiral symmetry can be 
restored at densities close to, or even lower than, the normal nuclear matter density. However, this 
defect can be related to the absence of vector mesons in our calculations. It is in fact well known 
that these vector mesons make the vacuum stiffer against the restoration of chiral symmetry with 
respect to density [^,^,^. We have however to check if these mesons do not change also too 
drastically the critical temperature. 



^^In figures ^ and we have kept a finite value for tlie current quark masses. This explains the long tail near Tc- 
The Tc value that we give is always the true value, obtained in the chiral limit so that the tail does not contribute 
(see the difference between figure ^ and figure This can readily be seen on the gap equation which gives, for 
a second order phase transition. 



The second term under the square root is due to the chiral limit which then reduces the critical temperature. 
This equation also shows clearly that the scaled model allows a reduction of the critical temperature compared 
to a pure NJL model, in the ratio Xc/xo- 

To keep things clear, we adopt the somewhat unconventional approach of calling a second order transition 
a transition which is continuous although, strictly speaking, we should have made the distinction between true 
second order and crossover. 
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F. Drawbacks of the models 



Effective Lagrangians are very attractive: they are quite simple while giving opportunities to have insights 
in an energy-momentum regime where lattice QCD is the only way to examine the true theory. There is 
a price for this simplicity. Physically, effective interactions are built in order to reproduce some effects at 
a given range of energy-momentum. This implies that their validity is restricted to this range. Even if 
calculations can be performed at another scale, there is nothing which enables to relate these calculations 
to the true theory. Mathematically, there are also several drawbacks which have to be handled if we want 
the model to be meaningful. In the NJL model, these drawbacks are the lack of renormalizability and 
the lack of confinement. 

• The lack of renormalizability comes from the point-like nature of the quark-quark interaction (see 
(^)). The model has to be regularized and this regularization cannot be removed by a renormaliza- 
tion procedure. Several regularization procedures are available (the references we indicate below are 
examples of application of these regularization procedure s to the NJL model) such as Pauli-Villars 



|35| , p2| , |9l[ , proper time ||35|, three dimensional cut-off ||81 , 100 , 101 1 , four dimensional cut-off 47 8§||, 
regularizators which depend only on the orthogonal component to the 4-momentum [ p9[ , ... An 
in-depth analysis of regularization procedures and more references to the corresponding literature 
is given in ||36| . Because there is no renormalization, it is clear that the regularization procedure 
is part of the model: to specify it requires to give both the Lagrangian and the regularizator. The 
hope is that observables are not too sensitive to a particular choice of the regularizator. For loga- 
rithmic divergences, this low sensitivity is quite understandable. However, this is a nontrivial task 
for quadratic divergences. However, fitting the regularizator to reproduce an observable such as the 
pion decay constant makes the results less sensitive to the particular choice. 

• The lack of confinement can be traced back to the fact that the quark propagator has a pole. A 
quark propagator without pole (with a p^-dependent mass) is not a necessary condition to get a 
confined quark but it is sufficient Jl9t since it ensures that the quark propagator cannot have any 
Lehmann representation. (Note, however, that the confined quark does not mean the absence of 
trouble to get meson properties. When calculating on-shell meson masses, we are looking for poles 
for time-like momenta < (in Euclidean space). The propagators have then to be continued 
analytically, which can introduce unphysical poles: the continuation is a rather hard task whose 
difficulty is related to the confining property. To quote Ripka |36| "Our ignorance as how to continue 
propagators in the complex plane reflects our ignorance of the confining mechanisrn!\) The lack of 
confinement of the NJL model means that quark loop s for phy sical meson processes will generate 



a qq production threshold at s = 4M^. Celenza et al. [102.103] still use the NJL quark propagator 
(constant quark mass). However they get loops without threshold by modifying the vertices in such 
a way that they cancel unphysical poles. Physically, these vertices are related to the confining part 
of potential models (this supposes that the confining part of a vertex can be separated from the full 
one, i.e. that it is possible to decouple the discussion of confinement from that of chiral-symmetry 
breaking). As a first step, we prefer to treat the NJL model as it is, trying to find a prescription 
to still give a meaning to it, even above threshold. We shall explain the prescription we use in 
section |H C| . Simply speaking, we can say that we use a phenomenological interaction where the 
confinement is introduced by throwing away, by hand, the unphysical part. This is justified as long 
as the unphysical decay width of the meson into free quarks is small compared to the mass. Since 
lattice calculations tend to show that the chiral phase transition and the confining transition are 
coincident |Q , we shall use the shortcut of speaking about the confined phase for the phase where 
the chiral symmetry is spontaneously, dynamically, broken, while we shall identify the unconfined 
phase with the phase where the chiral symmetry is restored (M„ — niu , Mg — nis ) . 

• Our model is an extension of the NJL model which takes into account the QCD trace or scale 
anomaly. Within the NJL model it is introduced in two steps [ [77|j79| ] : first a scalar x field is used to 
make the effective Lagrangian scale invariant (this includes also the treatment of the regularizator) 
and, second, an effective term is added which mimics the scale anomaly at tree level The 
vacuum is driven to a nonzero value of xo a-nd provides a mass for the dilaton excitations. The x field 
can be related to the trace of the energy- momentum tensor (|60| ) . One drawback of this approach 
is that it is assumed that the use of a single dilaton field is enough to saturate the anomaly, i.e. 
that the anomaly equation can be used to define one interpolating field for the giuball. Although 
working for the pion (PC AC), this idea of the interpolating field is however questionable for such a 
massive state as the glueball [Mgb ~ 1-5 GeV). The defect of our approach can be seen by looking 
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at the particle data book [104| where it is clear that the region near 1.5 GeV is populated by 
several scalar resonances. However, we take here the drastic assumption that all the scalars can be 
described with only three resonances 0: the /o(1370), /o(1500) and /j=o(1710). The robustness of 
this approximation has to be tested phenomenologically within the model. This will be performed 
when studying the scalar sector in section 

Without the anomaly, the scale invariance would be broken in two ways: i) explicitly, by the 
current quark masses and ii) spontaneously, leading to a Goldstone boson (massless dilaton field). 
However, the quantum breaking of the axial symmetry completely destroys this situation (massive 



dilaton field) and it can be argued that QCD does not scale 1105]. This is also shown in the model 



independent study [[75| and within the extended NJL model |63[. In ||76| it is explicitly shown in the 
extended NJL model that, to get a scaling of hadron masses and decay constants as suggested by 
Brown and Rho |68j, the dilaton mass has to be much smaller that the lightest (qq) scalar meson. 
Since this is not the case, the Brown-Rho scaling should not be expected in the data. Note that 
the model independent estimate made in | |75[ | concerns matter up to normal densities. It does not 
prevent the dilaton field to play a role at higher densities. In fact, our studies [ |63||8^ have shown 
that the gluon condensate only plays a role near the phase transition. This means that, although 
QCD does not scale (at least in the chirally broken phase), the gluon condensate can be relevant 
for phase transition studies. In our formalism, this can be traced back to the fact that the gap 
equations for the quark condensates are coupled with the gluon condensate (eqs. ( |68| , |69| , ^ ) . It 
is this coupling which allows our model to describe second order phase transitions as well as first 
order transitions. This also allows to get low critical temperature and densities. However, the 
most interesting cases (departure from the pure NJL case) show up only when xo is rather small, 
corresponding to bag models (the QCD sum rules estimate ||95| , p6| gives xo ~ 350 MeV, leading to 
an almost pure NJL behavior). 

We have given above arguments to claim that the gluon condensate is a meaningful quantity to 
describe the chiral transition, even if QCD does not scale. It is however clear that our approach 
could have a very strong shortcoming: we use the scale anomaly at zero temperature. When there 
is no quark in the theory {mu = rud — rug — *■ oo), eq. ( ^2| ) gives y — Xo for E^ny temperature or 
density^ However, the gluon condensate is expected to be flatQ below the transition, so that 
its temperature and density variations can be neglected, as a first approximation, in the chirally 
broken phase. 



^^In view of this, all the observed scalars near 1.5 GeV are manifestations of the same particle: the /o(1500). 
In the NJL model, the gluons have been eliminated to the benefit of an effective 4-quark interaction. In our 
scaled model, some effects of the gluons have been included (scale anomaly). However we still do not have a gluon 
sector. It is then not astonishing that removing the quarks does not lead to a gluon phenomenology: for example 
we do not have the ri ght n umber of gluonic d.o.f. 

*^See e.g. eq. (3) of |106] or lattice studies in the pure gluonic sector (e.g. figure 4 from [10'i|). 
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II. THE A-SCALED NJL MODEL AT FINITE TEMPERATURE AND DENSITY 



In the previous section, we have carefully described the motivations leading to the A-scaled NJL model and 
we have commented on the drawbacks of the approach. This section gathers results at finite temperature 
and density. We shall first focus on analytical results for the thermodynamics in section [I A. We shall 
then show our numer ical r esults concerning the condensates, the thermodynamics and the critical surface 
Tc(Mu,Ms) in section [IB . We shall end the section with a study of the anomalous decay of the neutral 
pion into two photons. This will give us the opportunity to extend the model above the threshold. 



A. Thermodynamics 



It is believed that QCD has two phase transitions: a deconfinement transition corresponding to going from 
a hadronic gas to a quark-gluon plasma, and a transition leading to a phase where the chiral symmetry 
is restored. Lattice calculations |g7,107| suggest that these two transitions coincide. In a purely gluonic 
theory, it seems that the transition is of first order (for Nc = 3). When quarks are introduced, the order 
of the transition depends on the number of light flavors. With Nf = 3 massless flavors, QCD has a first 
order chiral transition, not connected to the pure gluonic one. When the mass of the quarks is varied, the 
two first order transitions are separated by a region in which there is a crossover (to keep things simple, 
we shall call it a second order transition, however). According to the type of calculations - Wilson 
fermions [108| or staggered (Kogut-Susskind) fermions [109| - QCD (m„ « md ~ 10 MeV, iris ~ 200 
MeV) appears to be in the first order region, or in the second order region, respectively. It is then not so 
clear which case occurs and we have the freedom to play with the parameters in such a way as to allow for 
both types of transitions. Note that the distinction between pure gauge deconfinement and light quark 
chiral phase transition is of prime importance: the energy scales are different. Pure gauge transitions 
occur at a te mperature of abo ut 260 MeV while, with two light flavors, t he critica l temperature is around 
150 MeV [p7| JlO^JlO^ , pl0| , pTl[ , or even as low as 140 MeV according to p^ , pl3| |. 

The deconfining transition is hard to study in existing models. However, chiral symmetry breaking and 
its restoration can be studied in effective models such as NJL. Since both transitions are seen to be 
coincident on the lattice, we shall take the point of view that the study of the chiral phase transition 
can shed light on the deconfining transition. To simplify the discussion we shall call hadronic phase the 
phase where chiral symmetry is spontaneously broken while we shall call quark-gluon plasma phase the 
phase where chiral symmetry is restored!^ We shall also take the point of vie w tha t gluon confinement is 
linked to the gluon condensate Xs, although this hypothesis is questionable jl^,114]. There is a constraint 
between the quark and gluon condensates, which shows up into the form of three coupled equations for 
the three condensates. We can then get, according to the strength of this constraint, a second order 
transition (continuous transition from one phase to the other, see figure |^) or a first order transition 
(discontinuous passage from one phase to the other, see figure 

As already stated we work 
and Nikolov et al. l8^] go 



The preliminary results given in this section have been discussed in 
at the mean field level. The groups of Rostock , Heildelberg 

beyond this approximation, studying the first 1/Nc corrections and showing that they are not negligible 
at low temperature and density (because pions are almost massless) . However one can also take the point 
of view that an effective theory cannot be used above the mean field so that corrections should not be 
included. As a first approximation, we shall keep this argument. 



1. Pressure, energy density, entropy density, bag constant 



In this section we want to understand the equilibrium properties of a strongly interacting matter, i.e. 
determine the relations associated to the thermodynamics of a hot and dense system. The system is 
modelized by an effective action of the NJL type (free massive constituent quarks) with a dilaton field 
included^, eq. ( |6l| ) . In a grand-canonical system, the partition function is given by 



^^Strictly speaking, only the quark condensate goes to zero in our model. However we include the gluons in this 



plasma phase when the gluon condensate reaches its minimum value. 
*^For the thermodynamics of a scaled linear sigma model, see [ [t^ . 
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Z = exp(-/3f2), 



(110) 



where is the thermodynamical potential (or grand potential) . We exclusively consider a system in equi- 
librium: all the descriptions (micro-canonical, canonical, grand-canonical) are equivalent. However the 
grand-canonical description is the easiest. In the canonical formalism, the basic quantity is the Helmholtz 
free energy and the independent variables are the temperature, the pressure and the densities (one for 
each chemical potential). Since phase transitions occur at a constant chemical potential, and not at a 
constant density, it has always to be checked if a lower energy solution, obtained by separating the system 
into subsystems, exists [116|. Working in the grand-canonical formalism, where the independent variables 
are the temperature, the pressure and the chemical potentials, there is a d irect access to the solutions 
corresponding to the minimum of the thermodynamical potential. Eq. (110) leads to the identification 



Since 



we get 



116| 



d£l = —SdT — PdQ — pid^i [i — u, d, s). 



P 



an 



This implies 



iain(£)1 

/3 dn J^,^^ 



(111) 



(112) 



(113) 



(114) 



Physically, the pressure is not an absolute quantity. We have to consider it with respect to a reference 
system which is chosen to be the (non-perturbative) vacuum, of pressure Pq. Defining P — Pq = P' and 
replacing P by P, we then have 



P 



IdHZjZo 

/3 on 



(115) 



As explained for eq. ([73|), we can separate the action into a Fermi part and a Dirac part. Subtracting 
the vacuum, the lowest order of the action is 



11^^^ {M^,M,, Xs) + /fo,oo) {Mu, Ms,Xs) - Ilo,oo) {Ml xo)- 
The first term is the Fermi part, which does not need being regularized, 

i—u,s ^ 

X <; In (l + exp [-/? (Ei + fi,)] ) + In f 1 + exp {E, - fi. 



(116) 



where0 Ei = ^k^ + Mf, and a,, 
and is given by 



(117) 



2,as ~ 1, while the Dirac part (T = /i = 0) has to be regularized 



^When no confusion is possible between the 4-momentum k and the 3-momentum k, we use the notation k for 
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2iVc 
327r2 



Mr 



i(AXs)^ - ^(Axo)') + ((Af,Ax.)2 - (MfAxo)') 



E 



With the action (116), eq. (115) gives 



P 



ldIlff{Mu,Ms,Xs) 
^13 dVL 



which leads to 
P = 



1 



I^^P) (Mu, Ms,Xs)+ /fo.oo) {Mu, Ms, Xs) ~ 11^^^) {Ml, Ml xo) 



The physical meaning of this equation is that the grand potential is an extensive quantity: 

n = -pn. 



(118) 

(119) 

(120) 
(121) 



In section [D3, we have seen that, mathematically, the mean field approximation corresponds to fi nding 
the saddle-point solution of the equations (|63|,|65|), i.e. to finding the minimum of the action. Eq. (121 



shows the physics attached to this condition: the system chooses the phase where the pressure is a 
maximum. 

Quark densities can be evaluated from (112) and are given by (see ( [79| ) for the definition of the functions 
«j±) 



on' 



1 fdl. 



'eff 



I dk. 



(122) 



Since we work in the isospin limit (to„ — m^j) for a sy mmetric matter (/i„ = /z^), it is clear that p„ 
The entropy density can also be evaluated from (112) and is given by 



S 

n 



(M„,M.,x.). 



(123) 



Only the Fermi part appears in this formula since this is the only one which depends upon temperature. 
When T ^ the entropy density s goes to zero , in agreement with the third principle of thermodynamics. 
Finally, the internal energy is given by 11£ ] 



d 



d 



E=n + TS + ^i,p^= 1 + /37J73-M.T^P^= [l^- 15^^^1^]I!ffiMu,M„Xs) 



d(3 "^'d^i 



d 



d 



(124) 



This gives the energy density 
E 1 



1 

pn 



'(0,00 



) (M„, M„ X.) - /fo,oo) {K, Ml xo 



(125) 



Like the pressure, the energy density is a relative quantity : ( |125| ) gives the density energy of the system 
w.r.t. the vacuum energy. 
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2. Bag constant B 



Following [113, 117 1, we write 



P = ^Idcal gas — ^ (/3, ^t, , /i, 
e = Eideal gas + B (/3, , M. 



s] 1 



Ts — ^'idcal gas ^idcal gas l^iPii 



with 



-Pi deal gas 



(3n 



— {ni+ + Ui^) dk, 



^idcal gas 



1/9 1 d 



n \dl3 13"^ dfi. 



7r2 



(126) 
(127) 
(128) 



(129) 



(130) 



being quantities relative to a massive free quark system. The interaction measure, which is an indication 
of non-perturbative effects, is 



Nc ^ , r 



E, 



(131) 



In eqs. ([126D and (127), we have defined a temperatu re an d density dependent (through M„, M^, Xs) bag 



(132) 



constany'l i?. It depends on the Dirac contribution ( |118| ) to the pressure: 

B (/?, = ^ {/fo^^) [Mu, Ms, Xs) - /fo.oo) (^", A'/,", xo) } 



The definition (132) is different from | |63| , |70| , pl ,118| because two new effects are implemented: i) we take 
into account the explicit breaking of chiral symmetry (even if = 0, we have TOs ^ 0: the strange quark 
contribution can be non - negligible) ; ii) as already stated, eq. (1 132^ t akes into account the effects of the 
gluons. Moreover, ( |132| ) conceptually differs from the definition ||8l|] where the bag constant is zero in 
the chirally restored phase. Finally, it is different from the bag constant B' introduced in [E2l^5|j7^. In 



these references, it is only obtained at zero temperature and density, through the definition of B , being 
identical to the energy difference between the perturbative vacuum and the true vacuu mQ. 
In the perturbative vacuum, the chiral symmetry is restored and the gluon condensate vanishes. Then, 
according to 16^, 



B' 



^ {no,oo) (0, 0, 0) - /fo^^) {Ml Ml xo) } = Yq'^IlxI 



f3n 



(133) 



In SU(3), implementing uis = in the definition of the perturbative vacuum makes no sense. However, 
we sha ll see that the value of B'^^"^ is not so far from that of B^^'^ (132) (for T >>) so that the use 
of (|l33| ) in [^,^ 77 1 is verified a posteriori. This remark also applies for the effect of Xs which does not 
go down to ze ro. 

Note that in |118|, the authors study both the chiral symmetry restoration and the effects of the gluon 
condensate. They however define two bag constants, one associated to the restoration of scale symmetry, 
the other to the restoration of chiral symmetry. 



In 1 113 117 1, there is no glueball: the bag constant is purely chiral. Our bag constant is then a generalization 
of these references. 

''^One can also take the equivalent definition of considering it only through the glueball Lagrangian (decoupling 

1 4 1 2 2 

between the glueball and the other fields). We have B = - < >= — x = —moLX ■ 



16' 



16 
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3. High temperature zero density limit (T > Tc) 



In a phase where chiral symmetry is restored - in this section, we mean the phase where the constituent 
quark mass goes to the current quark mass {Mi = rrii, i — u, s) even if we arc not in the chiral limit - 
we have nis/T < 1: a high temperature expansion in ttLs/T is possible [115,119,120|. Calculations are 
lengthy and left to appendix taken from [p8| . We work at zero density. Our results are a generalization 
of 115, 119| , 120| ] where only a limited number of terms in the mg/T expansion have been retained while 
we are able to give here the full expansion, involving only elementary functions. Note that to describe 
the results in section II B, the first four terms will be enough (T > Tc, with Tc the critical temperature) 
so that we only keep them in the following. For the pressure, we get 



^idcal gas ~ J" ( gQ^^TT 



Nc ml 
12 



Stt^ T4 ttT 



3 „ \ m"! 



(134) 



where 7 is the Euler constant^. 

Note that the case of finite chemical potentials is much more complex and, to our knowledge, has never 
been treated to all orders in the fermionic case. (In the bosonic case, the constraint < Mi ( not presen t 
in the fermionic case) allows a high temperature expansion {mi/T and /x^/T ^ 1) to all orders |121| - |123| .) 
Above Tc, the bag constant B is temperature independent and writes 

1 



(0, ms,Xc 



(0,00) 



)} 



(135) 



where Xc is the gluon condensate above Tc- Note that, for any set of parameters (M°,xo)j we could not 
get Xc — 0. In our model there is never a complete gluon deconfinement. This is related to the fact 
that gluons are only poorly incorporated in our formalism (we have no explicit temperature dependence 
of the gluon condensate since the modeling of the gluon anomaly is through a temperature independent 
potential (|5^); we also do not h ave t he right number of gluonic d.o.f.). 
The energy density is given by (125) 

1 



^idcal gas ^ 



n dp 



(136) 



i.e., using (|lll| ) and ( |l2l| ) 
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gas 
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ideal gas ) 



With (134), the high temperature, zero density, expansion gives (appendix ^) 



•^idcal gas 
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Finally, (111) and (121) applied to (123) lead to 



ideal gas ^ 



or, equivalently, to 



2 9 d 

^ — P ideal gas ^ideal gas- 



The high temperature, zero density, expansion of the entropy density is then (appendix ^) 



15 



6 r2 



Nc m\ 



(137) 



(138) 



(139) 



(140) 



(141) 



The resuhs (134), (138) and (141) are used in section II B 



In |8l||, the massless free quark limit is unreachable, by construction: the d k regularization introduces a cut-off 
A for each Fermi or ideal gas quantity. These quantities, for example -Pjdeal gasi behave then at high temperature 
as A/r, decreasing to zero. This drawback is not present in this work, where we have chosen a d^fc regularization 
for the vacuum while the Fermi part is not regularized. 
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4- Low temperature zero density limit 



To simplify the discussion, we limit ourselves to the zero density case. If is of the order of 400 
MeV, the low temperature expansion (see appendix ^) is valid up to T « 100 MeV. Indeed, masses 



and condensates are not varying within this range of temperatures, see section [IB, and the expansion 
parameters /3M„ an d 0M s are then large enough - we have at worst PM^ « 4 - to allow a stationary 
phase expansion of (|l29|): 



-Pi 



ideal gas 



47V,/3- 



-5/2 



(27r 



,3/2 



3/2g-;3M. 



corrections. 



(142) 



where a^, = 2, = 1. To get the correction^^ it is better to work with ( |l2g| ) written in terms of K2, 
of which the asymptotic behavior is well known. The method is explained in appendix ^ which also 
contains the low temperature zero density expansion of the energy and entropy densities. 
The results relative to this section are given in section [I B 3 . They necessitate the knowledge of the 
behavior of the condensates as a function of temperature. 



B. Results: condensates and thermodynamical functions at finite temperature and density 



We show the variation as a function of temperature and density of the quark and gluon condensates in 
The results concerning the critical surf ace T{ nu, are given in section II B 2. We discuss 



section [I B 1 



the results relative to the thermodynamics in section II B 3 



1. Condensates 



The quark and gluon condensates are given by eqs. (|80|-82). The constituent quark masses M„ and Ms 
used in these equations come from the gap equations (|68|,69). When m„ — 0, = is always a solution 
of (|68|), so that it has to be checked if it corresponds to a greater pressure. We define a second order 
transition to be a transition for which the slope of the pressure as a function of the external parameters 
T, /z is continuous (this includes both true second order transitions and crossovers) ; otherwise it is said 
to be of the first order. We show results for four sets of parameters 



(i) M° = 450 MeV, xo = 350 MeV, 
{ii) Ml ^ 600 MeV, xo = 125 MeV, 
{iii) Ml = 300 MeV, xo = 80 MeV, 
{iv) Ml = 800 MeV, xo = 90 MeV, 



(143) 
(144) 
(145) 
(146) 



for which the model parameters have been obtained in table |, section [E. Different behaviors show 
up according to the chosen set and to the external parameters: we can have two first order transitions, 
coincident or not, or two second order transitions. We can also have a second order transition for one 
species of quark while the other experiences a first order transition. A discontinuous slope for the pressure 
means a discontinuous condensate: there is a mass gap. This is illustrated in figure ||. 
Since the pressure at point A is identical to the pressure at point D, the transition takes place between 
these two points. Only the region from B to C is unstable: there is metastability between A and B and 
between C and D. 

Note that looking graphically at the pressure to find the transition point is manageable only in the two 
degenerate fiavor model at zero density or temperature. In that case, the strategy is the following. Given 
M„, the corresponding Xs can be extracted from eq. (|8^) with the strange quark contribution removed 
(two-fiavor case). We can then use this couple {M^tXs) to extract the corresponding temperature from 
the gap equation, eq. (|68|). There is only one solution. (Should we have fixed the temperature, we should 
have to solve two coupled equations (|68|,p2|) with the problem that, for a first order transition, several 



Corrections to (142) are negligible only if /3Mi > 40. However their number is limited for /3Mi as low as 4. In 
appendix we quantitatively discuss the importance of these corrections with respect to the value of (3 Mi. 



34 



solutions are possible.) We then get the ABCD curve of figure ^(a). To get the transition point it is then 
enough to look at the pressure, figure §(b). 
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FIG. 3. Up constituent quark mass transition (with ruu 7^ 0^. 



Of course, as soon as more variables (/iu, /is, Ms) are introduced, this strategy is not anymore of interest. 
We have to solve the full set of coupled equations 69 18^) as a function of chemical potentials and 
temperature. Because of possible first order transitions, several local extrema of the pressure show up 
so that it is necessary to use a numerical algorithm searching for g lobal extrema. In our work ]115| ], we 
used the simulated annealing algorithm that we adapted from |124]. 

Figures H an d || are for the sets ( |143| ) and (145), respectively, and correspond to a second order phase 
transitioq |w.r.t. the temperature for vanishing chemical potentials. They are similar to results obtained 
in the two-flavor case p^,|65[|. The choice (143) almost corresponds to a pure NJL model, w hile the choice 



|145|) leads to a lower critical temperature, in agreement with lattice results [107 111-115] 



If we take the chiral limit analogue of figure ^, the critical temperature is found to be Tc ~ 193 MeV. 
The gluon condensate is almost flat, so that the quark and gluon condensates are almost uncoupled. The 
scaled model is then essentially a pure NJL. In figure ^ this coupling is more important and the chiral 
limit gives Tc « 140 MeV. Working in the three-flavor version of the model, these two pictures show also 
the strange quark condensate, for which we can make two remarks: 



< ss > decreases slower than < uu >, in agreement with [125|. This can be traced back to the 



See figure El section IE, for the order of the transition versus the parameters (M„, Xs 
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greater constituent strange quark mass compared to the up quark (see table^ which, in turn, is a 
consequence of a greater current strange quark mass compared to the upone|_| Note however that 
our results show a faster decrease of < ss > with temperature than in [126|; 



because of the coupHng between the condensates, the gluon condensate at the transition, Xc^ is 
smaller than in the two-flavor case. For the set ( |145D , we have Xc{SU (2)) w 0.8 while figure || shows 
Xc{SU{?,)) « 0.6. 
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FIG. 4. Quark and gluon condensates as a function of temperature for the parameters M„ — 450 MeV, xo = 350 
MeV ((ti pure NJL model). 
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FIG. 5. Quark and gluon condensates as a function of temperature for the parameters M„ = 300 MeV, xo = 
MeV. 



Figures H and ^ are the analogues of figures ^ and |5| for the sets of parameters (144) and (146), respectively. 
These two sets allow to get a first order phase transition |62 7^. With the second one, however, the critical 
temperature is lowered down to Tc ~ 140 MeV, in agreement with lattice QCD |107,111,113]. 



With these two sets of parameters, the coupling between quark and gluon condensates is so strong that 
all the condensates undergo the transition together. Although these figures suggest that a first order 



Note that, in the two-flavor limit where Ms, rUs oo, we have < ss > / < ss >o~ 1- 
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transition is connected to a low value of the gluon condensate, both and xo really determine it, as 
shown in figure ^ and confirmed in figure ^ 
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FIG. 6. Quark and gluon condensates as a function of temperature for the parameters M„ — 600 MeV, xo = 125 
MeV. 
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At zero density, the above analysis shows that: 

• we can reproduce both first and second order phase transitions. A first order transition is typically 
an effect due to the gluon condensate which then does not show up in pure NJL models; 

• a low critical temperature, as low as 140 MeV, can be reproduced. This is clearly related to the 
coupling between quark and gluon condensates. Pure NJL models are then unable to reproduce 
such a low temperature: they cannot go below Tc ~ 190 MeV. According to eq. ( |109| ), scaled 
models allow (for second order transitions) a reduction of the critical temperature in the ratio 
Xc/xo compared to a pure NJL model; 

• the coupling between quark and gluon condensates is mainly driven by the value of the vacuum 
gluon condensate xo- With large xqj the coupling is weak while the coupling becomes more and 
more important as we decrease xo- The quark or gluon condensate can then be considered as the 
order parameter for the phase transition. 
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It should however not be forgotten that all the above analysis is performed without vector mesons. 
Without them, our results show that high values of the gluon condensate are needed to get a transition 
above the normal nuclear matter density po- This is illustrated for two flavors in figure @. 
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FIG. 8. Quark and gluon condensates for the parameters M„ — 300 MeV, Xo = 350 MeV for niu = (plain 
line) and rriu 7^ (dashed line) as a function of the chemical potential (a); as a function of the density (b); (c) 
shows pu versus pu- 
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For = 300 MeV, the smallest required value of the vacuum gluon condensate, for a second order 
transition, is xo = 350 MeV. Even in this case the chiral limit leads to a first order transition. It is 
however clear from figure |8|(a) that, although sharp, the transition is of second order if niu is different 
from zero. Figure ||(b) is also instructive is the sense that it shows clearly that looking at the density is 
misleading to have insights on the transition. The dashed curve is a second order transition while the 
plain curve is a first order transition. The transition occurs from below the normal nuclear density to 
above the normal nuclear density (the phase between the two vertical bars is unstable). The fact that a 
small current quark mass can stabilize the system has already been noticed in [127| in a pure NJL. Since 
a pure NJL does not lead to a first order transition at zero density finite temperature, this indicates that 
the temperature has a stabilizing effect on the transition. 

Although the scaled models are not able to reproduce a transition above the normal nuclear density for 
a small value of the gluon condensate (which leads to a low critical temperature), it does not mean that 
they are inefficient. Indeed, it is well known that vector mesons make the vacuum stiffer against the 
restoration of chiral symmetry |^,|9j. Including these mesons will then correct what seems to be, at first 
sight, a drawback of the model. 



2. Critical surface Tc{iJ.u, iJ.s) 
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FIG. 9. Quark and gluon condensates as a function of temperature for the parameters = 600 MeV, xo = 125 
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MeV (b). 
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In this section, we add to the results of the previous section the effect of the density throug h the chemical 
potentials /i„ = Hd (symmetric matter) and ^s- We base our discussion on the work | 115 | which is 
a systematic study of the chiral phase transition, for three flavors, as a function of temperature and 
chemical potentials. To eliminate the long tail due to finite current quark masses, we work in the chiral 
limit for the light quarks (to„ — — 0) while we keep nis ^ which is the only source (apart from the 
chemical potentials) of chiral symmetry breaking. We also restrict ourselves to the sets of parameters 
(144, 145| ) which yield, at zero density, a low transition temperature of respectively first or second order. 
Adding chemical potentials to the temperature study is interesting because it shows new behaviors: 
separated transitions for up and strange quarks, both of first order or first and second order. This is 
illustrated in figures ^ and 1^. 

Figure ^(a) is quite similar to figure ^: both transitions are of first order and coincident. However, the 
addition of a chemical potential shows that the critical temperature is reduced, and eventually vanishes, 
if the density is high enough. This is illustrated in figure |l2|. Figure ^(b) gives a new behavior: the 
up and strange transitions are disconnected except through the gluon condensate x- In fact, when the 
up quarks make their transitions, eq. (js^ ) for the x field shows a discontinuity which also implies that 
the strange quarks feel the up ones making the transition. This is what happens around T w 140 MeV. 
However this is not the true strange quark transition: at T w 160 MeV, the strange quarks perform their 
own first order transition. 
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FIG. 10. Quark and gluon condensates as a function of temperature for the parameters A/„ — 300 MeV, xo ~ 80 
MeV in the light quark chiral limit ruu = rrid = for /Xu = /Xd = Ms = MeV (a); for fJ.u = fJ^d = fJ-s ~ 250 MeV (b). 
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Figure the analogue of figure]^ for vanishing light quark masses. It shows a second order transition 

for the up condensate while the strange condensate feels a weak first order transition. 
It is interesting to notice that the tail in figure || has completely disappeared, showing that the transition 
for the up quarks (for vanishing density) occurs at about 150 MeV (140 MeV, would we have only two 
flavors). It also shows that a small up current quark mass has a stabilizing effect for the strange transition. 
This is traced back to the mixing of the up and strange quark condensates through the gluon condensate. 
Figure |l^(b) is for the same set of parameters as figure |l^(a) but for nonvanishing chemical potentials. 
This explains why the transition temperatures are reduced. Note that we have two noncoincident first 
order phase transitions. Because the coupling between the condensates is smaller than in the case shown 
in figure |[ the first transition of the strange quarks, connected to the transition of the up quarks, leads 
only to a small mass gap. Its true transition occurs around T « 150 MeV. 

As we have seen, up and strange transitions can be completely disconnected. This is also illustrated in 
figures ^ and |l^ for the sets of parameters ( [144 ,145), respectively. Open circles represent the phase 
transition jigc as a function of /i^p. For figure ]ll[7 they correspond to a vanishing up quark condensate 
< uu > at a temperature (figure jll|) Tc = MeV (a), Tc = 100 MeV (b) and Tc = 161 MeV (c). The 
transition associated with the strange quark condensate < ss > is given by the crosses. 
We also observe on figure |ll| that, at low temperature, up and strange transitions are independent as 
in a pure NJL without 't Hooft determinant, except in a small region of the plane (/iscMfs) where the 
coupling with the gluon condensate is high enough to constrain the two condensates to vanish coincidently. 
This region is the only one to survive when the temperature is increased. At any point {Tc, fj-sc, fJ-us), 
the transition is of the first order, which is well illustrated on figure |^ in the case of a strong coupling 
between the condensates (figure ^(a): w 120 MeV, ^^c = Msc = 250 MeV) and in the case of a weaker 
couphng (figure |(b): T^c = 140 MeV, Tsc = 160 MeV, ^ 200 MeV, ^lsc = MeV). 
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FIG. 11. Chiral phase transition for the parameters M„ — 600 MeV, xo ~ 125 MeV in the light quark chiral 
limit mu = md = 0. ^lsc as a function of ^^c forT = Q MeV (a); for T = 100 MeV (b); for T = 161 MeV (c). 



Note that, in our model, a discontinuity in the behavior of one of the quark condensates imphes a 
discontinuity of the gluon condensate, in contradiction with the resuhs^ of |114]. Note also that, above 
the transition, Xs ^ 0: gluons are never completely deconfined. 
Figure O sums up the preceding results in the form of an "artistic" view. 



T 

c 




^ The authors of this reference also make the distinction between the two scales of deconfinement and restoration 
of chiral symmetry, arguing that valu es of the MIT bag constant has more to do with chiral symmetry than with 
quark confinement (see also [113 118|). 
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The space is divided into four regions: the internal region corresponds to quarks in the condensed phase 
or, in our terminology, in the hadronic phase; the external region is interpreted as the quark-gluon plasma 
phase; the last two regions, extending along the axis up to infinity, correspond to mixed phases with free 
up (strange) quarks in presence of strange (up) quarks in the hadronic phase. We notice an interesting 
prediction of the model: if non-strange matter (with a high enough up chemical potential) is heated, we 
observe the appearance of light free quarks in thermodynamical equilibrium with strange hadrons (of 
zero net strangeness, because /i^ = 0). These strange hadrons then transform into a free strange quark 
gas if the temperature is again increased. 

As already mentioned, a drawback of the model shows up when densities are calculated. When the 
transition is of first order, there is a density jump at the transitiorj^. The problem is that the jump 
occurs from before to after the normal nuclear density: in our terminology, the hadronic phase would 
be unstable. For example, at //„c — ^idc = tJ-sc = 250 MeV, we have p„ before/po — 0.362, ps beforo/po 
— 0.050 and p„ after/po = 2.646, ps aftor/po — 0.978, where po is the normal nuclear density of quarks 
(0.51/fm3). To cure this defect, one has to take a large value for the vacuum gluon condensate, while 
keeping a small value for the constituent quark mass in the vacuum M'^, see figure ^ This means that 
the model behaves as a pure NJL. In our model, scalar meson masses would then be badly reproduced 
(roughly speaking, Mmeson ~ 2Mu) a nd the critic al temperature would be large (Tc w 200 MeV) 
contradiction with lattice calculations |107, 



111 



, 113| |. However, as already indicated, one expects ||7|_ 
that the introduction of the vector mesons into the formalism will solve this problem, making the vacuum 
stiffer against the restoration of chiral symmetry. This would then allow to get larger densities before the 
restoration takes place, while maintaining intact the general behavior of the phase transition, since it is 
mainly related to the chiral mesons and their coupling to the condensates. 



The preceding results remain qualitatively unchanged when the set of parameters (145) is studied, except 
for the order of the transition: at low density, the transition is of second order (see figure p^(a)). 
Figures 13 



Figures 
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and n are the equivalent of figures |l] and |T^, respectively. 

and ^3| are qualitatively identical, except at high tem perature wher e th e quark condensates 
are always coupled, for any choice of (/x„c,Psc), for the set (144). For the set (145), the up quarks can 
only exist in the deconfined phase with a restoration of chiral symmetry for strange quarks at fj,sc ~ est. 
The results from figure O can be summed up in a 3-dimensional space, of which we give an "artistic" 



view in figure 14 



The results from this section form the basis to determine those on thermodynamics, which are given in 
the next section. 
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FIG. 13. Chiral phase transition for the parameters = 300 MeV, xo = 80 MeV in the light quark chiral 
limit mu = mci = 0. ^sc as function of ^uc forT = MeV (a); for T = 100 MeV (b); for T = 130 MeV (c); for 
T = 130 MeV (d). 
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FIG. 14. Chiral phase transition summing up results from figure for the set of parameters 



3. Thermodynamics 



In this section, we present results relative to the pres sure (equation of state), the energy density and the 
entropy density that we have adapted from [^|ll5 |. Because we would like to emphasize some points 
relative to fits, we first take the somehow unusual way to present these quantities as a function of T"* 
(pressure, energy density) or (entropy density). This allows us to separate curves corresponding to 
different parameters. Once these results will have been presented, we shall redraw some of our results 
in the usual way (pressure or energy density or T times the entropy density, over T*, as a function of 
T) , allowing us to make a more direct comparison with the general shape of these qua nti ties as obt aine d 
in lattice calculations. Pressure, energy and entropy densities are obtained from (120), (125) and (123), 
respective ly. The general behavior of these quantities can be understood (at vanishing densit y) from eqs. 
( 126| -128) with asymptotic behaviors given in appendix ^ and summarized in section II A 1. 



Pressure 



The behavior of the pressure is shown in figure |l^(a) for a vanishing density. As the strange quark mass 
iTis is different fro m ze ro, the behavior of the pressure above Tc is not the usual l aw: in addition to 
the bag constant ( |135| ) which is taken into account through the decomposition ( |126| ), the massive free 
gas part has the temperature expansion (134). 
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FIG. 15. Pressure as a function of for the set of parameters (144,145) for /i„ = /id = Ms = MeV (a); for 
fj.n = fJ.d = fJ.s = 250 MeV (b). 

Figure p^a) shows also the hnear fit in which is valid in the phase where the chiral symmetry is 
restored: 



for the set (144), and 



P w 3.40 - 2.097 10"^ 



P« 3.40 - 1.195 10" 



(147) 



(148) 



for the set (145). Note tha t the coefficient of T'^ is not equal to the coefffcient INcTt'^ /60 (=3.454) of the 



term T in the expansion ( 134 ) . Because of the limite d ra nge of temperatures investigated in figure |lf 



the corre ctive terms to the factor in the expansion (134) change the apparent slope, and the constant 
term in ( |l47|) cannot be identified with the bag constant B. When a fit is realized with all the terms of 
the expansion ( |l34| )p^, plus a constant C to adjust, this one equals C^^^ « 207 MeV for the set (144) and 
C^/^ w 179 MeV for the set ( jlisl ). These values are close to B^/^ « 209 MeV (set (pi) and S^/"* w 
183 MeV (set (|l45| )) obtained from a direct calculation of the exact bag constant (135). This shows 



the consistency of our numerical results and the fast convergence of the expansion (134), even if the 
expanding parameter nis/T is not so small! Note also that the numerical values of B^''^ extracted from 
the exact bag constant (135) or extracted from the fits are not so far from B'^^'^ given by the approximate 
equation (133). The latter is valid if we ignore the coupling between the quark and gluon condensates, 
and gives B''^^'^ = 201 MeV (set (pl^) and B'^^'^ = 161 MeV (set (|l4|y). Thisjustifies a posteriori the 
introduction of the approximate bag constant (132) in the references ||62| , p3|j7^ ]p^. 

This discussion shows that it could be quite dangerous to extract n ume rical v alues from fits: one could 
be tempted to identify the bag constant from the constant term in (147) and (148). The above analysis 
shows clearly that it would be wrong. 

Beyond the transition, the behavior of the pressure versus temperature depends on the order of the 
transition. For the set (144), /3AIi [i = u, s) remains quite large, so that the T behavior is described in a 
first approximation by (142). However, a careful analysis shows that this expansion is not well suited for 
(3Ms < 20, or even for (3Ms < 40, and that the expa nsion ( D35 ) should be used instead. The exponential 



behavior e ^ ' is however still correct. For the set (145), the phas e tra nsition is of second order so that 



Mu is progressively decreasing. This implies that the expansion ( D35| ) cannot be used over the whole 



^^By all the terms we mean the terms which are inc luded in (134). 
not have to take into account more terms defined in ( D32 ). 

56t 



Indeed the precision is sufficient and we do 



"It should however be stressed that it is better to work with the exact expression since an error in B 
amplified when going to B. 



1/4 
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83 



range T <Tc a nd, e ven more, that it is too crude to describe this behavior. We have however found that 
the expansion ( D34 ), with the sum hmited to n = 1 and n = 2, is well suited for temperatures below 100 
MeV. 

Note that for an analysis based on the 1/Nc expansion 



, it has been shown that the pious 

give the largest contribution to the thcrmodynamical quantities at low temperature^. In the chiral 
limit where the pions are massless, their behavior is in T"*, which effectively shows they have a bigger 
contribution than (142) based on the high constituent quark mass. This is an example where the l/Nc 
approach to the lowest order is not valid (see |]36|| ). 

The above analysis shows that an important ingredient is not included: indeed, the obtained for 
^ideai gas docs not coutaiu the gluonic d.o.f., as indicated in section [F. In the chiral limit, we should 
have 



Pidcal gas = ^^27V, + l7V,iVy4^T- 



5.2 



(149) 



with iVg = 8 if iVc = 3. 

Although the gluon condensate is in part due to these gluonic d.o.f., gluons do not contribute to the 
thermodynamics. To take into account the thermodynamics of a purely gluonic system, we should add 
to the Lag rangian ( ^8|) a temperature dependent potential V^{T). This has for exam ple b een noticed in 
| |70| , [76| , p"30| . The choice of this potential should be such that it leads to the behavior (14£). 
Figure]T5|(b) is the analogue of figure |l^(a) for the choice = fJ-s — 250 Me V. T hese values correspond to 
a strong co uplin g between the condensates < uu > and < ss > for the set ( |144| ) , and to a weak coupling 
for the set ( |145| ), as shown in figures |9|(a) and [lO|(b). With a chemical potential, there is no poss ible 
linear fit. This can be seen considering the chiral limit (m„ = md = nis = 0) of ( |12(S| ) with (129). This 
leads to 



60 2 ^ 47r2 



(150) 



Because ^ 0, eg. (15C) has to be modified. Figure ^(b) shows that the part is not modified 
compared to (150), while the part is slightly smaller. Once again, this shows that one has to be 
very careful when makin g fit s. The coefficient of the term is not identical to the one of (147): the 
supplementary terms in ( 134 ) , and terms coming from the chemical potential, have a repercussion upon 
all the terms of the fit. 



Energy 



The behavior of the energy density versus temperature, as given by (127) and (13C), is represented in 
figure |l^. At high temperature and for vanishing density, the whole expansion ( D37 ) can be restricted 
to the first four terms (138), which perfectly describe the curves above the chiral transition. For the sets 
(|14|) and (|T4|), we get the linear fits 



and 



10.269 + 1.620 10" 



10.271 + 0.877 10" 



(151) 



(152) 



respectively. 

It is useful to stress once more the difficulties one encounters to extract meaningful information form fits 
such as ( |151[ ) and (152), si nce t he co nstan t terms in ( |151| ) and ( |152| ) are not the opposite of (147) and 
( |148| ) , while the exact eqs. (126) and (127) show clearly that they should. In fact, all the terms from the 
expansion ( |138| ) contribute to the determination of the coefficients of the fits in the restricted range of 
temperatures investigated. 

The expansion ( |138| ) is perfectly adequate since a fit from its different terms plus a constant C to be 
adjusted gives « 208 MeV for the set ( pi] ) and C^'^ k. 182 MeV for the set ( |l45|) , in excellent 

agreement with the results obtained from the behavior of the pressure. 



This is in agreement with the results |128 based on chiral perturbation theory | 129 
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Figure shows that the order of the transition and the nature of the coupUng (strong or weak) between 
the quark condensates are very well visualized with the help of the energy density curves: there is a jump 
at Tc if the transition is of first order, while there is a change of slope quite visible if it is a true|^ second 
order transition. 




This is the case for the set (145), see figure ^(a): the transition is of second order for the up quarks 
while the strange quarks feel a first order transition. There is then a change of the slope of e{T) at T^c 
and an energy jump at Tgc, see figure p^(a). 

The behavior below the transition can be understood from cq. ( D39 ), with the same restrictions as in 
the case of the pressure, for the corresponding set of parameters. 

Figure [l6|(b) is the analogue of figure [l^(a) for = /^s = 250 MeV. Chemical potentials introduce a 
dependence in the simpler case of chiral limit. We have, with fi = fiu — fJ-s, 



20 2 



4 



B. 



(153) 



We mean a transition which is not a crossover, i.e. we are in the chiral limit. 
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This equation comes from ( |130| ) or, more directly, from ( |125| ) together with (111) and (121). Figure 16(b) 
shows also that having nis different from zero does not affect too much the second term of ( |l53| ), and 
introduces a constant supplementary term. All the remarks concerning the fits are also valid here. 



Entropy 



For a massless free quark gas, the entropy density behaves like T^. Since the strange current quark mass 
does not vanish, the high temperature expansion has correcting terms. Eqs. ( |128| ), (134) and (13S) give 
the first three terms of Ts (temperature times the entropy density) in powers ofT^, leading to eq. (pll| ). 
The complete expansion is gi ven b y ( p43| ). 

Figure |l^(a) shows that eq. ( |l4l|) can be approached by a linear expression in T'^ in the given range of 
temperatures: 



and 



13.868 - 0.136 10" 



s« 13.839 T^- 0.912 10~^ 



(154) 
(155) 



for the set of parameters (144) and (|145|), respectively. 



Although the numerical results seem to fit exactly t he relation ( 123), the coefhcients of (154) and 
(|155|) do not correspond to the combination of eqs. (147,148,151,152). 
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FIG. 17. Entropy density as a function of for the set of parameters (144,145) for fiu = Hd = f's = MeV 
(a); for tin ^ fJ-d = fJ.s ^ 250 MeV (b). 
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This remark confirms once more what we claimed in the pressure and energy density case: it is dangerous 
to extract information from fits if we do not take care of possible corrections. 

In our case, fits in T* (for pressure and energy) and in (for ent ropy ) do n ot hel p to get the bag 
constant or the number of excited d.o.f. The more complete forms (134), (138) and (141) have to be 
considered. 



The behavior of the entropy below the chiral transition can be understood using ( |D45| ) with the same 
warning as in the pressure case, for the corresponding set of parameters. 

Figu re |l7| (b) shows the entropy behavior for nonvanishing chemical potentials. Using eqs. (128), 



and (153), and the density (122) in the limit where all masses are vanishing 



(156) 



we get 



Ts 



7 

15 



(157) 



which has to be modified in order to ta ke in to account the finite value of m^. It is worth noticing that 
rus only slightly affects the term of (157). 



As a conclusion to these results, we can mention that, above the chiral transition, all the d.o.f. of the 



model are excited. However the vanishing mass limit of QCD is not reached for two reasons |127|: 



• the gluon d.o.f. are not included at high temperature, 
potential V^{T,fi) should be taken into account; 



A temperature and density dependent 



the strange quark mass is not negligible. We need going to very high temperatures in order that 
the lowest term in subsists in the expansion mg/T. 



We should also notice that the 1 /Nc corrections in references [ p0| , piy85| -|85|| show that the low temperature 
thermodynamics is driven by pion motion, pions being much lighter than the constituent quark mass. 
Since the quark loop contribution takes into account thermal excitations of quarks with mass M > 300 
MeV (whose probability is reduced by the Boltzmann factors exp{—(3M)), the low temperature thermal 
excitations are completely dominated by the almost massless pions. To obtain the effects of pions in our 
results, we should integrate over the meson fields in the path integral formalism, which is however beyond 
the scope of this paper. 

Even with the above mentioned limitations, the scaled NJL models have important new features: some 
gluonic effects are included through the gluon condensate x which couples the up and strange quark 
condensates. Thanks to this coupling, our model allows simultaneous transitions for the up and strange 
sectors (strong coupling), even though they tend to remain uncoupled for high chemical potentials. The 
coupling then allows first order transitions as a function of temperature while, within a pure NJL, they 
are always of the second ordeif^ (e.g. p8,131|). 



4- Comparison with lattice QCD 



To make a comparison with lattice QCPf^ it can be advantageous to normalize P, e, Ts and the 
interaction measure (e — 3P) to T**. The interaction measure gives the non-perturbative contribution 
to the thermodynamics: it vanishes in the Stefan-Bolzmann limit. It is also interesting to plot 3P and 
s in the same picture to see how sharp is the increase of the corresponding thcrmodynamical function. 
The origin of the coefficient 3 in front of P compared to e c ome s from the coin cidence of their respective 
asymptotic behavior (see the comparison between eq. (134) and eq. (138)). In the same spirit, one 
can normalize the entropy density by a factor 3/4 (see eq. (141)). In this way, 3P, e and 3sT/4 have 



The pure NJL model does, however, allow first order transitions w.r.t. density, e.R. p' 



Because lattice QCD has only turned recently towards finite density, see e.g. |132|, we restrict ourselves to 
/iu = Ms = 0. 
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the same asymptotic value TNc'k^ /2Q^ which is a direct consequence of the number of d.o.f. which enters 
the model. Note that the quantities we examine are relative to the quarks. In our simplified model, the 
glueball only enters through the bag constant. 

In figure ^ we show both the pressure an d the energy density of the A-scaling NJL model versus T/Tc 
(Tc = 150 MeV) for the set of parameters ( plsj ) {M^ = 300 MeV, xo = 80 MeV), with to„ = 0. Here, we 
have taken the critical temperature corresponding to the chiral symmetry restoration connected to the 

up quarks. 

Several interesting points have to be mentioned. One expects from lattice studies, e.g. [^,107 111 |, 
that the thermodynamical quantities are almost vanishing below Tc, then increasing. This increase is 
very sharp for e and Ts, while the pressure approaches the Stefan-Boltzmann limit very slowly. Lattice 
calculations show also that e/T^ has a peak^] just above T^, then approaching its asymptotic value from 
above. Finally, they also show that e — 3P ^ above Tc. Our results, summarized in figure ^ show that 
the model is in qualitative agreement with lattice results. The quantitative difference can be understood 
in the following way: lattice calculations show a rapid variation of the entropy density in a narrow region 
of T 10 MeV), which is traced back to the liberation of quarks and gluons. It seems then quite trivial 
to relate this fast increase to the confinement-deconfinement properties, which are not included in our 
model. This is clearly seen in the entropy density calculated with our model where the entropy is already 
increasing (although not as fast as near Tc) for T as low as 0. 2Te. Once this entropy curve is understood, 
the general behavior of P and e can also be deduced, see e.g. |133]. It is explicitly shown in that reference 
that, starting with a sharp entropy density, the energy density has a peak, and that the pressure increase 
above is low. In fact, would the entropy be approximated by a step, we should have the exact result 



PjT) 
Psb{T) 



(158) 



which gives P/Psb = 50% (90%) for T/T^ = 1.2 (1.8), independently of the details of the model. Since 
we are far from a step for the entropy, P/T"^ has an even weaker T dependence. This is shown in the 
general model for the entropy |133| and is confirmed by our particular m odel. On the same ground, it 



can also be shown that the interaction measure (e — 3P)/T , given in eq. (|131| ), has a peak above 
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FIG. 18. Pressure, energy density, entropy density and interaction measure for the set of parameters 
M° = 300 MeV, xo = 80 MeV. 



We have seen that the general behavior can be understood from the analysis of [133| which, together 
with the lack of confinement of our model, explains the quantitative disagreement between lattice gauge 
calculations and scaled NJL ones. However, our figure shows nice features not discussed extensively in 
the literature. If we concentrate on the energy density, it is clear that the peak has its slope broken in two 



This is not the case for a pure gauge theory. 
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places. These broken slope points coincide with the temperature where the chiral symmetry is restored. 
Since the current up quark mass is zero, the transition corresponding to the up quarks has no tail (see 
figure |l^(a)), leading to the first slope discontinuity while, because the transition of the strange quarks 
is of first order, there is in fact a jump in the energy density. Since this jump is small, it looks like a 
discontinuous slope. To get a nice peak, one then has to consider only crossovers (second order transition 
with nonvanishing current quark masses). Note also that a gap in energy only transforms into a change 
of slope for the pressure, while a change of slope in the energy plot is almost invisible in the pressure. It 
is evident that the energy density is the adequate quantity to be investigated in order to have insights 
on the order of the transition, and for extracting the critical temperature^. 

Figure |l^ illustrates that the broken peak of figure |l^ is due the combined effect of a second order phase 
transition for the up quarks (in the chiral limit = 0) and a weak first order transition for the strange 
quarks. We have taken the set of parameters (M° = M° = 400 MeV, xo = 350 MeV). In that case, 
there is only one critical temperatur e, an d the transition is of second order (see figure |^) , with the critical 
temperature given precisely by eq. (10£). 

All we have said for the energy density remains valid for the interaction measure (two peaks in figure |l8| 
which, with niu ^ and a second order transition for the strange quarks, would lead to a single peak). 
This gives, in the limit of three degenerate fiavors, the interaction measure of figure |l9|. Note that, in the 
chira l limit, the interaction measure just gives 4i?/T^ (for T > Tc), B being the bag constant, see eq. 



12 
11 
10 
9 
8 
7 
6 
5 
4 
3 
2 
1 




: A-scsJiiLg 










(3/4)s:t3 i 




spit'' 




^^^^ (E-3P);T* : 



FIG. 19. Pressure, energy density, entropy density and interaction measure for the set of parameters 
Ml = = 400 MeV, xo = 350 MeV. 



C. Anomalous decay ttq 77 



Until now we did not pay too much attention to the fact that the model is not confining, allowing the 
unphysical decay meson — > g + g as soon as the meson is more massive than the sum of its constituent 
quarks. If the unphysical width is small compared to the meson mass, it is tempting to argue that 
the model is still meaningful. We just need extending it above threshold, performing calculations, and 
checkingthat the unphysic al wi dth is small. Although there exist calculations which can be considered 
as exactQ above threshold [134|, w e pre fer to work in the small width approximation. The calculations 
we give here are adapted from |88 135]. As an application we show, in the 2-flavor version of the A- 
scaling model, the variation as a function of temperature of the (physical) decay width ttq 77. This is 



These informations can of course be obtained from the reconstruction of the quark condensates. 
®^By exact, we mean calculations which take into account the full unphysical width, for any value of its 
magnitude. 
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interesting for two reasons: i) it is an anomafous proces^^ and ii) we shall show that we expect a strong 
enhancement of the decay width near the deconfining transition. 

The calculation of the two-photon decay width of the neutral pion TT^^y^ requires the knowledge of 
the pion mass and the pion coupling constant to the quarks gTrqq. The lack of confinement above 
the threshold m.„{T) > 2M„(T) implies imaginary parts for both the mass and the coupling constant. 
Because of the two different scales connected to QCD and QED (as and a, respectively), the decay 
width r7r^77 has to be calculated non-perturbatively for the strong part, as it has been done previously, 
while the electromagnetic part can be treated perturbatively. To make things simple, we only consider a 
second order phase transition. Because we want a low critical temperature, we shall keep the set (145): 
(M° = 350 MeV, xo = 80 MeV). Taking a larger xo does not change the conclusions, except that the 
transition occurs at a higher temperature (see figure 0). 



1. Definition of the model above the threshold n 



qq 



Although there is an exact study |134|, we work in the small width approximation, in the same spirit as 
in [ pmj . 

To remove any ambiguity about the sign of the widthp^ F^^q^, we make a continuation of the model 
over the whole real line {q^ s] — oo, -I-cxd[), before searching for poles in the Bethe-Salpeter equations (or, 



equivalently, in the definition of the pion mass as described in section [ D 4 ) . Other ways can be found in 
the literature. For example. 



the authors of [136| define the mass and width above threshold by looking at the weight position 
and the half-height of the imaginary part of the propagator, respectively. Below threshold, this 
imaginary part is reduced to a Dirac distribution whose position is identified with the pion mass. 
How ever this approach has been criticized in [137| in the vector sector. In the more recent reference 
[134 1, the authors show why there is a discrepancy and give a way to cure it; 



• the authors of |]38| simplify the problem, identifying the expression above threshold with its Cauchy 
principal value. Although very simple, this method completely ignores the process tt — > which, 
while unphysical, is nevertheless present. 



The method we have developed in |gS,135| is the following: let the inverse propagator be G~^{qQ) (we 
work in the static limit g = 0), in Euclidean space. The pion propagator^ is then defined in the range 
- 4:M^,+oo[. We extend it in the range ] - oo, -4M^] by defining 



hm G-\q^o - = MqI) + ^B{ql), 

£—►0 



(159) 



where A{q^) and B{ql) are real functions of the real variable q^. Below threshold, B{ql) vanishes identi- 
cally. With this definition, the pion propagator is defined for real Qq. The pion mass and its unphysical 
decay width are obtained from this expression, by searching for the zero of the propagator in the complex 
plane: 



G-^(-(m^ - iF/2)2) ^ A{~{m^ - i^f^f) + iB{-{m^ - iT/2)'^) = 0. 



(160) 



Since qg is now a complex variable, the quantites A an d B are also complex. We write A = Ab. + iAi 
and B = Br + iBj, so that looking for the zero of (16C) implies 



Ani-im^ - ^F/2)2) - 6,(-(™- - *r/2)2) = 0, 
Aii^im^ - iF/2)2) + Bni-im^ - iT/2f) = 0. 



(161) 
(162) 



We apply this method to the inverse pion propagator (eq. (CI)) and we write lim Z^{q^ — ie) = A{q^) 



iB{qpi). The analogue of (|160|) is then 



It is here a QED anomaly which does not have to be introduced by hand as the strong axial anomaly or the 
scale anomaly. Here, the anomaly shows up as a result of the introduction of the electromagnetic field. 
^^We do not have any it as in Minkowski space. 
^^For the kaon, the threshold is at % — -~{Mu + Msf' . 



53 



+ i{Ai{-{m^-iT/2f) + BR{-{m^^iT/2f)] ] +a'x'j^ = 



(163) 



In the limit T/2 << m^, we can neglect Bj and Aj in front of and B^, so that this equation can be 
rewritten in the form 



We just need solving^ 



(m. - ^^/2r {Ai^mD + ^B{-rr^l)) + a\l'^ = 0. 



(m. - rYl2fZ-{~ml - zO+) + a\l'^ = 0. 



(164) 



(165) 



This equation is valid in the chirally broken phase, where we made use of the gap equation (68) to obtain 
(CI). The inverse propagator which is valid in both phases is given by 



Below threshold, we then hava^ 



- ?>NcgM^,l3.Si - m^Z^i-m^) + hir seal = 0, 

while, above threshold, in the small width approximation r/2 << m^^, we have to look for 



8iV,5A/„,^,^ - (m, - ir/2)2z^(-m2 - ^0+) + 



0. 



(166) 



(167) 



(168) 



Eq. ( |168| ) is a system of two coupled equations for two unknowns and V. They are however decoupled 
because the half- width is neglected in Z". These equations are 



B(-ml) 



ml + ml 



A^{-ml) + B^[-ml) 
A{-ml) 



A^{~ml) + B^{-ml) 
4 



(^A'^{-ml) + B'^{-ml) 



(8iVc5M„,/3,p-a'x') 

(8iV,5M„,/3,M-a'xD' = 0. 



(169) 



(170) 



In figure we show the pion mass and twice the constituent quark mass versus temperature. An 
unphysical width appears as soon as the threshold = 2M„ is reached. 

The usual definition of the pion to quarks coupling consta nt g^ po starts from the polarization operator 
( J^'^ in the notation of 1 135 1 or 11 in the one of Klevansky |101]) 



g-rrqq 



(171) 



uj—mTr—ir/2 



In our formalism, the coupling constant can directly be expressed from the pion propagator through the 
identification 



1 



g-rrq 



(172) 



^^Our small width approximation is not exactly the same as in [L34]. They are however equivalent by definition 
of a sr nall width approximation. 
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In [135], we work in term of a polarization operator to facilitate the comparison with Klevansky |101|. 
keep here the same notation as in the previous sections. 



We 
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which is consistent with the identification of the on-shell coupHng constant through the action. 



> 

(3 



0.7 
0.6 
0.5 
0.4 
0.3 
0.2 
0.1 




M^ = 300MeV,Xp = 80MeV 



1 1 1 1 1 1 r— 

A-scaling 


— 1 — 1 — 1 — 1 — 1 — 1 — 1 — 1 — 1 — 1 — 1 — 1 — 1 — 1 — 1 — 1 — 1 — 

2M '■ 






\ m ±r /2 z 




\ h ^ 


" ■ , ■ , 1 , , 


' 







0.05 



0.1 0.15 
T(GeV) 



0.2 0.25 



FIG. 20. Twice the constituent quark mass (plain line), the pion mass (dashed line) and the unphysical width 
Ttt (dot-dashed line) versus temperature. 



Indeed, withg 

a Taylor expansion around = (below threshold) gives 



G'^\q^)^G-\~ml) + 

By definition of the pion mass, G^^(— m^) = 0, we then have 

d{q^Z:{q^)) 



[q^ + ml) + 



(173) 



(174) 



G-\q' 



{q'+mi) + ..., 



(175) 



I.e. 



dq^ 



g^qq "^{-ml) = 



(176) 



whose generalization above threshold is obvious: 

-2, , ^{q^Z:{q^-^Q+)) 

gTrqq (-(TOtt - «r/2) ) = 



dq^ 



This is nothing else than (small width approximation) 



g2=_(m^_ir/2)2 



5.,-,-'(on-shell) = ^ [Z^i-ml - iO+) + Z^iO)] + ^{m^ - iT/2fM-ml - iO+), 



(177) 



(178) 



'Because g = 0, we have f{q) —* /(go,0) — > /(go) foi' function /. To simplify the notation, we write 



2 2 

Qo ^ q ■ 
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where 



1 

■(fc2+M2)[(fc-g)2+Af2]2 



(179) 



Note that, with the approximation sr^gg(— m^) « oWqCO) = l/\/.^^(0) valid at low temperature, where 
pions are essentially massless, we have, with eq. (108) 



(180) 



which is the Goldberg er-Treiman relation. With A/^ = 300 McV and U{T = 0) = 93 MeV, eq. ( |l80| ) 
gives Q-Kq aj—rn^) « 3.23, in excellent agreement with an exact numerical calculation, as illustrated in 
figure ^ This figure gives the behavior of the real and imaginary parts of the on-shell coupling constant 
g-Kqq given by (178), together with its modulus]^. Near threshold there is a huge effect: the coupling 
constant goes to zero. This is explained by the fact that ImJa, and then Im gTTqq~'^, are infinite at this 
point. 
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FIG. 21. The pion to quarks coupling constant g-^qq versus temperature (plain line), its real part (dotted line) 
and its imaginary part (dashed line). 



2. Disintegration ttq — + 77 

The process ttq 77 can only be explained by considering it as an anomalous process. Because the 
gluons are not present in the model, the strong axial anomaly has to be included by hand. In that way, 
this quantum breaking of the axial symmetry can be treated at tree level. However, the electromagnetic 
field is included in the formalism. That means that we can have access to the full quantum effects 
without having to add a term by hand. Note that we have to go beyond the tree level approximation. 
The electromagnetic field can be studied as a perturbation {a w 1/137). From a formal point of view, the 
amplitude for the process ttq 77 can be obtained starting from the bosonized action (|6l|). The term 
Tr In is the fermion determinant and contains all the quantum effects of the fermions. The derivative ^ 

has to be replaced by the covariant one Ip = !^ -\- iQj^ , where Q is the charge quark matrix diag(e,j, e^) 
2 1 

— diag(— e, — e). An expansion of the logarithm has to be performed until it includes the considered 



To simplify the terminology, we define from now on the coupling constant as being its modulus. 
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process|_J. This is very tedious so that it is much more efficient to work with Feynman diagrams. Two 
of them have to be examined: the direct diagram (fig. ^(a)) and the crossed one (fig. |2|(b)). From the 
action (|6lh, it is clear that the pion vertex isj^ g-K^qilb^i (i = 1, 3). Vertices connected to photons are 
given by[j e*^{ki)Q-ff, {i = 1, 2). 
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FIG. 22. Feynman diagrams for the decay width ttq 77; direct diagram (a); crossed diagram (b). 



The flavor trace over the amplitude associated to figures p3(a) and P3(b) is 



tr 



/ TTo %/27r+ 
\ \/27r^ -TTo 



en 




e« 
ed 



ell 



(181) 



which shows explicitly that only the neutral pion has a two-photon width (as it should from electric 
charge conservation) . 

The color trace gives a factor (global gauge invariance of the model) , so that only the Dirac trace and 
the loop contribution have to be determined. The transition amplitude for the process ttq 77 associated 
to the direct diagram ^(a) gives then (in the Euclidean space and in the isospin limit m„ = to^) |101 135]: 



-trD 



{75(^2 +^-^«)7mW3 + ^'^«) 



X 7, 



,(^i+M„)} 



1 



{ql + Ml){ql + Ml){ql + Ml) 



(182) 



which has the general form, as shown by tr/^, 

%,^^^ cx e^'"'^el{ki)k^,el(k2)k2p, (183) 

which could have been guessed on the basis of Lorcntz, parity and gauge invariance ipsl. The summa- 
tion over photon polarizations in (|182|) leads to 



•"■0^77 I 



where 



Finally, 



JM = 



5127r2a2_| \g^^^Mlj^{-m^2 - iO+)\ 



1 



AT) 



(27r)4 [(93 - k,y + M2][(g3 + k^f + Ml){ql + Ml) ' 

d^ki <Pk2 



(184) 



(185) 



1 1 



2 2m^ (27r)2 J 2\kA 2\k 



l'^o-*77 



I 5{q-„ - fci - k2) 



1 



327rm7r 



1^0^77! ' 



(186) 



^^Electromagnetic corrections to a given order are expected to be small because a is small. The strong corrections, 
which should be included because the strong coupling constant is not small, are eliminated on the basis of the 
1/A'^c ordering. 

^^The pion to quarks coupling constant enters the vertex because we have normalized the pion wave function 
in order to get a usual kinetic term for the pion propagator (see eq. (175)). This corresponds to a canonical 
normalization. 

^^e^ is the photon polarization. 
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where the factor 1 /2 takes into account the presence of t wo id entical particles in the final state 
The temperature behavior of r7ro^77 given in figure 
unphysical process ttq — > qq. The plain line gives TTro^-y-y 
imaginary part of J4 and g-nqq" ■ The pion mass is still given by ( 170 ) 



The dashed line takes into account the 
where we artificially removed by hand the 



Note that a pure NJL would give the same figure except that the peak would be shifted to a larger 
temperature (around f90 MeV instead of 140 MeV). This is due to the fact that scaled NJL models 
incorporate the coupling of quarks to the gluon condensate. The zero temperature result is however not 
changed from a pure NJL because xo enters into the game only through the combination Axo, which is 
fixed by the pion weak decay constant fj^- 

The numerical value of F^^^-y^ is very sensitive to the actual pion mass, being proportional to m^. A 
small pion mass error on J4(— m^) and g^^qq can also have some importance. A great numerical precision 
of is then required, and also for quantities which depend on itP|. 
At finite temperature, we get the table || where the experimental value is from |104]. 



TABLE II. Decay width r^o_ 





Exp. 


NJLa=0.9 GcV 


NJLa=oo 


Non-local 


A-scaling 


r^o^^^ [eV] 


7.7± 0.5 


5.50 


8.75 


7.22 


8.75 



We also compare our results with an instantaneous 3-dimensional non-local NJL model [ B4 , 100 139| . 
An important feature has to be noticed: each value given in table || is relative to a cut-off put to infinity 
in the convergent integrals, except for the third column. As one can see, this prescription is indispensable 
to reproduce the experimental result (see also |140[ ). One can wonder about the prescription "putting a 
cut-off only for diverging quantities" . Does the consistency of the model not require to keep the cut-off 
everywhere? However, it is shown in the work [141| that a sharp cut-off is responsible for this fact. It 
then seems reasonable to accept the prescription. One can also argue [|l^,^ that the diagrams of figure 
are the only ones to contribute to the anomaly if the cut-off is infinite, as in case of a renormalizable 
theory. In a non- renormalizable one, the ratio M^/A is finite and other diagrams have to be included 
(see appendix A of 133]). 
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FIG. 23. Decay width r,ro^77 with the unphysical process tto —> qq included (dashed line) or not (plain line). 



Because of the behavior of g,rqq, the width goes to zero at the threshold, then increases sharply. This is clearly 
an artifact of crossing the threshold, which would be removed provided we would smooth gT^qq around it. Because 
of this artifact, we do not draw this part of the width. 
'^^This remark implies that we should take into account the mass splitting between ttq and tt* in the determination 



of 

r,ro^77- Table O has been obtained with — 139 MeV, as is usually done in the literature, e.g. jlOll. One 



expects a better agreement with experiment if the value m-,, 
the case. 



135 MeV is chosen. We have verified this is actually 
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If the cut-ofF is kept for each integral in the NJL model, these diagrams should then have to be taken 
into account. One can imagine that these diagrams would compensate for the discrepancy between the 
third and fourth column. However, it is simpler to only consider the diagrams of figure ( p^ and to take 
an infinite cut-off wherever possible. The functions J3 and J4 are evaluated in this way. 
Figure (|23|) shows a sharp peak. It is present both in the scaled m odel and in a pure NJL model 



|101, 135[142| . At high temperature, it goes to zero. We have shown in |135] that the peak is sharper in 
the scaled model than in pure and non- local NJL models|3- It is located at the pion threshold temperature 
which is almost coincident with the critical temperature (it would be exactly the critical temperature in 
the chiral limit). The peak is mainly due to the behavior of J4. 
As a conclusion to this section, we can write that: 

• taking the cut-off to infinity whenever possible (table ^ gives similar results in pure NJL and A- 
scaling models for r7ro^77 zero temperature. There is a 14 % di screpanc y with the experimental 
result 1 104 1 . The non-local NJL model of the Rostock group [184, 1001 , 1391 is in better agreement 
(6 %). The footnote (75) suggests, however, that the discrepancy between theory and experiment 
i s connected to the pion mass we take in our calculations {m^ — m^± — 139 MeV). Indeed, with 
( [l80| , p^ , p^ and M-ml) at zero temperature and density 



J4{-ml)f3^ 

we straightforwardly obtain 

r 



1 



■■oo,fi,—0 



2AL 



647r'^/2 Y TO 



/2M„ . TO^ 
arcsm- 



2M„ 



(187) 



(188) 



With the parameters from this section, eq. 188 gives TT^o^-yy = 8.65 eV, v ery close to the exact 
theoretical value (8.75 eV). This shows that the Goldberger-Treiman relation (180) is a very accurate 



approximation. If we had chosen to^ = 135 MeV, eq. ( |188| ) would have given r^,-,^ 
which is now very close to the experimental value given in table ^ 



►77 



eq. (188) can still be more simplified. In the limit TO7r/(2Af„) 

rv m'^ 



0, we have 



•"•o— »77 



64^3/2 ■ 



= 7.91 eV, 



(189) 



identical to the result quoted in [E2[. For m-^ — 135 MeV, this gives 7.64 eV; 



• in the A-scaling model (and in the non-local NJL model) |135|, the decay widthTVj,^^^ has a sharp 
peak near the critical temperature. This is similar to a pure NJL model [101 142| , except that the 
peak of the latter is rejected to a higher temperature; 

• above threshold, the unphysical process ttq — > qq is possible, since the models are non-confining. 
The pion acquires an unphysical width, gT^qq and J4 become complex quantities. Considering this 
process consistently in the model leads to a non-negligible influence on r^n^-y^ versus temperature, 
as shown by the comparison between the plain and dashed curves in figure E3^. Note that, taking 
this unphysical process into account, we still have a peak, in contradiction with [101[; 



the constituent up quark mass which enters the calculation of F. 
through the pion mass) is extracted from the gap equation (p 
the electromagnetic field, this equation has to be modified [[143 144 



j-y (either explicitly or implicitly 
In the case of a coupling with 
It is however shown in these 



references that it is justified to limit oneself to the usual gap equation for a result valid up to order 
0(a); 

• it should also be mentioned that our results are at complete variance with those of Pisarski et al. 
[145|,146{, who obtain a continuous decrease of the decay width with the temperature. 



^^We have shown in the B-scaling NJL model |135| that the two- photon decay width of the neutral pion goes to 
zero at the transition. This is connected to the fact that, in this model, the pion mass goes continuously to zero 
below the chiral transition, and vanishes at the transition. The width r^o^.y^ then vanishes at T = Tc. However, 
this is true only in the two-flavor version of the B-scaling model, see |63|. 
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III. MIXING BETWEEN SCALAR ISOSCALAR MESONS AND THE GLUEBALL 



In the previous section, we focused on the conventional pseudoscalar channel, plus the effect of a scalar 
glueball. Because of its scalar nature, it can be interesting to consider this channel. The scalar glueball 
is an isoscalar. It can then couple to similar d.o.f. In our three flavor models, this implies a coupling 
between three particles: the glueball and two qq-\ike mesons. In our formalism or in GCM type models, 
these mesons are associated to the canonical vertex 1 in Dirac space. This is probably a very crude 
approximation for the scalar channel. However, as a first approximation, it can give insights into the way 
of treating these mesons. 

The interest in scalar mesons mainly lies in the search for the scalar glueball, a particle which should 
show up in QCD because of the color charge of the gluons, and which is seen on the lattice (see below). 
In our three- flavor model, the scalars form a set of ten particles (singlet, octet plus glueball). The 
identification of the scalar nonet is anything but clear, both experimentally and theoretically. Does it 
have to include the ao(1450), /o(1300), /o(1590), (1430) as suggested by Montanet ||l4^, leaving the 
ao(980) and /o(980) as KK molecules and the /o(1500) as a solid candidate for the glueball? Or should 
the nonet include the /o(980), as proposed by Palano [ ]148[ (in that work, this /o(980) is also said to 
be mainly a ss state)? Tornqvist [ |f 49| presents the mesons /o(980) and /o(1300) as two manifestations 
of the same ss state. Lindenbaum and Longacre |150| claim that the 0"*"^ data can be reproduced with 
the four mesons /o(980), /o(1300), /o(1400) and /o(I7IO). In the latter case, it is also assumed that the 
61(1710) = /o(1710) has a spin 0. 
Theoretically, the situation is also very obscure. 



For example, Miinz, Klempt et al. |151 , 152| | calculate 
masses and decay properties of the scalar nonet in a relativistic model with an instanton-induced interac- 
tion and a linear confinement. They give support to identify the nonet with ao(1450), iCg (1430), /o(980) 
and /o(1500). (There is no glueball in their mode l.) L attice QCD, in the UKQCD Collaboration, pre- 
dictsa scalar glueball of mass around 1500 MeV |153|. This favors its identification with the /o(1500) 
[ ^ , |154[ . However, the valence (quenched) approximation of Sexton et al. gives a glueball with a mass 
of approximately 1710 MeV, suggesting its identification with the /j(1710), provided its total angular 
momentum is J = 0. The coupling with qq states would then increase the glueball mass of about 60 MeV. 
The fact that th e glueba ll disperses over three resonances is also shown in the K-matrix fit form alism of 

Anisovich et al. [ |155 ,156| which, however, favors a glueball at 1500 MeV, as also shown in |157|.^ 

This scalar resonance /o(1500) has been discovered at LEAR by the Cry stal Barrel Collaboration | 15^ , 159|1 
and also seen|^ by this group in the decay /o(1500) KK MeV [ 16_l[ 
between the glueball and the neighboring qq states, as shown in [ j4|. 



Note that this requires a mixing 
As already mentioned. Sexton et 



al. |93| have also noted a mixing bet ween the glueball and the qq states. In their work, and in subsequent 
researches from one of the authors |162 ,|l63 |, the glueball is the /j(1710). Its mixing with /o(1370) and 
/o(1500) in turn requires that J = fgj. 

Although we have varied the parameters of the scaled NJL model in order to examine both t he / o(1500) 
and the /o(1710) as possible glueball states in the study [164|, we shall mainly focus on ]165| ] in the 
following, and restrict ourselves to the /o(1500) as the scal ar gl ueball. 

We concentrate on the minimal list of light scalars given in |104]. It includes the /o(400 — 1200), /o(980), 
/o(1370), /o(1500), ao(980), (W1450) and i^o(1430). We also consider the heavier /o(1710) as a ss 
candidate. We give in section ID 4 and appendix ^ the masses of the various scalars within our model. 
It is clear that, without coupling with the glueball, the ag and /o,uti are degenerate and, in the chiral 
limit, of mass m = 2Af". The /o,ss has then a mass m = 2M^ and tuk* = + M°. To make contact 
with our model, the /o(400 — 1200) has to be rejected without explanation: there is no room left for it 
in our formalism. One has also to choose between the sets (/o(980),ao(980)) and (/o(1370),ao(1450)). 
(Because of the meson masses in the chiral limit as given above, it is not possible to have, for example, 
/o(980) with ao(1450).) To compare our work with Q, we choose to keep the set (/o(1370),ao(1450)), 
noticing th at th e /o(980) and the ao(980) can be interpreted as KK molecule^^ 



In section HI A 



we describe how the scalars are treated in the scaled NJL model. In section HI B| , 
we study the mixing angles and meson to quarks coupling cons tants. Finally, results concerning the 
two-photon decay width of the scalar isoscalars is given in section HI C . 



In the Mark III data on the radiative decay J/tp 4n + a sharp peak is se en in the 47r spectrum, having 
the quantum numbers 0"*"^ and a mass and width compatible with /o(1500) MeV |160|. 
For the alternative choice, see [164 . 
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A. Hybrids: from unphysical pure nn, ss and glueball fields to physical particles 



In appendix ^ we give the general structure of the matrix propagator of both pseudoscalar P and scalaif^ 
r particles. In the following, we are interested in the latter. The cases of the qq and the Kq are trivial 
since these particles are unmixed. However, there are three coupled states in the isoscalar sector: two 
members of the nonet (the singlet and the eighth part of the octet) and the gluebalQ. Going from mixed 
to physical states can be treated just like an eigenvalue (e.v.)/eigenvector (E.V.) problem and is discussed 
in some detail in appendix This is generally a non-trivial problem because of the relativistic nature 
of the theory: a relativistic theory does not allow one to describe observable states as superpositions of 
other states with different masses. Mathematically, this is traced back to the energy dependence of the 
r matrix ele ments. O btaining the physical masses is however simple: one has just to look for the zeroes 
of the r e.v. |164,165|. The only difficulty consists in attaching the right e.v. to the right particle. This 
can be done by varying (slowly enough to follow the evolution of the e.v.) the coupling to the glueball 
from zero to its actual value. We can also use the analytical solutions for the e.v. Since we have a 3 x 3 
problem, we have just to look for the solutions of a cubic equation^ Once we have identified to which 
e.v. corresponds each state (this is done by putting to zero the couphng to the glueball), we can switch 
on this coupling without having to care about the e.v. ordering^. We proceed as follows: for a vanishing 
coup ling to the glueball (matrix elemen ts (Z!12 ,G13) put to zero), the glueball mass can be taken from 
eq. (|Cll|) . We can then identify from (19C-192) the e.v. to attach to it. The remaining 2x2 system 
(eqs. (|C8|- C1C )) can be exactly diagonahzed, giving for one of the /o the same mass as the oo, eq. (C6). 
The other /o ha s the sam e mass expression, except for instead of Af". We can then also attach 



one of the e.v. (19C - 19p ) to it. Since we know which particle /o(1370), /o(1500), /o(1710) is attached 



to each e.v. in ( 19Cl| - 192| ), we are then able to study quite easily the mass evolution with respect to the 
parameters (M^,xo). Because the meson to quarks coupling constants can be written as an eigenvector 
(E.V.) problem (see next section), the identification of the e.v. to their corresponding state is also of 
prime importance to determine these couplings. 

In the preceding sections, we used a glueball mass of 1300 MeV in order to find the parameter in a 
sigma model limit of the scaled NJL model. This was done to obtain without having to think about 
the scalar sector. It is clear from this section that this is only an approximation. In fact, as this study 
shows, one has to impose the value of the glueball mass. Once its e.v. has been identified, one can look 
at the val uef^ of the glueball mass without coupling. Since it has to coincide with the mass extracted 
from eq. (Cll), we can then deduce b^. 



^^In a 3-fla vor N JL model with a 't Hoo ft determinant, scalars have already been investigated by Hatsuda and 



Kunihiro [p8Ul6f| and by Bajc et al. 167 1 



We are working in the isospin limit, which explains why all the components of ao and Kq are unmixed and 
why the third component of the octet does not contribute to the isoscalars. 
*^The cubic equation x'^ + ax^ + 6a; + c = has the three real solutions 

a;i = -2VQcos(^)-^, (190) 

X2 = -2yQcos(^±^)-|, (191) 

x^ = -2^cos{^-^)-^, (192) 

(193) 

with e = arccos(i?/ y/q3) and Q ^ {a - 36) /9, R = (2a^ - 9a6 + 27c) /54, provided that < and R, Q are 
real. One has then to find a one-to-one correspondence between each of these solutions and the physical states. 
This is performed by comparing these solutions to the known mass of the mesons when the coupling with the 
glueball is removed. 

*^In the purely numerical case, this is not possible because the e.v. are not ordered once for all. 

*^With the approximate method of the previous sections, we get 6 = 4.69 for (M^ = 450 MeV, xo ~ 350 MeV), 
see table ^ In the correct procedure developed here, we have b = 5.02. Should we have taken a glueball mass of 
1500 MeV in the preceding sections, as done here, we would have got 6 = 5.16. This shows the consistency of all 
the performed approximations. Note that, although playing almost no role on the thermodynamics and on the 
variation of the condensates as a function of temperature and density, this difference in 6^ is important for the 
scalar sector and explains why we have to take the exact expression in this section. 
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Another subtlety shows up in the scalar sector: scala rs are always unbound, lying above the unphysical 



threshold fo ^ Q + Q- The technique of section P C 1| has to be applied. Note that this makes sense only 



if this unphysical width is small. In our studies, this is always the case except for the /o(1710). However, 
it can be argued (see below) that meaningful quantities, such as two-photon decay widths, can still be 
extracted. 

In the following, we want to calculate the two-photon decay widths of the scalar isoscalars. Close, Farrar 
and Li have used general ideas to estimate the relative strength of these decay widths. They consider 
two schemes for the mixing between /o(1370), /o(1500), /o(1710). Firstly, they consider the case where 
the /o(1500) has a large glue content whi le /o(1710) is mainly a ss excitation. This scheme allows to 
understand the /o(1500) KK data |l6l[| . T he secon d scheme consists in considering that the glueball 
lies above the scalar nonet, as sugg estedin [p p2p^ : /o(1710) can then be the glueball while /o(1500) 



is mainly a ss excitation. Within these two schemes. Close et al. |^ have estimated the relative strength 
of the 27 widths for the three /o states: 

/o(1370) : /o(1500) : /o(1710) w 12 : 1 : 3 (scheme 1), (194) 

/o(1370) : /o(1500) : /o(1710) w 13 : 0.2 : 3 (scheme 2). (195) 

For any mixing, the decay width of /o(1500) is always the smallest. This suggests that the experimental 
estimation of the 27 widths might be a good test of the general idea of qq and glue mixing. 
We sha ll restric t ourselves to scheme 1. Note that, although derived using very general assumptions, the 
ratios ([l9^ , p^ ) are obtained assuming M° = M° (51/(3) symmetry) and that loop phase-space effects 



can be ignored (see eq. (6.2) of |9j]). We show in the following that these effects are in fact crucial, at 
least for the crude model we are using. 

Considering that the three scalar isoscalar mesons we want to investigate are /o(1370), /o(1500), /o(1710) 
leaves little choice for the value of M^: without coupling with the glueball, and in the chiral limit, we 
would get wjp(i37o) = 2M^. The state oq has no coupling with other particles and its mass is exactly 
2M ° in the chiral limit. This leads us to take ao(1450) as a member of the scalar nonet, rejecting ao(980) 
(and (/o(980)) as KK molecules. We then take Af° = 725 MeV. The mass of the isoscalars is dependent 
on the choice of xo- However, we immediately get TOao(i450) = 1450 MeV and mx'(u-io) = 1600 MeV, 
the latter being quite poorly reproduced. We shall not focus on this problem here. We show in figure ^ 
the Xo behavior of the mass of the three /o states. 
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FIG. 24. xo behavior of the mass of the three states /o(1370), /o(1500) and /o(1710). 

For large values of xo j the states take the mass they would have without coupling (m (1370) = 2Af " = 1450 
MeV , m/j,(i7io) = 2M° = 1728 MeV in the chiral hmit). The mass of /o(1710) is quite stable, while an 
acceptable value for the mass of /o(1370) limits the value of xo to a domain between 300 and 450 MeV. 
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B. Mixing angles and coupling constants as an eigenvector problem 

We explain in appendix ^ how to go from mixed states to pliysical states. Tfie derivation sliows how to 
obtain the masses (they are the zeros of the e.v.) and how to define the physical states. This enables 
to get the meson to quarks coupling constants. Similarly to what has been performed in the ttq —^ 77 
case, the decay of scalar mesons into two photons is represented as a triangle diagram. Since the glueball 
field is not pure (it mixes with quarks states), two types of diagrams have to be investigated, as shown 
in figure ^ 




FIG. 25. Two-photon decay of the scalars. 

In the limit of a large xo: the glueball decouples from the nonet. This is then the vacuum gluon condensate 
which fixes the mixing^: the smaller xo^ the larger the mixing. This is confirmed in the mass plot, figure 
|2^ . This plot shows that, in the range of xo investigated to reproduce the meson masses, the /o(1710) 
remains mainly a ss excitation, while the strength of the coupling to the u and s quark has about the 
same order of magnitude for /o(1370) and /o(1500). This has of course to be checked, looking at the 
meson to quarks coupling constants. This is the subject of this section. 

In our model, the three scalar states have a glue content. As mentioned, the /o(1500) is the one which 
yields a pure glueball if no mixing. This state has the larg est glue content even if that of /o(1370) can 
become important. We then work with scheme 1, eq. ( |194[ ). 

A mixing between the three glueball and scalar fields modifies the transition amplitude. Without coupling, 
the /o(1370), for example, can decay into two photons only through the production of a uu pair. As 
soon as a mixing with the other states is allowed, the production of a ss pair has to be considered. The 
coupling constants can be obtained from appendix R. We recall here the main r esultpl. 



k,l 



with {g^ = 0- V* is the E.V. matrix which allows going from mixed states to physical states and 



g is the mixed scalar to quarks coupling constant matrix. Eq. (E38) is the basis to obtain the physical 
couphng constants to quarks. Indeed, the functional (^) contains the term (g(/3"r")g, where denotes 
the mixed states. For the scalar isoscalars, this gives qaaXa/'iq. Using the relation going from the mixed 
states to the physical ones, 



it is clear that the meson to quarks coupling constants can be identified fromf^ 



q{V4-mj)g4-mf}j:)^^q, (196) 



**This is also the o nly p arameter at our disposal since we have fixed to reproduce the ao(1450). 
Note that in ref. |165| we use a prediagonalization for the system fo(1 370), /o(1710) which leads directly to 
the physical states if no coupling with the glueball. In that case, eq. ( E37 ) can be used 



86 



The conventions of appendix |e| are used, where we denote by A^q^ = — m^) a matrix A of which each element 
is evaluated at the same q^, in opposition to ^(—771^) which indicates a matrix whose first column is evaluated at 
—mi, the second at — m2, the nth at — m^. 
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which imphes that the meson to quarks couphng constants are (no summation over i) 



(197) 
(198) 



where g.{-m% = g.{-mj)u and z = 1 = /o(1370), i^2 = /o(1710), i = 3 = /o(1500). 
This expression clearly shows that, as soon as the E.V. matrix and the meson masses are obtained, 
the couphng constants can be algebraicaUy constructed. They are plotted, as a function of xo, in figures 
H and 1^, for g*"" and g*^"*, respectively. 
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FIG. 26. xo behavior of the coupling constant to up quarks, g"^ , of the three states /o(1370), /o(1500) and 
/o(1710). 
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FIG. 27. xo behavior of the coupling constant to the strange quarks, g"" , of the three states /o(1370), /o(1500) 
and /o(1710). 



We see from picture that the coupling constant of /o(1710) to the up quarks is small, vanishing quickly 
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as xo is increased. This shows that this meson can be considered as an almost pure ss excitation^. Note 
also that, because this coupling constant is small, the processes involving up quark diagrams for this meson 
will be suppressed. This is interesting since the threshold coming from 2M° = 1450 MeV < 1710 MeV 
implies that the unphysical decay width into free quarks could be important, so that i) the small width 
approximation would not be correct; ii) unphysical effects would radically change physical quantities. 



Because 



'•/o(1710) 



) is small, it is however clear that these unphysical processes play no role. 



This picture also shows that, although figure ^indicates that the masses attain almost their asymptotic 
value at large xo (no coupling of the quark states with the glueball state), the quark content of the 
glueball is not negligible. This is seen from the nonvanishing of 9 {—'m''fg (i^oo) 



) and in picture 27, The 



latter also shows that the /o(1370) has still a ss component for xo as large as 1250 MeV. 
The conclusion we can draw from these plots is that the glueball to quarks coupling constant is a slow 
decreasing function of the vacuum gluon condensate xo, showing that the glueball has a non- negligible 
quark content. This is particularly true for values of xo in the regime giving "good" masses for the scalar 
mesons {xo G [300, 450]Mey). Moreover, for this set of parameters, the /o(1370) has also an important 
strange quark content. 

Having the coupling constants, it is now easy to determine the two-photon decay width. 



C. Two-photon decay of scalar isoscalars 



With the meson masses and the coupling constants to quarks, we are in a position to calculate the two- 
photon decay of scalar isoscalars. Without coupling with the glueball, we would just have to evaluate 
diagrams similar to figure 22. They are given in figure Eg, where we omit the crossed diagram. 
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FIG. 28. Direct Feynman diagrams for 2'y decay o//o(1370) (left) and o//o(1710) (right). 

Since, in the scaled NJL model, the mesons /o(1370) and /o(1710) are not anymore pure uu or ss 
excitations, these diagrams have to be supplemented by similar ones with a permutation of the quark 
flavor. We then arrive to the following amplitudes: 



^/^(laTO) = N,{el + 6^)5^1370) <(fci)4(fc2)C/ + ^ce^<?;„7i37o)<(^i)4(fc2)5, 



(199) 
(200) 



for the /o(1370) and /o(1710) states while, having some qq excitations, the glueball can also decay to two 
photons: 



>o(i5O0) = N,{< + e^)g;.7i500)<(fci)<(fc2)f/ + N,eig}ll^^^^el{k^)el{k2)S. 



(201) 



In eqs. (199-201), we have written 



(f/,5) 



MlMis + < J7.(^i + Ml,)} 



ill + {Ml,Y){ql + {Ml,Y){ql + (MO J2) 



(202) 



'^This confirms what was only intuitive in figure 
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which is similar to the trace quantity appearing in eq. ( |182| ) except for the 75 which is due to the 
pseudoscalar nature of the pion and which leads t o the e^""'^ structure of eq. (183). In the scalar case, 
the amplitude still leads to the function J4 of eq. (185) that we write J„ and Jg according to the flavor 
which enters the equations. We finally get, after a summation over the photon polarizations (and for an 
incoming meson at rest) 



2 2m, (27r)' 



2 ki 2 h 



327r rrii 



(203) 



with 



2 _ 12800^2„,2 



9 



-TT a 



{ml - AMf)Mlgf^U^mi) + ^(m^ - 4Mf )Af0gr"V.(-m2) 

5 



(204) 



where i = /o(1370), /o(1710), /o(1500). 

With these equations, we are now able to plot the two-photon decay width of the scalar mesons, see figure 
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FIG. 29. xo behavior of the 27 decay width of the three states /o(1370), /o(1500) and /o(1710). 



The three widths decrease with increasing xoi as it should since, in the limit xo ^ cOj one of the 
corresponding coupling constants vanishes. fQ(niQ) is always the smallest: the first term of the r.h.s. 
of (204) is small due to the small value of ^^"("710) (i-e. this meson is mainly a ss excitation), while the 
second one nearly vanishes due to the fact that to/o(i7io) ~ 2M^. Note that, in the chiral limit, and 
for a vanishing mixing between the scalar isoscalars, the model predicts that the decay widths vanish 
exactly. For the qq states, this can be traced back to the quark loop which gives a phase-space like 
factoif^ irnf — 4(M° j,)^). For the glueball state, it is clear that the absence of coupling implies g* = 
(no quark content), so that the glueball has also no decay width. This vanishing of the decay widths 
of /o(1370) and /o(1710) in the pure NJL model (no mixing) in the chiral limit is a co nsequence of the 
chiral symmetry in the scalar sector. This has already been shown by Bajc et al. [167|. It means that, 
even if the meson to quarks coupling contant is large, the decay width into two photons can be small. 



Note that this phase-space factor can be negative. In that case, the coupling constant is also negative. The 
products (7j*"" X (m; — 4M^ ) and g*"" x (m^ — 4Mg ) are however always positive: the triangle diagrams with up 
quarks and strange quarks always add. Because of this change of sign of the coupling constant, we have plotted 
their modulus in figures BGI and ETl 
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We can conclude that this decay width may not be a good observable to claim whether a meson is glue 
rich or not. 

Figure also shows that the meson and glueball widths can va ry fr om one to three order of magnitude 



and that their relative magnitude is xo dependent (see fig. 7 of |165|) 



Our results are in complete disagreement with those of Close et al. |94| ], summarized in eq. (194) for 
scheme 1. The relative strengths are not reproduced, even qualitatively. In fact, we have, for xo — 350 
MeV, 

/o(1370) : /o(1500) : /o(1710) - 5.5 : 1 : 0.08. (205) 

Although being a glue rich particle, the /o(1500) has a two-photon decay width much larger than that 
of /o(1710), which is a quark rich particle. It arises from the fact that the latter has a mass almost 
corresponding to the chiral limit without coupling (to/o(i7io) ~ 2Afg), with small coupling to the up 
quarks (5*"" ~ 0), while the decay width of the glueball is not suppressed by the factors (^71-/^(1500) ~ 

29 . , ,n2. 



4M° , 771 (1500) ~ 4M° ). Moreover, as can be seen on figures 26 and g7[ the glueball to quarks coupling 
constant is not so small. 

One can wonder where is the origin of the discrepancy with the results |Q. Close et al. give very general 
arguments, while our results are probably dependent on the crudeness of the model. However, they show 
clearly that chiral symmetry effects {mf — 4(M°g)^) are important and that mass effects should also 
be considered (see eq. (6.2) of where the authors do not consider the it — s mass difference in the 
meson mixing; in the present model, these effects are considered). Should we have worked with the same 
hypothesis as in @ (M° = M° in the mixing matrix and in the chiral symmetry reduction factor^ | 
(mf - 4(M° J^), we would have obtained 

/o(1370) : /o(1500) : /o(1710) = 10.7 : 1 : 5.7, (206) 



in surprisingly good agreement with eq. ( |194| ). It seems that, working with the same hypothesis as |94[] , 
our model would be able to reproduce the same order of magnitude for the decay ratios. Because our 
results are very sensitive to the difference M° ^ M°, our model could then suggest that the M° — M° 
{SU{3) symmetry) hypothesis is a very drastic one and has to be rejected. We have however to keep in 
mind that this might be an artifact of the crudeness of our model. 

Finally, it is worth notic ing that the two-photon decay width of /o(1370) is always too small compared 
to the one given in |l04[| (5.4 ± 2.3 keV). Even restricting to the lower bound of the experimental value, 
one needs considerering low gluon condensates. This gives rise to a too small mass for /o(1370) (of the 
order of 1.2 GeV). 



Note that because the ao(1450) is not mixed to other channels, its mass is always very close to the chiral limit 
value 2M°. Due to the reduction factor, this implies that its width into 27 is very small, being exactly zero in 
the chiral limit. 
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IV. CONCLUSIONS 



We have reviewed the 3-fiavor A-scaled Nambu-Jona-Lasinio model both in the vacuum and at finite 
temperature and density. In section |, we have mainly described the features of the model and its 
connecti ons to QCD. Starting from a symmetry-anomaly point of view, we have used Fierz identities 
(section [ C ) to rewrite the local 4-point color-color interaction into the N JL form. For the sake of 
completeness, we have derived Fierz identities for a general covariant gauge, as well as for diquark modes. 
Although not needed in the present study, these generalized identities (compared to a Feynman-like 
gauge) are important to make contact with Global Color Models or Schwinger-Dyson equations, or any 
model which incorporates a covariant gauge. In particular, the Landau gauge is included in our Fierz 
identities. To reproduce the strong axial and scale anomalies, we have motivated the introduction of 
a new parameter (^) and of a dilaton field (x), respectively (section [D). After having bosonized the 
model, we have derived the gap equations, an equation for the dilaton field and the meson propagators. 
Apart from the drawbacks that we have indicated throughout the text, we have also included a specific 
section to discuss this important subject. 

In section U we have applied the model to study the thermodynamics of the 3-flavor A-scaled NJL 
model. We have investigated the quark condensates, the pressure and the energy and entropy densities. 
For vanishing densities, we have given both a low and a high temperature expansion. We have carefully 
defined a bag constant and have made warnings about the way to ext ract it from fits. We have also 
included a section to compare the NJL behavior with lattice QCD (section [I B 4 ) , as w ell as an application 
of the finite temperature formalism to the ttq — > 77 anomalous process (section II C). 
Finally, we have included a section on scalar mesons. The study, limited to zero temperature and densities, 
is interesting due to the coupling between the scalar isoscalar members of the qq nonet with the glueball. 
In order to obtain the physical states from a 3-particle coupled system, which is far from trivial because 
of the relativistic nature of the theory, we have indicated a way to extract masses and meson to quarks 
coupling constants. The procedure we have developed consists in studying the energy-dependent mass 
matrix. As an application, we have determined the two photon decay width of the scalar isoscalars. 
Although the model has to be taken with caution, especially in the scalar sector where qq mesons are 
unbound due to an unphysical threshold coming from the non-confining nature of the NJL type of models, 
we have shown that these decays might be not a good choice for the determination of the glue content of 
scalar particles. 
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APPENDIX A: CONVENTIONS 



The units h = c = ks = where h is the Planck constant normaUzed to 2tt, c the speed of Ught 
and fcs the Boltzmann's constant, are used throughout this review. 

In Minkowski metric, the metric tensor is 

10 

_ _ I -1 

5^1^-5 - I -1 

-1 

so that for = (t, x, y, z), we get x^ — g^^x'^ = (t, —x, —y, — z) and x'^ ~ x^x'^ = i^ — (x^+y^ + z^). 

• In Euchdean metric, we define x^ — — {it,x,y, z), which impUes x^ = x^x^ — x^x^ = —t^ + 

2 , 2 , 2 

X +y + z . 

• Dirac matrices are taken from |^^. In this reference, 75 = i7"7^7^7^. We define = (7°, 7) 
and we work in the Dirac representation, which means 7" = /3 = diag(l, 1,-1,-1). In Euchdean 
metric, we define 70 = 75 being unchanged (i.e. defined through /3 instead of 70). The relation 
{7^,7''} = 2gf"' from Minkowski metric becomes {7^,7''} = -26^"". 

• When not explicitly indicated, summation over repeated indices is used. 

• The totally antisymmetrical tensor g^'^"'' is defined by the relation: 

+1 {/i, ly, a, /?} even premutation of {0,1,2,3}, 
_ ^ 2 odd permutation, 
otherwise. 

Indices are lowered with the help of the metric tensor. Note that £0123 = — 1- 



APPENDIX B: COLOR QUARK CURRENT- CURRENT INTERACTION 

We start from the definition ( |l^ ) of the gluon propagator, for which the interaction term of the GCM 
model (Kw is q(x) (Ag /2)^^q{x)Duu (x, y) q(y) (A g /2)"f''q(y). We note that the spin part involves only the 



vector part of the F matrices (sse section ICS ): F" =7^, F^/ = j^. 

Making use of eq. (p5|), one can project 7^F'^ and 7,yr_a onto the basis F. The calculation can be 



performed usingf^ (together with the vanishing of the trace of an odd number of 7-matrices) 

tr(7'7^) = 0, (Bl) 

tr(7'^7'') = 45^^ (B2) 

tr(a^'') = 0, (B3) 

tr(7^7'^7'') = 0, (B4) 

tr(7^7''7''7'^) = Aig^'^gP'^ - g^Pg'"^ + g^^g'^P), (B5) 

tr(7^7^7"7''7") = -iie'"''"' , (B6) 

which allows one to obtain the following decomposition (for 7^ on F"): 



""when using eq. (|25|), we simplify the summation by noting that: 

6 

1 J2 rTtr(xr^) = \Y.T. '^^Mxan -llT. ^^^Mx^n 

T— 1 /.[ i/>/^ fi,iy 
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U 1/ 



(B7) 
(B8) 

(B9) 

(BIO) 

7/^^75 = -^757''. (Bll) 

With these equations, we are now able to evaluate (we leave aside the color and flavor d.o.f., reintroducing 
them when appropriate) 



1 



q{x)^^'q{x)D^,{x, yMyh^qiy) = -^q{x)rT^qiy)D^,ix, y)qiy)-,''TfMx), 



(B12) 



where the tensor summation is evaluated as described in the footnote (|9(]|). We are then left with the 
evaluation of 

-ii?^,(x,y)| (g(x)7''g(y)) [qiyh^qix]) + {qix)n'qiy)) {mYl^qi^)) 

+ {q{x)rr' q{y)) [qiyh^if.'qix)) - [qix)^^'^' qiv)) {qiy)Yi'if.'qix)) (B13) 

+ i(g-(.T)7^a^''''g(y)) (g-(y)7'^a,v.9(.T)) |. 
The first t wo terms are already on the F basis. The evaluation of the last three terms can be performed 



using eqs 



10|). Using the property D^^{x,y) — D^^{y,x), one can then show that 
iD^^(a;,2/)(g(a;)7''7^'g(y)) (qiyh^Jt^'qix)^ = 

-^D,4x,y)i^g'^''[qix)q{y)) {q{y)qix)) ~ [qix)a'^^' q{y)) [qiy)a\, q{x)) (g^^^ 

+2i{q{x)qiy)) {qiy)<y'^'' q{x))\ , 



\d^,{x, y) ( - g(a;)7^757'''9(2/)) [qiyh'l'^lt.'qi 

-\D^,{x,y)\^~g^^^(q{x)^'>q{y)) (q{y)l'q{x)) - (q{x)a^^' q{y)) (q{y)l\{x) 

-i (g-(x)a''"^"g(y)) (q{y)a^^"' q{x) 



"f^" (B15) 



, „ ..MM , ,„ ,,,1' 



and 



-il?^,(x,2/)i(q-(x)7^a^''''g(y))(g(y)7'^^MV'<z(: 

-iz)^.(x,y)| -g^''{qix)r'q{y)){qiyh^'q{x)) + {qix)j'' qiy)) {q{y)rq{x) 

+2z{q{x)r'q{y)^ [q{yh'j'''q{x))e,,^^''^, + ^{qixh^r" qiv)) {qiyh'^'" q{x) 



(B16) 



With t^^'P'^e^!' c " = -det(5"" ), (a; a') = (^, p, ct; p', a'), eq. ( pT5| ) gives 

- J^M-^Ca^' y) ( - 9(a^)7^7^7''''?(j;)) (9(y)7"7^7M'9(a^)) = 

-\D,,{x,y)\^~g^''(q{x)^\{y)) (q{y)l'' q{x)) - [qix)a^^' qiy)) {q{yh'q{x) 

+ l{qix)a,>.,q{y)) (q{y)a^''''q{x))gP'' - (q{x)a\,q{y)) (q{x)cj^^P' q{y) 



(B17) 
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whereas e'^'^'^'e^/ " = -2{gPP g"" - gf" gP applied to eq. ( |B16D gives 

-^D^,{x,y)l -g'"'[qix)r'qiy)){qiyht.'qix)) + (q{x)l''q{v))(q{y)rq{2 



g(a;)7^7^'g(y)) (^qiyH 1 f.' q^x)) g'"" + [q{x)l'' l"" q{y)) (q{y)l^l^q{^ 
+2i(qixhP'q{y)) (q{yWY' q{x)y,r 



(B18) 



Gathering these results gives the final answer for an arbitrary propagator D^y{x,y) within a current- 
current interationp^ 

q{x)-iPq{x)D^^{x,y)q{y)Yq{y) = -^D^,^{x,y) 



< js^"^ {q{x)q{y)) (q{y)q{x)) + g^^'^ (q{x)i^'' q{y)) [q{y)h'q{x) 

+{q{xh''qiy)){qiyh"qix)) + (q{x)i''q{y))(q{y)i''q{x) 



gP''(q{x)r'q{y)) {q{yh^.'q{x)) + (qixhS^'qiy)) {qiy^'l^qi^ 

+ {q{xh^i''q{y)) [qiyW^Y^qix)) - g'"'{q{xh^i'''q{y)) {q{yh''l^.'q{x) 

+2t{q{x)q{y)) [q{y)aP'' q{x)) + ^[q{x)a,,,,qiy)) (q{y)aP' q{x)) gP^ 
-{q{x)a\,q{y)) (q{y)aPP' q{x)) - (q{x)oPP' q{y)) (q{y)a\, q{x)) 

+2i{q{x)r' q{y)){q{y)l^l''' q{x))e,.P\, +i[q{x)aP 
Using the invariance of the propagator under fi v, we are left with 



(B19) 



To write down this Fierz-transformed action, we used the property D^^{x,y) = D^^{y,x). Without using 
tills property, one can find the Fierz transformation of the product (7'')rs(7'')tu, which is the generalization of 
the usual form 

{-/.)rs{intu = Wrumts + (^5 ( 475 ) + {i^)ruii^)ts + {i^)ru{i^)t. , 

and is given by (7^)rs(7'')tu = (7^r'')ru(7T/3)ts- This leads to 



4(7^)rs(7'')tu = gP^'SruSts + Qp" {il^)ru{i'y^)ts + (7^) r« (7" ) t s + (7" )™ (T** ) is - d'"' {l^ )™(7M')i3 

+ (7'7n-(7%nt. + (7'7'')-(7%'')t. - 5''^(7'7^')™(7'7M')t. 

+ 25^''('^'' " )ru(o-pv)ts ^ {^'^'^ )ru(o-''^/)ts - ((t"^/ ) ru (o"'^^ )ts 



For diquarks, we start from {"/,^)rs{'y'')tu = (7M)rs((7'')^)iit = {'yti.)ra{C'y'' C)ut = CuaCbt{'yP)rs{Y)ab, C = 17^/? 
being the charge conjugation matrix. This allows to write, using the previous Fierz transformation. 



4(7^)rs(7'')t" = g'"'CrtCus+gP''{il''C)rt{Ci'i^)us + {^P C)rt{C-f'') us + (7'' (^7" " <?^''(7'' C)rt{C-i^,)^ 



Hl^CUiCYYU + (7^"'C)..(C77^7^).. -P^''(7 7" CUiC-^^i^,). 



+i 



-(i'y-'C)rt{CaP ) 



+ {1" CUiCY-yP Ue^'.'P- - (77^ C)rt{C-f- Ue^,, 
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x< 5 



+2(qix)rqiy)) (qiyh^qi^)) ~ .9^"^ (g»7^'«(y)) {qivh^^'qi 

+2(g-(:r)7S^(y)j (g-(y)7''7''9(x)) - g'^" {q{xh'l^' q{y)) {q{yh'l^'q{x) 
+ i.9^'^(g»^M'."Z(2/)) (9(y)a^'^'g(x)) - 2{qix)a%,qiy)) (q{y)a^^^' q{x) 



(B20) 



With the Fierz transformations for color and flavor (see section [ C 1 and [ C 2 ), we have - taking then 
into account the appearance of diquarks- 

q{x)^ Yq{x)D^,^{x,y)q{y)^ q{y) = -~^D f,^{x,y)^ 



- 2 
-4 
+ 



^ (q{x)G'^K"^q{y)) (q{y)G^Kq{. 

a,e 

[q{x)GH^y{y)) (q{y)GHYq{x)) ~ 2(q{x)GH-i^rq{y)) {q{y)GH-i''Y q{x) 
(g-(x)G^-^g(y)) (q[y)G^ — q{x) 



-2( 
-2( 



^Y.[q{x)G\/-^K'^Cf{y)) (f{y)CG\/-^Kq{ 
(q{x)G%ir'-^Cf{y)) (f {y)C G%'-^i^'^ q{x)) 



[q(x)G'X'^t^'rCq^{y)) (q" {y)CG\'-^i^''^^ q{x) 



V2 



4 (q{x)G^s '-^^Cf {y)) {q^{y)CG% '-^^q{x) 



(B21) 



where we have defined the shortened notation^ K'"' = (1,17^,17'^ ,^7^7'^ ,0''^ /v^)- Note that a cofor 
triplet diquark state is antisymmetric in the exchange of colors; then - this is the Pauli principle - it must 
be antisymmetric under the exchange of the other coordinates (both spatial and internal) . This implies 
that G')^ has to be associated to (1, ij^, 17^7'' ) while Gg has the Dirac counterpart if '° = (27^ , cr^ / 
Mathematically, these relations can also be shown to hold in the path integral formalism, where the 
quark fields are anticommuting variables. This implies that (see eq. (Elh for the explanation of this 
term) (fCV^q — —q^{CT°'Y'q which, in turn, implies CT" = -~{GT^^ . Because the color part is 
antisymmetric, it is clear that the expected relation shows up: {GK'°'G\ 5)^ — {GK'°'G\ g). 



APPENDIX C: PROPAGATOR MATRICES 



The propagator matrices indicated in eqs. (^,^) are defined below for symmetric matter (/i„ = /id 7^ /^s 
in the isospin limit (m„ = 7^ m^, Afi = fi^ — /ig). 



P^^=[q^Z^{P,^,)+a\l'-^ 



(CI) 



^The factor l/\/2 for the tensor part comes from the convention defined in the footnote (po|). 
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APPENDIX D: THERMODYNAMICS 



1. Pressure 



For a vanishing chemical potential, eq. (129) reads 



Pi 



ideal , 



Nc 
) 



oo ,4 



1 



Eul + ePE, 



-dk 



1 



1 + ePi^' 



-dk 



Only one of these terms has to be analyzed. We define 



2Nr 



3/2 



dy, 



(Dl) 



(D2) 



where the infinitesimal quantity e will allow to regularize the summation that will be encountered in the 
following. Expanding, we get 



Pi 



ideal gas 



is) 



2N, 
3^ 



M.^lun^C-ir^e— / (y^-i; 

^ n=l 



(D3) 



which can be expressed in terms of the modified Bessel function of order two [168| 
Pidcalgas(s) = — hm^(-l) e 



(D4) 



This series has also been investigated in |62 . Because of the fast decrease due to the factor an d 

because of the asymptotic behavior of K2, it can be numerically more advantageous to use (D4) than (D3). 



a. High temperature zero density expansion 



The converging factor e in (D3) and (D4) is necessary to obtain nondiverging quantities in t he high 
temperature expansion. In such an expansion, only the first two terms are finite. Following |168|, 



K2 (z) = 2z- 



1 °° 



■ip{k + l) + ^j{k + 3) fz 



k=0 



k\ (2 + fc)! 



-t) -^^7;h{z)., 



(D5) 



where tp is the Digamma function (defined as d\nT{z)/dz with T{z) the Euler Gamma function) [168|, 

n ^ 

V^(n+1) = + with ^(l) = -7, (D6) 



fc=i 



and where l2{z) is the other modified Bessel function of order two |16S] 



/c=0 



fc! (fc + 2)! 



(D7) 



In (D6), 7 is the Euler constant. 

The logarithm ln(z) in K2{z) implies a singularity at the origin. There is no Taylor expansion arou nd it . 
This problem, and the way to circumvent it through the converging factor, has been established in |ll9|] 
for a fermionic gas (only the first few terms of the expansion are given) and in |l2l| |123| for a bosonic 
case where a full expa nsio n, valid also at non zero density, has been given. 
Combining (p5|) and (D7), we obtain 



4 -T - ? E 



fc=0 



A:!(fc + 2)! 



^(V'(fc + l) + V(A: + 3))-ln(0 



(D8) 
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This leads to 



/3)2 \{nMsPY 2 



4^^^ 1)"^'^^ "'fc!(fc + 2)! ( 

n=lfc=0 \ ' ) \ 



1 ({nMslif 



i(^(fc + l) + ^(fc + 3))-ln('i^^ 



(D9) 



The first two terms are easy to determine and coincide with the two non-diverging terms from the Taylor 
expansion around Ms — 0. They give 



27V, 

TT 



°° 11 



2 ' 1 M4/34 ^^^^ ^ 2Af4/52 



n=l J 



(DIO) 



They are tabulated in |168| and can be expressed through the use of the Riemann zeta function^ 



^ oo 

11=1 

C(-2rn) = 0, m=l,2,..., 



(_!)«+! 



Re(z) > 



C(2m) 
which leads to 



(2m!) 



-, with 771=1,2,... and B the Bernoulli numbers. 



180 12 
The contribution of the term fc = is rewritten in the form 



47r2 
Using 



'""^Mtl hmi 



f;(-l)"+ie— (1(^(1) + 7^(3)) - ln(^)) - f;(-l)"+ie— Ini 

n— 1 n— 1 



(Dll) 
(D12) 
(D13) 

(D14) 
(D15) 



lim X](-l)"+'e— 77^ = (1 - 2i+^)C(-z), 



(D16) 



we have 



— V(-l)"+ie-""77^ = V(-l)"+ie-"^77^ ln77 = — ((1 - 2^+')a-z)) 

dz ^ — ' — ' rl-y ^ 



dz 



(D17) 



which leads to 



^(-l)"+ie-"=77^ ln7^ = -(ln2)2i+X(-z) - (1 - 2^+^)C'{-z), 



(D18) 



where C'(~-z) = -rC{^) 
dz 

C(0) is obtained from 



1^ p-ne 1 

V(-l)"+ie-"'' = = — when e ^ 0, 

n— 1 



(D19) 



^For 'Re{z) < which we shall need in the following, a converging factor e "'^ is needed. 
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while C'(0) is given in [168 



C'(0) = --ln2^. 



With (Die) and (D2C), (D15) gives 



(D20) 



(D21) 



We still need the A: > terms of the expansion. With (D16) and (D18), they are 



lim < > 

£ — ^ 



2k 



Eq. ( |D6D implies 



fc!(fc + 2)!4'= 
+ (ln2)2i+2'=C(-2fc) + (1 - 2i+2'=)C'(-2fc) 

V'(l) = -7, 
V'(3) = -7+| 



1 (V(fc + 1) + V'(fc + 3) - 2 In(^) ) (1 - 2i+2'=)C(-2fc) 



(D22) 



(D23) 



so that ( ^T^JD14^2i1 , |D22|) lead to 



7 iV 

Pideal gas(s) = ^^^c^^T^ - ^M^T^ 



(D24) 



It is clear that the converging factor has regularized the summationf^ of the expansion. Note from (D24) 
that the separation into a logarithmic term and a constant one is arbitrary: one can always write ln(a6/3) = 
ln(a/?) + ln6 (a and b being dimensionless constants) and put In 6 into the constant term. This has some 
importance for the interpretation of the high temperature results. 

To write (D24) into a form involving only elementary functions, we still need to know C,'{—2k){k > 1), 
i.e. 



1 



(1 - 2i+2fc) ^ 

^ ' 71=1 



^(-l)"+ie-"^n2'= Inn, 



because of (PIBD and (|D12D . Using |68| 



21-zj. 



(z)C(z)cos(^) =7r^C(l-^) 



and 



we have 



r(z)r(i-z) 



sin(7rz) ' 

2--^CWcos(f)=.^r(i-.)C(i-.), 



(D25) 



(D26) 



(D27) 



(D28) 



Summations of this kind are called Euler sums [ 16£ ] 
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so that 



With 



lini . , , az) = {-\f2-'^^-^-K-'^^V{\ + 2fc)C(l + 2k). (D29) 
z^-ik sin(7rz) 



li^ ^L^CW = = ^qZ^^C'(-2fc), (D30) 
,^_2fesin(7rz)^^ ^ 7rcos(-27rfc) ^ ^ ^ ' 



we finaUy obtain 

1 



C'(~2fc) = ^(-l)'=(27r)-2'=r(l + 2/c)C(l + 2k). (D31) 



Eq. (D24) is then rewritten into the form 



7 

Pidcal gas(s) = J^N^^^T^ - J^M^T^ 

+ E - 2^+^^-)^(-l)'^(2-)-^'=r(l + 2fc)C(l + 2k), 



which only necessitates the evaluation of known functions. 



(D32) 



b. Low temperature zero density expansion 



We can search for a low temperature expansion, /3 oo, starting from (Dl) or (D4). The last one is 
better suited because of the well known asymptotic expansion of K2 |l7Ct| 



K,{z)^^ -e-^-\l + — 

" 2z I 8z 



- 1)(422 - 9) - l)(4i2 _ 9)(422 - 25) 



2!(8z)2 



Combining (g) and ( p3^ ), we have 

2A^, 



Pi, 



deal J 



1 + 



15 



3!(8z)3 



(D33) 



15-7 



15- 7- (-9) 

8nA/s/3 ^ 2!(8nMs/3)2 3!(8nMs/?)''' 



(D34) 



When /? is large, we are in the chirally broken phase where the quark masses are the constituent masses. 
Since the expanding para meter is PMi, i — u,s, the approximation (D34) becomes better as /3Mi is 
increased. In section |II B| , it is shown that the mass variation is low for T ^ 100 MeV. For a constituent 
quark mass of about 400 MeV (at T = 0) (3Mi is, at least, 4. The expansion (D34) is then perfectly 
justified. In that case, the first term n = 1 is enough and the pressure is 



-Pidcal gas {^) 



(27r)3/2 

15 15-7 15 •7- (-9) 

8Ms(3 ^ 2!(8Ms/?)2 ^ 3!(8Ms/3)3 



(D35) 



We have checked that for the set of parameters (Af° = 300, xo = 80) MeV (see section 11 B ) this expansion 
is not well suited. In that case, the second term n = 2 in ( D34D , as well as the three corrections to "1" 
for both n = 1 and n ~ 2, are necessary to reproduce results valid up to 100 MeV. 
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c. Finite density, zero temperature 



For finite density at vanishing temperature, eq. (129) can be exactly integrated. We have ni+ = 0, Ui 
— 9{fii — Ei) and /i^ = \ kjp, + Alf, where kp. is the Fermi momentum of the ith flavor, so that 



P. t \ 



37r2 



kp Us 



^F, 3 



■In 



fee 



M, 



(D36) 



2. Energy density 

High temperature zero density expansion 



Combining eqs. (137, D32) immediatly gives 



7 N 

Eidcal gas(s) = ^^c^T^T* - j^M^T^ 



167r2 M ' 2 



TT 



N 

0^2 ""s 



fc=l 



^fff^''^j 2fc + 1)(1 - 2i+^^-)i(-l)'=(27r)-2'=r(l + 2A:)C(1 + 2k). 



(D37) 



b. Low temperature zero density expansion 



Eq. (137) is not well-suited because it would imply taking the derivative of a truncated series (see 
eq. (D34)). It is bett er to search for the expansion of the exact solution in terms of the modified Bessel 
function obtained in |115 



r,^-3P (,)^^^^M^\-( iV'+iM^ 



(D38) 



Its low temperature asymptotic expansion is obtained from eq. (D34) for the flrst term and from eq. (D33) 
with i — 1 [170 1 for the second one. We then have 



^ideal gas (^) 

X <! 1 + 



15 



15-7 15 •7- (-9) 



2N, 



8nMs(3 2!(8nA'f,/3)2 3!(8nM,/?)3 

cx: ^ 



X n + 



nMsf3 V 2nM,/3 

3 -(-5) , 3 -(-5) -(-21) 
3!(8nM,^)3 



(D39) 



Once again we can limit ourselves ton = lorn = l,2, depending upon the chosen set of parameters 
(M°,xo)- However the first two or three corrections to "1" are necessary. 



Finite density, zero temperature 



As for the pressure, eq. ( 130 ) can be exactly integra ted. It is however more judicious to use the vanishing 
nature of the entropy at T = in order to get (eq. (128)), 
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^idcal gas ^ -^idcal gas ~^ f^iPi^ 



where pi is given by eq. (|122| ), i.e 



Using ([D36D , (|D40|) and (|D4l|) , we obtain 



Nr. 



3. Entropy density 

a. temperature zero density expansion 



Using eqs. ( |l40| , p3^ ) gives 



- S^'^^fi WTW^^'^' 2i+-)l(-l)'^(2.)--r(l + 2A:)C(1 + 2k). 



b. Low temperature zero density expansion 



We can obtain this expansion starting from eq. (128) with /i^ — 0: 

s(s) = /3(Pidcal gas(s) + Eidcal gas(s)), 



SO that, using eqs. ( D34 , D3£ ), we have 



15-7 15- 7- (-9) 



15 



2iV, 



8nMs/3 2!(87iM,/3)2 3!(8nM,^)3 

n+l ^ 



^M4/3^(-1) 

n=l 

3 



X n + 



3 -(-5) , 3 -(-5) -(-21) 



8nAf,/3 2!(87iAf,/3)2 3!(8nM,/3)3 

Once again we can limit ourselves ton = 1 orn = 1,2, depending upon the chosen set 
(M°,xo). However the first two or three corrections to "1" are necessary. 



c. Finite density, zero temperature 

It is clear that 

s{s) = 0, 

in agreement with the third principle of thermodynamics. 
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APPENDIX E: MASSES AND ON-SHELL COUPLING CONSTANTS FOR MIXED 

PARTICLES 



In this section, we define the masses and on-shell couphng constants for mixed particles. In appendix 
the scalar isoscalars form a set of three coupled particles, whose inverse propagator is given by the 3x3 
subset of the matrix T. Searching for the physical particles means that this matrix has to be diagonalized. 
In the following, we show how these physical states can be obtained. This is a non-trivial procedure 
because of the energy dependence of the normalization functions Z . The following demonstration is valid 
for a set of n coupled particles. 

In order to shorten the notation, we use the following convention throughout this appendix: we denote 
by A{q^ = —mf) a matrix A of which each element is evaluated at the same in opposition to A{—mf), 
which indicates a matrix whose first column is evaluated at —ml, the second at —m\, the nth at 



1. Generalities 



Let r be the n x n meson matrix propagator 

ru^A,q^ + B,,rj,^r,j. (ei) 

It is better to first consider Ai, Bi, ^ i) as being independent of q^ . In this case, an exact diagonal- 
ization (valid for each g^) can be performed. We proceed in the following way: 

a^Ta = a^g-^V'g-^a, (E2) 

where (a is the vector denoting the coupled fields) 

r-. = 9' + x,r^ = and g-' = diag(V^, v^). (E3) 

A.i yJAiAj 

Writing T' = g^l + T", which defines T", we deduce 

det(r' - A'l) = det(r" - (A' - q^)l) = dct(r" - A"l). (E4) 

Since F" is (/^-independent, so are the eigenvalues (e.v.) A". Writing A'/ = m^, it is easy to see that the 
e.v. of r' are 

\[^q^ + ml (E5) 

Moreover, it is clear that the eigenvector (E.V.) matrix associated with F' is identical to that of F". This 
E.V. matrix is denoted by V . It is g^-independent (it can then be evaluated exactly) and is orthogonal: 
= V^ . In short: 

cr^FcT = CT^y- Vdiag((j2 + m^, q^ + ml)V-^g-^a, (E6) 

i.e. 

a^Fcr = I]^diag(g2 + ml, q^ + m^)!], (E7) 

with 

S = V-^g-^a, (E8) 

the vector denoting the physical fields. 

On the other hand, a direct diagonalization implies 

a^Ta ^ a^V,{q^)Amg{\i, ...,X^)V-\q^)a. (E9) 

At this level, everything is exact (because this is the same q^ for each E.V.), and we have the property 

V-\q') = V^q'). 

Writing 
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A(g2) = y-i(q2)rv;(q2) =diag(Ai, A„), (ElO) 

and using 

±^{V-\q^)VM^))^Q, (Ell) 



we have 



^[^^^^^V4q%A{q')] + V-\q^)g-'V4q'). (E12) 

The physical masses are clearly obtained from the zeroes of the e.v. Since the trace of a commutator is 
identically zero for a finite matrix, we get the following relation between the exact quantities^ 

tr5*"'(<z') = trg-^. (E13) 

Another exact relation comes from 

a^V,{q^)A{q^)V-Hq^)<y = g-'Vdiagiq^ + ml q'' + ml)V-' g-'a, (E14) 

which implies 

det(K(g2)A(g2)^-i(^2)) ^ det(g- Vdiag((z2 + ml + ml)V-'g-'), (E15) 

i.e. 

det(A(g2)) = det(g" Vdiag((?^ + ml q^ + ml)V^^) = det{g^^) det(diag(g^ + ml q^ + ml)), 

(E16) 

i.e. 

Let us once more emphasize that this relation is only valid as far as the e.v. are all evaluated at the same 

2. First order Taylor expansion, as a function of , around the respective zeroes of the 

eigenvalues 



If each diagonal component of eq. ( E12| ) is evaluated at the zero of the corresponding e.v., the commutator 



gives a zero contribution. The latter can be established from 

g.-^(g^ - -raj) = [ ^^- y V .(g^), A(g^)] +V-\q'^~m^)g-^VAl'--m!). (E18) 

Since each quantity is here evaluated at the same q^, the inverse E.V. matrix is identically equal to its 
transposed. Taking out the j,j component on both sides of the equality, we obtain (summation over k 
and I, no summation over i) 



^^We call "exact" the quantities which are all evaluated at the same q^; approximations only play a role when 
quantities must be evaluated at their respective mass shell. 
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jk 



A{q 



2^ f dv-Hq '),.,^.2^ 

\ dq^ 



-V.{q^) 



2\ — 2t/ / 2 2\ 



(E19) 



Since A(g^) is diagonal, it is clear that the contribution of the commutator vanishes. Anyway, if we take 
j = i, this commutator gives also trivially zero since it implies the evaluation of Xi{q^ = —mf) which is 
zero by definition. Now, let us effectively take j = i. By definition of the transposed of a matrix, we then 
have the relation 



E 

k 



9kk = ~mj)k^y 



(E20) 



which sho ws tha t the zth E.V. plays a role only for the component z, i of ^(q^ = —mf). Then, we can 
generalize ( E20D to 



(E21) 



where, according to the conventions of this appendix, each column of V and ^ is now evaluated at its 
respective zero. Although we have demonstrated eq. (E21) for diagonal components, it can be extented 
to 



1„2 



i-mf) = g:'V-'g 



(E22) 



It includes the previous one for the diagonal components of g^^. But eq. ( E22| ) can also be shown to 
imply the vanishing of the off-diagonal elements. This is seen as soon as we can prove that 



V-\-m^) ^ g.{-m^)V-'g-' 



(E23) 



i s tru e, with being g^^ the diagonal matrix taken from the defin ition dia.g{d\i{q^)/dq^). The relation 
(E23) is suggested by a comparison between eqs. (E7) and (E9). Let us notice that it is the inverse 
V^''{—mf) which mu st be considered, not V'^ {—mf)(F,.Y . eva luate d at different q^). 
In order to prove eq. (E23), we start from the E.V. definition (ElO): 

r{q^)K{q^)^K{q')A{q'). 

Taking the component i, j, evaluated at q^ = ^fn'^, we obtain 



(E24) 



T{q^ 



-'m^)ikV*{q^ = -m^Akj = T4(<7^ = -m^A^kHl'^ = 



i.e., (no summation over j) 

T{q^ = -m^kV^iq^ 



-m^j)k] = V^{q^ 



(E25) 



(E26) 



It is clear that the on ly E .V. which plays a role is associated with the jth e.v. Moreover, by definition of 
rn^, the r.h.s. of eq. (E26) is identically zero. We then have: 



T{q^ 



(E27) 



It is useful to stress that eq. ( E27 ) has to be verified component by component since the q used for T 

(i) 

is t he on e of the corresponding E.V. in K. If we denote by rf* the E.V. associated with the jth e.v., 
eq. (E27) is then identical to 



T{q' 



-mj}.,kd^. 



^ T{q^ 



0. 



(E28) 



The demonstration of eq. (E23) consists in starting from the relation (E27) and in showing that ( |E23D is a 
possible solution thereof. It is then sufficient to show that the normalization is correct for the uniqueness 
of the solution. (Note : the last relation does not imply the existence of n orthogonal E.V. associated 
with q^ = —m?j- There i s ind eed only one E.V. associated with the e.v. A(— m|) = 0.) 
If wc take the inverse of (E23), i.e. 



V*{-m^)=gVg:\ 



(E29) 
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and if we insert it on the l.h.s of (E27), we have 



i.e. (see eq. (|E2D) 

r(g2 = ~7n%k{9Vg:%, = {g-^V {^m^^)V)ag:}y 
By definition of we have 

T'V ^Vdmg{q^ + 7vl...,q^ ^ml), 
so that (there is no summation over I) 

{g-^V'V)u = {g-^V),kAi^g[q^ +ml...,q^ +ml)ki = {g-'V)uiq^ + mf). 
Then, inserting ( p33| ) in ( |E31| ), we get 

r(q2 ^ -m^^MgVg;\, - + 

and then, because g~^ is diagonal, 



(E30) 
(E31) 
(E32) 
(E33) 

(E34) 

(E35) 



We have thus shown that (E2£) is a solut ion o f (E27). In fact , ( |E29| ) is exact since the normalization is 
correct (because, from ( E29| ), we deduce (E22) and then (E21) which was shown to be true). 
As a conclusion to this appendix, we note that with the E.V. matrix K(— m^), the physical fields are 
given by 



S = V-'g-'a = g-'i-mi)V-'i-mf)a. 



(E36) 

This definition will be kept even if the coefficients Ai, see eq. (|E1|), are q^-dependent (this is for example 
the case of the NJL model). Nevertheless, the previous relat ions are not all valid since t he m atrix g^^ 
becomes also g^-dependent. However, we can conclude from ( |E19[ ) that the equivalent of ( E20 ) is: 

9:Hq' - -m^h -Y.alu{<l' = W - -^lU" -Y.9kk{<l' = -™?) (K(-m,^),.)' , (E37) 

k k 

as far as the coeffic ients Bi and the elements Tij, i ^ j are g^-independent. In case they are q^-dependent, 



(E38) 



it is obvious from ( |E19| ) that 

9^\q^ = -m,2),, = Vf (q2 = ^r7?,Ugkhq^ = -^i^WM^ = -m^)u 



k,l 
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